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Abstract 



Probes of the Early Universe 
by 

Max Erik Tegmark 
Doctor of Philosophy in Physics 
University of California at Berkeley 
Professor Joseph Silk, Chair 



One of the main challenges in cosmology is to quantify how small density fluctuations 
at the recombination epoch z ^ 10^ evolved into the galaxies and the large-scale structure 
we observe in the universe today. This thesis discusses ways of probing the intermediate 
epoch, focusing on the thermal history. The main emphasis is on the role played by non- 
linear feedback, where a small fraction of matter forming luminous objects can inject enough 
energy into the inter-galactic medium (IGM) to radically alter subsequent events. The main 
conclusions are: 

• Early structures corresponding to rare Gaussian peaks in a cold dark matter (CDM) 
model can photoionize the IGM early enough to appreciably smooth out fluctuations 
in the cosmic microwave background radiation (CBR) , provided that these early struc- 
tures are quite small, no more massive than about 10^ Mq. 

• Typical parameter values predict that reionization occurs around z = 50, thereby 
reducing fluctuations on degree scales while leaving the larger angular scales probed 
by COBE relatively unaffected. 

• For non-standard CDM incorporating mixed dark matter, vacuum density, a tilted 
primordial power spectrum or decaying r neutrinos, early reionization is not likely to 
play a significant role. 

• For CDM models with CI < 1 and A = 0, the extent of this suppression is quite 
insensitive to CIq, as opposing effects partially cancel. 

• It is still not certain that the universe underwent a neutral phase, despite the new 
COBE FIRAS limit ?/<2.5xl0~^on Compton ^-distortions of the cosmic microwave 
background. 

• The observed absence of a Gunn- Peterson trough in the spectra of high-rcdshift quasars 
can be explained without photoionization, in in a scenario in which supernova-driven 
winds from early galaxies reionize the IGM by 2; = 5. 
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• It is possible to place constraints on cosmological models that are independent of 
the shape of the primordial power spectrum — the only assumption being that the 
random fields are Gaussian. As an example of an application, the recent measurement 
of bulk flows of galaxies by Lauer and Postman is shown to be inconsistent with the 
CBR experiment SP91 at a 95% confidence level regardless of the shape of the power 
spectrum. 
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Preface 



As this thesis is laden with cosmology and astrophysics jargon, it is recommended that 
the non-cosmologist reader flick through Chapter 2 before attempting to read anything else. 
This chapter is by no means intended to be a comprehensive treatise on the foundations of 
cosmology. It is merely an attempt to provide the basic knowledge necessary for a reader 
with a general physics background to be able to follow the other chapters. Cosmologists are 
advised to skip directly from the introduction to Chapter 3. 

During my graduate studies here at Berkeley, I also did some non-cosmology work 
(Tegmark 1993, Tegmark & Yeh 1994, Tegmark & Shapiro 1994, Shapiro & Tegmark 1994), 
but that research will not be discussed in this thesis. 
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Chapter 1 

Introduction 



One of the main challenges in modern cosmology is to quantify how small density fluctu- 
ations at the recombination epoch, about half a million years after the Big Bang, evolved 
into the galaxies and the large-scale structure we observe in the universe today, some 10^*^ 
years later. This thesis focuses on ways of probing the interesting intermediate epoch. The 
main emphasis is on the role played by non-linear feedback, where a small fraction of mat- 
ter forming luminous objects such as stars or QSO's can inject enough energy into their 
surroundings to radically alter subsequent events. Specific questions addressed include the 
following: 

• Which cosmological models predict that the universe will become photoionized early 
enough to measurably affect spatial fluctuations in the cosmic microwave background 
radiation (CBR)? 

• Which cosmological models are ruled out by the constraints from the COBE FIRAS 
experiment on the spectral distortions of the CBR? 

• Which cosmological models satisfy the Gunn- Peterson test, i.e. predict the intergalac- 
tic medium (IGM) to become highly ionized by z = 5? 

• Is reionization required to make bulk flows consistent with degree-scale anisotropics? 

Thus this thesis provides a quantitative link between cosmological models and extragalactic 
observational astronomy, showing how current observations constrain the allowed parameter 
space of cosmological models. 

The next two sections of this introduction are aimed at placing the work in a larger 
context, and the final section gives a brief description of how the thesis is organized into 
chapters. As mentioned in the preface, the non-cosmologist is advised to flick through 
Chapter 2 first, before attempting to read further. 

1.1 Cosmology towards a New Millennium 

Like the universe, the field of cosmology is expanding. Thanks to the rapid technological 
developments of the last few decades, wc humans now have better eyes than ever before with 
which to scrutinize our surroundings, both from the ground and from various spacecraft. 
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On the ground, the recently inaugurated ten-meter Keck telescope on Hawaii provides 
record-breaking light-gathering power in optical wavelengths. Interferometry between radio 
telescopes on separate continents provides angular resolutions that our ancestors could only 
have dreamed of. Semiconductor technology has replaced photographic plates by much more 
sensitive electronic detectors (CCDs), and also, in the form of computers, enabled previously 
unheard of capacity for data processing. A notable example is the MACHO project, looking 
for dark matter gravitational microlensing events in the halo of our galaxy, a search that 
involves looking at some one million stars every night (Alcock et al. 1993). Another example 
is the Sloan Digital Sky Survey, whose goal is to measure the redshifts of a million galaxies 
before the end of the millennium. To place this in perspective, when Sliphcr obtained the 
the 18 data points that Hubble used to support his famous expansion law (Hubble 1929), he 
required weeks of observations with a by contemporary standards good telescope in order 
to obtain the redshift of a single nearby galaxy. 

In space, scores of detectors have been observing the cosmos unhindered by the as- 
tronomer's worst enemy: the earth's atmosphere. After being repaired, the Hubble Space 
Telescope now obtains optical images with ten times better angular resolution than obtain- 
able from the ground. The Cosmic Background Explorer satellite (COBE) has made an 
all-sky map at millimeter wavelengths (Smoot et al. 1992). The IRAS and ROSAT satel- 
lites have done the same in the infrared and x-ray wavelength bands, respectively, and a 
number of satellites have explored gamma rays and the ultraviolet as well. Many of these 
wavelengths are almost completely absorbed by the atmosphere, and thus inaccessible to 
ground-based experiments. An overview of our knowledge of the cosmological photon spec- 
trum in all wavelength bands from 10~^^cm to 10^ cm is given by Ressell & Turner (1990). 

With all this data pouring in, it is a daunting task for theorists to keep up and try 
to make sense of it all. Although it cannot be overemphasized how much remains to be 
done, the advances (interrupted by occasional false clues) of the last few decades have been 
quite astonishing by any standards. As a more detailed example, let us take something 
that is directly relevant to the material in this thesis: the cosmic microwave background 
radiation (CBR). Thirty years ago, we did not even know that it existed. Confirming 
theoretical predictions of Gamow and others, it was serendipitously discovered by Penzias 
and Wilson in 1964 (Penzias & Wilson 1965) and given its cosmological interpretation by 
Dicke, Peebles, Roll & Wilkinson (1965). In Figure 1.1, the best limit on the Compton 
y-parameter is plotted as a function of year. This parameter is a measure of how much 
the spectrum deviates from a Planck spectrum, and when it was first defined (Zel'dovich k, 
Sunyaev 1969), one was still arguing as to whether the CBR had a Planck spectrum or not, 
so the constraints on y were of order unity. About a decade later, the limit had dropped to 
about 0.25 (Woody & Richards 1981). Almost another decade later, it dropped to around 
0.1 thanks to a Berkeley-Nagoya rocket experiment (Matsumoto et al. 1988). Here a fairly 
common cosmo-sociological effect occurred: an apparent deviation from a Planck spectrum, 
now sarcastically referred to as "the Berkeley Bump", led to scores of theoretical models 
explaining it before it was officially admitted to be merely an artifact caused by problems 
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with the detector. In 1990, the FIRAS experiment aboard the COBE sateUite suddenly 
pushed the hmit down by a full two orders of magnitude (Mather et al. 1990). And the 
exponential progress continued, the current limit being 2.5 x 10~^ (Mather et al. 1994). 

This mind-boggling progress with CBR spectral distortions has been paralleled by im- 
proved measurements of the spatial distortions. The fairly isotropic radiation field was 
observed to have a dipole anisotropy of order 10^'^ (Conklin 1969; Henry 1971; Smoot et 
al. 1977), which is normally interpreted as a Doppler shift due to our motion relative to 
the CBR rest frame. A few years ago, the COBE DMR experiment detected smaller scale 
fluctuations as well (Smoot et al. 1992), at a level of about 10~^, which are believed to 
be of cosmological origin. Since this time period, there has been an enormous surge in 
balloon-based and ground-based CBR experiments that probe smaller angular scales, and 
now in the spring of 1994, new results are coming in at a rate of almost one per month. 

It is precisely this type of data, in combination with improved constraints on the matter 
power spectrum from galaxy surveys, that is allowing us to start pinning down the thermal 
history of the IGM since recombination, the main topic of this thesis. 

1.2 Reheating 

It is now widely accepted that the universe underwent a reheating phase at some point after 
the standard recombination epoch at redshift z 10^. The absence of a Gunn- Peterson 
trough in the spectra of high-redshift quasars has provided strong evidence for the reheating 
occurring at a redshift z > 5, since it indicates that the intergalactic medium (IGM) was 
highly ionized at lower redshifts (Gunn &: Peterson 1965; Steidel & Sargent 1987; Webb et al. 
1992). The smallest baryonic objects to go non-linear in a standard cold dark matter (CDM) 
model are expected to reionize the IGM at a redshift somewhere in the range 20 < z < 100 
(Couchman & Rees 1986). In the most recent models with baryonic dark matter, reheating 
and reionization are predicted to occur at an even higher redshift, typically in the range 
100 < z < 1000 (Peebles 1987; Gnedin k Ostriker 1992; Gen et al. 1992; Cen et al. 1993). 

A reheating epoch would have at least two interesting classes of effects that may be 
measurable today: effects on subsequent structure formation and effects on the cosmic 
microwave background radiation (GBR). 

Subsequent structure formation would be affected in a number of ways. First of all, the 
heating of the IGM up to a higher adiabat would raise the Jeans mass, thus suppressing 
the formation of small objects. For instance, an IGM temperature of lO^K at a redshift of 
a few would suppress the formation of galaxies of mass below lO^^M©, thus alleviating the 
ubiquitous problem of theories overpredicting the abundance of faint galaxies (Blanchard et 
al. 1992). Secondly, if the objects that reheat the IGM also enrich it with heavy elements, the 
ability of gas to cool would be greatly enhanced in the temperature range lO^K <T < lO^K, 
presumably facilitating future structure formation. 

The GBR would be affected in at least three ways: 

1. Hot ionized IGM would cause spectral distortions (Kompaneets 1957; Zel'dovich &; 
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Sunyaev 1969; Bartlett & Stebbins 1991) which might violate the stringent hmits on 
the the Compton y-parameter (Mather et al. 1994). 

2. Spatial fluctuations on angular scales below a few degrees might be suppressed, while 
fluctuations on larger scales would remain fairly unaffected. This effect is particularly 
interesting since some degree-scale microwave background experiments carried out at 
the South Pole (Meinhold &; Lubin 1991; Shuster et al. 1993) and at balloon altitudes 
(Devlin et al. 1992; Meinhold et al. 1993; Shuster et al. 1993) appear to detect 
fluctuation amplitudes lower than those expected in a CDM model normalized to the 
the COBE DMR experiment (Smoot et al. 1992), which probes larger angular scales. 
The interpretation of this remains unclear, however, as many other experiments {e.g. 
de Bernardis et al. 1993; Dragovan et al. 1993; Cheng et al. 1993; Gundersen et al. 
1993) have detected higher levels of fluctuations. 

3. New spatial fluctuations will be generated on smaller angular scales, through the so 
called Vishniac effect (Vishniac 1987, Hu, Scott & Silk 1994). The current upper 
limit on CBR fluctuations on arcminute scales {e.g. Subrahmanyan et al. 1993) places 
constraints on some reheating scenarios. 

With the recent surge in CBR experiments and the many new numerical, theoretical and 
observational results on structure formation, the thermal history of the universe is now 
coming within better reach of our experimental probes. In view of this, it is very timely 
to theoretically investigate the nature of the reheating epoch in greater detail. It is also 
important to use observational data to place firm constraints on the thermal history of the 
universe, constraints that do not require any assumptions about uncertain entities like the 
primordial power spectrum. 

1.3 A Sneak Preview 

In this thesis, we will do the following: 

In Chapter 3, early photoionization of the intergalactic medium is discussed in a nearly 
model-independent way, in order to investigate whether early structures corresponding to 
rare Gaussian peaks in a CDM model can photoionize the intergalactic medium early enough 
to appreciably smooth out the microwave background fluctuations. We conclude that this 
is indeed possible for a broad range of CDM normalizations and is almost inevitable for 
unbiased CDM, provided that the bulk of these early structures are quite small, no more 
massive than about 10^ Mq. Typical parameter values predict that reionization occurs 
around z = 50, thereby suppressing fluctuations on degree scales while leaving the larger 
angular scales probed by COBE relatively unaffected. However, for non-standard CDM, 
incorporating mixed dark matter, vacuum density or a tilted primordial power spectrum, 
early reionization plays no significant role. 

In Chapter 4, the results of Chapter 3 arc generalized to CDM models where Q < 1 and 
the cosmological constant A > 0. Such models often require early reionization to suppress 
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degree-scale anisotropies in order to be consistent with experimental data. It is found that 
if the cosmological constant A = 0, the extent of this suppression is quite insensitive to 
J7o) as opposing effects partially cancel. Given a erg-normalization today, the loss of small- 
scale power associated with a lower is partially canceled by higher optical depth from 
longer lookback times and by structures forming at higher redshifts as the universe becomes 
curvature dominated at z ~ ^o^- The maximum angular scale on which fluctuations are 
suppressed decreases when Qq is lowered, but this effect is also rather weak and unlikely to 
be measurable in the near future. For flat models, on the other hand, where Aq = 1 — 
the negative effects of lowering Qq dominate, and early rcionization is not likely to play a 
significant role if f^o ^ 1- The same goes for CDM models where the shape parameter T is 
lowered by increasing the number of relativistic particle species. 

In Chapter 5, we find that it is still not certain that the universe underwent a neutral 
phase, despite the new COBE FIRAS limit y < 2.5 x 10~^ on Compton y-distortions of 
the cosmic microwave background. Although scenarios where the very early {z ~ 1,000) 
ionization is thermal (caused by IGM temperatures exceeding 10*^K) are clearly ruled out, 
there is a significant loophole for cosmologies with typical CDM parameters if the dominant 
ionization mechanism is photoionization. If the ionizing radiation has a typical quasar 
spectrum, then the y-constraint implies roughly h^^'^QhQQ^^''^^ < 0.06 for fully ionized models. 
This means that baryonic dark matter (BDM) models with ~ 0.15 and rcionization at 
z ~ 1,000 arc strongly constrained even in this very conservative case, and can survive the 
y test only if most of the baryons form BDM around the rcionization epoch. 

In Chapter 6, a model is presented in which supernova-driven winds from early galaxies 
reionize the IGM by 2; = 5. This scenario can explain the observed absence of a Gunn- 
Peterson trough in the spectra of high-redshift quasars providing that the bulk of these 
early galaxies are quite small, no more massive than about 10^ Mq. It also predicts that 
most of the IGM was enriched to at least 10% of current metal content by z = 5 and 
perhaps as early as z = 15. The existence of such early mini-galaxies violates no spectral 
constraints and is consistent with a pure CDM model with b < 2. Since the final radius 
of a typical ionized bubble is only around 100 kpc, the induced modification of the galaxy 
autocorrelation function is negligible, as is the induced angular smoothing of the CBR. Some 
of the gas swept up by shells may be observable as pressure-supported Lyman-alpha forest 
clouds. 

In Chapter 7, a formalism is presented that allows cosmological experiments to be tested 
for consistency, and allows a simple frequentist interpretation of the resulting significance 
levels. As an example of an application, this formalism is used to place constraints on 
bulk fiows of galaxies using the results of the microwave background anisotropy experiments 
COBE and SP91, and a few simplifying approximations about the experimental window 
functions. It is found that if taken at face value, with the quoted errors, the recent de- 
tection by Laucr and Postman (1994) of a bulk fiow of 689 km/s on scales of 150/i^^Mpc 
is inconsistent with SP91 at a 95% confidence level within the framework of a Cold Dark 
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Matter (CDM) model. The same consistency test is also used to place constraints that are 
completely model-independent, in the sense that they hold for any power spectrum whatso- 
ever — the only assumption being that the random fields are Gaussian. It is shown that the 
resulting infinite-dimensional optimization problem reduces to a set of coupled non-linear 
equations that can readily be solved numerically. Applying this technique to the above- 
mentioned example, we find that the Lauer and Postman result is inconsistent with SP91 
even if no assumptions whatsoever are made about the power spectrum. 

These chapters are all self-contained and fairly independent (with the exception of 3 and 
4), so they can be read in any order. 



1.3. A SNEAK PREVIEW 



7 




Year 



Figure 1.1: Limits on the Compton y-parameter. 

The observational upper hmit on the Compton y-parameter is plotted 
as a function of year. 
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Chapter 2 

A Cosmology Primer 



This chapter is by no means intended to be a thorough and systematic treatise on the 
basics of cosmology. It is merely an attempt to provide the bare-bones knowledge necessary 
for the reader with a general physics background to be able to follow the other chapters. 
Cosmologists are advised to skip directly to Chapter 3. 

2.1 Where to Read More 

For the reader interested in the basics of modern cosmology, there are many excellent books 
to chose from. Two early classics are 

• Physical Cosmology by Peebles (1971) and 

• Gravitation and Cosmology by Weinberg (1972). 
More recent texts include 

• The Early Universe by Kolb and Turner (1990), 

• Modem Cosmology by Dolgov, Sazhin and Zel'dovich (1990), 

• Introduction to Cosmology by Narlikar (1993) 

• Structure Formation in the Universe by Padmanabhan (1993) and 

• Principles of Physical Cosmology hy Peebles (1993). 

There is also a large number of good review articles, for instance in conference proceedings 
and printed lecture series. A very incomplete list includes the following: 

• Vatican Study Week: Large Scale Motions in the Universe, edited by Rubin & Coyne 
(1988), 

• Development of Large-Scale Structure in the Universe by Ostriker (1991), 

• Physics of the Early Universe, edited by Peacock, Heavens and Davies 1990 (see es- 
pecially the articles by Efstathiou and White), 
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• Observational and Physical Cosmology, edited by Sanchez, CoUados and Rebolo 1990 
(see especially the article by Bernard Jones), 

The detailed reference information for these publications is in the reference section at the 
end. Essentially all the results in the rest of this chapter can be found in books such as 
these, so very few explicit references will be given. 

2.2 So what is cosmology all about? 

Cosmologists study the history of the universe, from its birth some 10^° years ago up until 
today. Since the discovery of the cosmic microwave background radiation, the Big Bang 
model has been generally accepted by the scientific community. Although we still know 
very little about what happed before the universe was about a hundredth of a second old, 
the Big Bang model is widely believed to give a correct basic picture of what happened after 
that. Here are some of the cornerstones of this standard cosmological model: 

1. The universe is homogeneous and isotropic on large scales. 

2. The universe is expanding. 

3. The universe used to be dense and hot. 

4. The universe is full of microwaves. 

5. Structures formed through gravitational instability. 

6. The universe contains dark matter. 

We will now discuss these issues one by one. 

1. By large-scale homogeneity is meant that if we average over sufficiently large volumes, 
the density of the universe is roughly the same everywhere. This of course implies that the 
density is isotropic as well, i.e. that it appears the same in whatever direction we choose 
to look^. Not only the density, but all other fields as well (the velocity field, the metric, 
etc.) are assumed to be homogeneous and isotropic on large scales. There is now strong 
observational support for 1. Apart from this, we also tend to believe in homogeneity and 
isotropy because this greatly simplifies our calculations... 

2. All distant galaxies that we see in the sky are moving away from us, and the further 

away they are, the faster they tend to be receding. The homogeneity implies that on large 

scales, everything is indeed receding from everything else. In terms of classical physics, this 

is analogous to the way that all the raisins in an expanding raisin cake recede from one 

another, and just as in this analogy, the speed with which to objects are receding from 

each other is proportional to the distance between them. The proportionality constant H is 

called the Hubble constant. Its value is such that objects separated by more than about 10^*^ 

^Notc that homogeneity imphcs isotropy only for scalar quantities (such as the density), not vectorial or 
tensorial quantities such as the velocity or the metric. 
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light-years are receding from each other faster than the speed of hght. Doesn't this violate 
special relativity? Yes it does, but it's OK in general relativity... As will be discussed in the 
next section, rather than picturing the matter expanding through space at enormous speeds, 
a more natural way to think of this is to imagine the the matter as merely sitting still while 
space itself is expanding. The mathematical description of the expanding universe is called 
the Priedmann-Robertson- Walker (FRW) metric, and will be introduced in the next section. 

3. Extrapolating this expansion backwards in time, one concludes that the universe used 
to be much denser and hotter in the past. Indeed, the FRW metric tells us that some 10^° 
years ago^, the universe was denser and hotter than any object we have ever been able to 
study. Since we have no experimental knowledge about how matter behaves when it is so hot 
and dense, there was an early phase of the universe about which we cannot predict anything 
without making speculative extrapolations. After the universe was about a hundredth of 
a second old, however, the temperatures were no greater than those that we can study in 
stars and in laboratories, so the physical processes that operated after that arc fairly well 
understood and can be used to make detailed quantitative predictions. One such prediction 
is that most of the mass of the universe will be in the form of hydrogen, and about 25% will 
be helium. The fact that the observed helium fraction is indeed around 25% is one of the 
great triumphs of the Big Bang model. 

4. Another such prediction is that the universe is full of microwaves. These microwaves 
have roughly the same wavelength as those that heat our breakfast hot dogs, but their 
density is much lower and their spectral distribution is almost exactly that which would 
be emitted by a black object whose temperature was a few degrees Kelvin, i.e. just a few 
degrees above absolute zero. The discovery of this cosmic microwave background radiation 
in 1964 was another great triumph of the Big Bang model, arguably the greatest of them 
all. 

5. By studying these microwaves, we can infer that the density of the universe was 
almost perfectly uniform about half a million years after the Bang. Today, it is full of 
dense clumps, ranging in size from superclusters of galaxies to the planet we live on, not to 
mention the author. The main mechanism responsible for turning the primordial ordered 
state into such a mess is believed to be gravitational instability, whereby the gravity from 
regions slightly denser than average attracted surrounding matter, gradually creating ever 
larger lumps. This phenomenon is usually modeled as perturbations to the FRW metric, 
and will be discussed in Section 2.4. 

6. By studying the rotation of different parts of our galaxy, we can calculate how much 
mass there is at various distances from the galactic center. The surprising conclusion is 

^An amusing fact is that when Edwin Hubble first attempted to measure the constant that now bears his 
name, a false assumption led him to calculate a value more than five times larger than the current estimates. 
This implied that the universe was "only" about 2,000,000,000 years old. Since radioactive dating methods 
had showed certain meteorites to be about twice as old as that, the Big Bang theory came under fire from 
people wondering how rocks could be older than the universe itself. This age discrepancy was one of the 
reasons that Hoyle and others developed a competing (and now dead) theory, the so called steady state 
theory, in which there is no Big Bang — and indeed the reason that Hoyle condescendingly coined the term 
"Big Bang". 
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that there is much more mass present than can be accounted for by the matter we see. 
The mysterious invisible substance (or substances) is referred to as dark matter. There 
are also other indications of dark matter. Figuring out its nature and its density is one of 
the hottest problems in modern cosmology (Sadoulet &; Cronin 1991), as its gravitational 
presence affects the formation of structure mentioned in 5. 

The only parts of cosmology and astrophysics with which the reader should be familiar 
to be able to follow the following chapters are 

• The standard model of space and time in an expanding universe, including quantities 
such as the redshift z, the density parameter Q and the cosmological constant A 

• The basics of structure formation, including quantities like the random field (5(x) and 
the matter power spectrum P{k) 

• The microwave background, including quantities such as the the radiation power spec- 
trum and window functions 

• Rudimentary hydrogen chemistry 

• Some basic cosmology jargon, such as units like Mpc and fudge factors like h 



These will be the topics of the following five sections. 
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2.3 Bare Bones General Relativity and the FRW Metric 

In this section, wc will say a few words about general relativity and the large-scale structure 
of space-time, as well as list a number of handy formulas that will be used in subsequent 
chapters. 

2.3.1 General Relativity 

Journalist: "Professor Eddington, is it really true that only three people in 
the world understand Einstein's theory of general relativity?" 
Eddington: "Who is the third?" 

It would obviously be absurd to try to squeeze a real presentation of general relativity, 
one of the deepest physical theories ever to be invented by man, into a subsection of a 
Ph.D. thesis. In fact, I have been asked how so many cosmologists can use a theory as 
mathematically complex as "GR" in their work without even knowing basic differential 
geometry. I think that the answer is that some aspects of GR are quite intuitive, so that 
one can acquire a working feeling for what is going on even without being able to tell a 
Lie derivative from a Killing field. In addition to this, the FRW metric is so much simpler 
than the most generic solutions that what happens on distance scales much smaller than 
10^° hght -years can often be adequately described by classical physics plus some geometric 
considerations. In fact, in some ways GR resembles classical physics more than it resembles 
special relativity. 

2.3.2 The raw equations 

So how do you calculate things with GR? What is the basic mathematical structure of the 
theory? The basic objects are some fields that live on a four-dimensional manifold, and 
satisfy a bunch of messy nonlinear second order partial differential equations. The Einstein 
field equations are^ 

Rixu — -j^Rdiiv = ^irGTfj^i,, 

where 

R = g^^R/j.u, 

r> — r>a 

r>oc -pa -pQ; i -pT j^ct p7 T^oc 



and 5'^^ is the matrix inverse of g^n, i.e. 



number of different notational conventions can be found in the literature. Here we use the sign 

convention (+ — ) for the metric. In some texts, the signs of the last two terms in the definition of of 

R^^fj are reversed, which corresponds to a different convention for the index placement. 
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Here repeated indices are to be summed over from to 3, commas denote derivatives as 
in the standard tensor notation, and G is the gravitational constant. Throughout this 
section, we will use units where the speed of light c = 1. In the Einstein field equations, 
the dependent variables are the two tensors g^i^ and Tfj,^. They are both symmetric, and 
thus contain ten independent components each, g is called the metric tensor, and describes 
the structure of spacctime at each spacctimc point x^. T is called the stress- energy tensor, 
and describes the state of the matter (what is in space) at each point. The quantities F^^^, 
R'^iyfj and R^i, are named after Christoffel, Riemann and Ricci, respectively, and will not be 
further discussed here. 



2.3.3 The FRW metric 

Truth is much too complicated to allow anything but approximations. 

John von Neumann 



The differential quantity 

ds^ = g^i.dx^dx'^ 

is called the line element, and for the Priedmann-Robertson- Walker (FRW) metric takes the 
simple form 



ds^ = dt^ - a{tf 



dr^ 



+ r'^de^ + sin^ 9 dip^ 



1 — kr'^ 

where A; is a constant that is either zero or ±1, and a is a function to be determined later. 
For reasons that will soon become clear, we have renamed the independent variables 

x^ = {t,v), 

and changed to spherical coordinates 

r = {x\,X2,x^) = (sin cos (p, sin sin (p, cos 0)r. 

2.3.4 Interpreting the FRW coordinates 

The Newtonian equation of motion for a point particle moving along a trajectory x(t) in a 
gravitational field g is 

where i = 1,2, 3. The analog of this in General Relativity for a point particle moving along 
a trajectory x^{t) in spacetime, t now being any parameter, is the geodesic equation of 
motion 

xt" = T';;^x''x^. 

For the FRW metric, a cumbersome but straightforward calculation shows that 



= {t,ro), 
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is a solution to the geodesic equation for any constant tq. This means that if we interpret 
t as time^ and call the vector r the comoving spatial coordinates, then the FRW metric can 
be interpreted as follows: an object at a comoving position r = tq that is at rest (defined as 
r = 0) will remain at rest. The time t is simply the time that such an object would measure 
if it were a clock. By definition, the physical distance between two comoving points ri and 
r2 at some time t equals 

where the integral is to be taken along the geodesic, the curve from {t,ri) to {t,r2) that 
minimizes this quantity. A quick look at the expression for ds^ above tells us that this 
distance will be proportional to a{t). Thus objects that are at rest relative to the comoving 
coordinates nonetheless move relative to each other, the time-dependence of all distances 
scaling as the function a{t), which is named the scale factor. When our universe is modeled 
by the FRW metric, distant galaxies are assumed to stay at fixed comoving coordinates, 
and the universe expanding simply corresponds to a{t) increasing with time. The Hubble 
constant is evidently given by 

a 

(Note that calling it a "constant" is a misnomer, as it generally changes with time.) 

2.3.5 Curvature 

Let us for a moment keep t fixed and say a few words about the constant k. If k = 0, then 
the geodesic (shortest curve) between two points ri and r2 in will be the conventional 
straight line. This is not the case when k ^ 0. If we chose three points ri, r2 and 
and connect them in a triangle by the three geodesies, the sum of the three angles will be 
greater that 180° if A; = 1 (in which case space is said to be positively curved), equal to 180° 
if A; = (in which case space is said to be flat), and less than 180° if /c = — 1 (in which 
case space is said to be negatively curved or hyperbolic). If A; < 0, space is infinite at all 
times^, the coordinate r is allowed to range over the entire space R^, and the universe is 
said to be open. If k = 1, however, space has the topology of the 3-sphcrc, the universe 
is said to be closed, and its volume is finite for any t. In this case, the comoving spatial 
coordinates arc restricted to the range |r| < 1. (This coordinate patch covers exactly half 
of the three-sphere.) 

2.3.6 Redshift 

Let us now turn to the time-dependence. One can show that the momentum of all free 
particles decreases as the inverse of the scale factor a as the FRW universe expands. Thus 
massive particles gradually slow down and become comoving, "go with the flow" , ending up 
at rest with respect to the comoving coordinates. For photons, which cannot slow down, the 
momentum loss corresponds to a lowering of their frequency, i.e. they become redshifted. 

''This is a very fortunate feature of the FRW metric, since for generic metrics, it turns out to be impossible 

to define a global time that observers in different parts of the universe can agree on. 
^Except in unusual topologies that will not be discussed here 
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Thus the wavelength of a given photon is proportional to a, and we can whimsically think 
of the corresponding electromagnetic waves as being "stretched out" at the same rate that 
space is being stretched out. Another way to interpret this redshifting is as a Doppler shift, 
since due to the expansion of the universe, the distant objects emitting the photons that we 
see are receding from us. 

If we observe spectral lines in the spectrum of a distant object, they will all be shifted 
towards longer wavelengths by the same factor. Calling this factor (1 + z), the redshift of 
the object is defined to be the number z. For objects receding by a velocity w <C c, we have 
to a good approximation that 

V 

Z —. 

c 

If we assume that the recession velocity v is mainly due to the expansion of the universe, 
we can immediately estimate the distance d to the object as 

d fa -^z lO^^cm X z. 
H 

(If z is of order unity or larger, this formula is replaced by a more complicated one.) 

In cosmology, redshift is often used as a substitute for time. Indeed, the horizontal 
ax;is on most of the plots in the subsequent chapters is labeled z, not t. This is mainly 
because z corresponds much more directly to what astronomers actually measure. When we 
look out into space, we are looking into the past. If we observe a distant quasar, we know 
that we are observing events that occurred at some earlier time but we have no direct 
and accurate way of measuring or how long ago the quasar emitted the light that we see. 
However, we can measure its redshift simply by observing its spectrum and measuring where 
the spectral lines are. Thus "at redshift z" in cosmology jargon means "at the time when 
photons would have had to be emitted to have a redshift z by now." There is of course an 
explicit mathematical relationship between the two variables t and z, and this will be given 
later in this chapter. 2; = oo corresponds to f = (the Big Bang), and 2; = corresponds to 
today. 

2.3.7 The Friedmann equation 

So far, we have only discussed properties of the FRW metric. For it to correspond to 
something physical, it has to satisfy the Einstein field equations. If we take the stress- 
energy tensor T^i, that corresponds to comoving matter of uniform density p{t), substitute 
it into the Einstein field equations, and grind through a not particularly amusing calculation, 
it is found that the FRW metric is indeed a solution providing that the scale factor a satisfies 
the ordinary differential equation 



/a\2 SttG k 
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The solution to this equation, known as the Priedmann EquationP, clearly depends both on 
the value of k and on the equation of state, the way in which the density p depends on a. If 
the density is dominated by nonrelativistic matter (particles whose kinetic energy is much 
less than their rest mass), just "matter" for short, it will obviously drop as the inverse of 
the volume as the universe expands, giving 

p oc a~^. 

This has been the case for quite a while, for redshifts z <C 10^. If the density is dominated by 
photons (or other highly relativistic particles), which become redshifted as discussed above, 

-4 

p OC a 

since the energy of each photon dwindles as 1 /a as the universe expands. This was indeed 
the case very early on. A highly speculative third type of mass is vacuum energy, which if it 
exists has the property that its density stays constant, independent of the scale factor. The 
presence of such matter is cqTiivalent to the presence of an infamous "cosmological constant" 
A in the Einstein field equations. 



Let us define the critical density 



and the dimensionless variable 



n = 



Pcrit 

Then for the matter-dominated case (to which we limit ourselves in the remainder of this 
section), the Priedmann equation becomes 



- — Ho 
a 



1/2 



where if fc 7^ 0, we chose the scale factor today to be 

/ k _ 3000/i-iMpc 



which can be interpreted as the current radius of curvature of the universe. If = 1, then 
oo has no physical meaning and can be chosen arbitrarily — it will always cancel out when 
some physical quantity is computed. (Here and throughout cosmology, the subscript zero 
on a variable denotes its value today, at z = 0. For instance, Hq is the vahic of the Hubble 
constant today.) Aq is the above-mentioned cosmological constant, believed by many to be 
zero. The constant k was eliminated by using the fact that 

k = H^{Xo + no-l). 



^This is an example where general relativity is in a sense more "classical" than special relativity. The 
Friedmann equation is identical to that describing the time-evolution of the radius o of a self-gravitating 
cloud of uniform, pressureless dust in classical physics. Note that c, the speed of light, appears nowhere in 
the equation, so with this classical interpretation, nothing prevents the dust particles from moving faster 
than the speed of light. 
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This relation shows that the universe is fiat (A; = 0) if Aq + J^o = !> and closed or open 
depending on whether Aq + exceeds unity or not. In the standard cold dark matter 
(CDM) scenario, one assumes that = 1 and Aq = 0. For this simple case, the Priedmann 
equation has the solution 

a{t) oc t^/^, 

and n = 1 at all times. If the universe is open, it is clear that unless Aq < 0, we will always 
have d > 0, and the expansion will continue forever. If the universe is closed and Aq < 0, the 
expansion eventually ceases and the universe contracts back together again, ending with a 
hot Big Bang in reverse, a big crunch. For the Aq = case, the graph of the function a{t) is 
a cycloid, and the Priedmann equation above describes only the increasing half of the curve. 
Closed universes expand forever if (Glanfield 1966) 



Ao > 4f2o cos^ 



1 /r^-i .\ 47r 

- arccos [^q - IJ + — 



An recent review of properties of the Priedmann equation is given by Carroll et al. (1992). 

2.3.8 Some handy formulas for getting by in curved space 

This section contains a collection of useful formulas, many of which are used in the following 
chapters. It is intended as more of a reference section, and can be safely skipped by the 
general reader. 

Friedmann solutions, conformal time 

Since a^/a = 1 + z, the Priedmann equation (neglecting radiation) can also be written 

H = Ho^Xo + (1 + z)^l - Ao + ^oz) 

or, alternatively, 

dz 



— = Hoil + z)^\o + (1 + ^)2(1 - Ao + noz). (2.1) 

In the rest of this section, we will set Ao = unless the contrary is explicitly stated. This 
means that the Friedmann equation can be solved analytically, the SI > 1 solution being the 
ubiquitous cycloid. 

It is often convenient to introduce a new time coordinate r called conformal time, defined 

by 

/■* dt' 

In terms of the conformal time, the solutions a{t) to the Priedmann equation can be written 
in the parametric form 



Hot 



^^j^^(sinhT-r) forS^o<l, 
ife)' for Ho = 1, 

I 2(f)o-i)3/2 ' 



(r-sinr) for J^o >1, 
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a 
ao 



1 + z 



< forl7o = l, 

When ^Iq ^ 1, T is often called the development angle. When J7o = 1) '^o = ^/{Hoao), where 
the constant ao can be chosen arbitrarily. Again, this is because the radius of curvature, 
which fixed oq when Qq ^ 1, is infinite when the universe is flat. In all cases, r = 
corresponds to the Big Bang. Inverting these last expressions, we obtain 

2sinh^^i 



n^IT^ for ^0 < 1, 

for = 1) 



.2arcsiny^5j^ for J^o > 1- 



Thus the conformal time today is 



To 



2sinh-Yf^o ^ - 1 for ^^o < 1, 
2arcsin y''l — Qq^ for Oo > 1- 

Although one cannot give an explicit expression for a in terms of t in general, 

\ 2/3 



ao 



for t < Hq\ 



Causality, horizons, angles 

A useful property of r is that during a conformal time interval dr, a light ray always travels 
the same comoving distance dr. If $7 = 1 so that space is flat, the situation becomes 
especially simple when viewed in (r, x)-space: all light rays (null geodesies) simply move 
in straight lines, and always make 45 degree angles with the r axis. Thus after changing 
variables to comoving spatial coordinates and conformal time, we can use all our intuition 
about light propagation and causality in good old non-expanding Euclidean space. This is 
illustrated in Figure 2.1. 

Parts of the universe receding from us faster than the speed of light are not in causal 
contact with us. No signals have ever been able to reach a comoving object in the universe 
from objects that are outside of its backward light cone. Such objects are said to be outside 
of its horizon. Thus defining the horizon radius as ar and combining some of the previous 
expressions, we obtain 



rh = ar 



(i+z)7i=r^ 

(i+2)v^niFT 



cosh 



1 + 



2(i-no) 



cos 



-1 



2(no-i) 

^0(1+2;) 



if JIq < 1, 
if flo = 1, 
if Qo > 1. 
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If an object at redshift z has a proper diameter D (perpendicular to us), then one can show 
that it win subtend an angle on the sky given by 



2 tan - = < 
2 



' ng(l+z)sinh[l(l+z)v^]^ifoD] -r Q < 1 

v/r^[no^-(2-Qo)(Vi+n^-i)J ^ ° ' 

rif ifOo = l, 

Vn^[no2-(2-Qo)(Vi+n^i^ ° 



Thus if we look at an object at a redshift z, its horizon radius subtends an angle 6 in the 
sky that is given by 

2tan^ = ^TV^- 

2 noz-{2-no){VT+u^-i)' 

for all values of CIq. For = 1, this takes the particularly simple form 

e 1 

2 tan - 



v^TT^-l' 




whereas for large z, it reduces to 



This last expression is a rather poor approximation for for typical values of JIq, unless 
z » 100. 

Age of the universe 

The current age of the universe, to, is given by 



to 



dt 



{aH)-\z) 



da 
dz 



dz 



dz 



{l + z)H{z)' 



Performing this integral for Aq = gives the age of the universe as 

■ cosh~^ (^-^ - ij for Qq < 1, 
for = 1, 



Hoto = < 



1 ^0 

l-fio 2(l-f2o)3/2 
2 
3 



Op 



cos 



I. 2(Qo-l)3/2 

For the fiat case where k = and J7o = 1 — Aq, the result is 



1/2 



Hoto 



,1/2 



In 



I + Aq 
(l-Ao)V2 
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2.4 Perturbing the Universe 

In this section, wc will give the rudiments of the theory of structure formation. First we 
discuss how small ripples in the FRW metric grow larger as the universe expands. We 
then briefly mention the theory of Gaussian random fields and the so called Press-Schechter 
model for galaxy formation. 

2.4.1 Linear perturbation theory 

Let us assume that demand is an arbitrary function of price; D = aP + h. 

Anonymous economics professor 

In this section, everything will indeed be linear. As mentioned in the previous section, 
it is generally futile to look for analytic solutions to the Einstein field equations of GR. 
Fortunately, the actual metric is fairly uniform on large scales, indicating that the relative 
deviations from the FRW metric are rather small. For this reason, cosmologists attack the 
problem of structure formation with linear perturbations. One perturbs the FRW metric, 
throws away all terms of second order and higher, and after Fourier transforming ends up 
with a simple equation for the time-evolution of the various perturbation modes. 

As mentioned in a footnote a few pages back, the Friedmann equation can be derived 
using classical physics alone. This is no coincidence. A result known as Birkhoff 's theorem 
states that if we take an FRW metric and remove all matter from a spherical volume, 
then space will be flat inside of this sphere, i.e. all gravitational effects from the matter 
outside cancel out. If we put the uniformly expanding matter back into the sphere, then 
its gravitational behavior will be correctly described by classical physics. As long as the 
gravitational field is weak (which it is in the present problem, as opposed to say near a black 
hole) and the sphere is much smaller than the horizon scale lO^'^ light years, it is easy to show 
that classical physics will apply even if the matter is not uniform. With this motivation, let 
us write down the classical equations governing a gas of density p, velocity v and pressure 
p and the corresponding gravitational potential (f). These are the continuity equation, the 
Euler equation of motion and the Poisson equation of classical gravity, respectively: 

p + V ■ (/9V) = 

V + (v V)v = -V((/) + 2) (2.2) 

A simple solution to these equations is 

>o(i,r) = ^, 

< vo(t,r) = |r, 

<^,(t,r) = 2ZE|Po^2^ 

where a is a function of time that satisfies the Friedmann equation, and since the density 
is independent of r, this classical solution clearly corresponds to the unperturbed FRW 
solution in GR. From here on, let us write these fields as functions of the comoving position 
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X rather than the physical position r = a(t)x. Let us now expand the fields as 

' p{t,^)=Po{t)[l + 5{t,^)], 
< v(t,x) = vo(t,x) + vi(t,x), 
. x) = (/)o(x) + 4>i{t,x.), 

and assume that \6\, |vi|, 0i ^ 1. Substituting this into equation (2.2), dropping all non- 
linear terms, Fourier transforming everything with respect to x, and doing some algebra, 
one obtains the second order ordinary differential equation 

t(k) + 2^t(k) + 1^^^ - AirGpo^ S{k) = 0. (2.3) 

Here the equation of state for the gas enters only through the sound speed, 



\dpJadiab ~ V 3mp' 

where the last expression applies if the gas consists mainly of atomic hydrogen (rrip denotes 
the proton mass). If wc substitute an expanding universe solution for a, equation (2.3) tells 
us a number of interesting things. 

The coefficient of 6 has two contributions with opposite signs: the destabilizing effect of 
gravity competes with the stabilizing effect of gas pressure. On small scales (corresponding 
to large k), pressure dominates and perturbations do not grow, merely oscillate. On large 
scales, on the other hand, self-gravity dominates, and fiuctuations grow monotonically until 
S becomes of order unity and the entire linear approximation breaks down. The border 
between these two regimes is when the perturbation has a physical wavelength 

_ 27ra 

called the Jeans wavelength 



X.T = V, 



TT 



for which this coefficient vanishes. 

The coefficient of S, twice the Hubble constant H, acts like a kind of friction: in the 
absence of this term, we could get exponentially growing modes, but as will be seen, the 
growth is never that rapid in the presence of this term, i.e. in an expanding universe. 

If the density field is smoothed on a scale A ^ Aj, then we can neglect the ^-dependence 
and trivially Fourier transform back to real space. Thus all contributing modes grow at the 
same rate, and a perturbation retains its (comoving) shape as the universe expands, simply 
increasing in amplitude. For the simple case where O = 1, Aq = and a oc we obtain 

which has the most general solution 

<5(t,x) = A(x)<5(io,x) B{x)S{to,^) (^^) 
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for some time-independent functions A and B. Ignoring the physically uninteresting de- 
caying mode, and writing 5 as a function of redshift rather than time, we thus have the 
extremely simple result 

where 5o is the linearly extrapolated field of density perturbations today. 



2.4.2 Random fields 

The mathematical theory of random fields (sometimes known as three-dimensional stochastic 
processes) is a very useful tool when analyzing cosmological structure formation. A random 
field is simply an infinite-dimensional random variable, such that each realization of it is a 
field on some space. In cosmology applications, these fields tend to be S, v or (f), and the 
space is physical space at some fixed time t. As every quantum field theorist knows, it is a 
nightmare to try to define a nice measure on an infinite-dimensional space, so random fields 
are defined by specifying the joint probability distribution of their values at any n points, 
n = l,2,3,..., thus circumventing the need to define a probability distribution on the space 
of all fields. Hence, to define a random field 5, one must specify the 1-point distribution 
(the 1-dimensional probability distribution of 5{xi) for all xi), the 2-point-distribution (the 
2-dimensional probability distribution of the vector [(5(xi), (5(x2)] for all xi and X2), etc. 
In cosmology, the random fields are always assumed to be translationally and rotationally 
invariant, i.e. homogeneous and isotropic. Hence the 1-point distribution is independent of 
xi, and the 2-point distribution will depend only on the scalar quantity x = |x2 — xi|. 



Ergodicity 

A random field S is said to be ergodic if one can use ensemble averaging and spatial averaging 
interchangeably. The ensemble average of a random field (5 at a point, denoted ((5(x)), is 
simply the expectation value of the random variable S{x). Thus for an ergodic field, 

(5(xi)) = lim (iTTRA ^ [ 54x)d^x 

il-»oo Vo / J\x\<R 

holds for all points xi and for all realizations 5*(x) (except for a set of probability measure 
zero). Ensemble averages are completely inaccessible to us, since we have only one universe 
to look at, namely the particular realization that we happen to live in. So as cosmologists, 
we are quite happy if we have ergodicity, since this means that we can measure these elusive 
ensemble averages by simply averaging over large volumes instead. 



Gaussianity 

A random field is said to be Gaussian if all the above-mentioned probability distributions arc 
Gaussian. This is a very popular assumption in cosmology, partly because, as we will sec, it 
greatly simplifies matters. A first nice feature of this assumption is that all homogeneous and 
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isotropic Gaussian random fields are ergodic^. Taking the spatial average of the definition 
of S, for instance, ergodicity implies that 

(<5(x)) = 0. 

Let us define the r.m.s. fluctuations 
and the correlation function as 

('5(x2)<5(xi)) 



t2 



(Note that because of the homogeneity and isotropy, a is independent of x and the correlation 
function depends only on x = |x2 — xi|.) Since the n-point distribution is Gaussian, it is 
defined by its mean vector ((5(x„)) (which is identically zero) and its covariance matrix 

Cmn = (<^(Xm)(^(x„)) = 0-^^(|x„ - X„|). 

Thus the Gaussian random field S has the extremely useful property that it is is entirely 
specified by a and its correlation function. 

The power spectrum 

If we Fourier expand 6 as 

'^^"^ = (2^ / ^^'•^^(k)'^'^' 
we see that if its Fourier transform 



m = I < 



is a Gaussian random variable for any k, then S will automatically be a Gaussian random 
field (since a sum of Gaussians is always Gaussian). Cosmologists like to postulate that the 
complex numbers (5(k) have random phases, which implies that they are all uncorrelated. 
One postulates that 

(?(k)*^(k')) = (27r)=^,5(k-k')P(k), 

where the function P(k) is called the power spectrum and 6 is the Dirac delta function 
(which will hopefully not be confused with the random field 6). This implies that even 
if (5(k) does not have a Gaussian distribution, the random field S, being an infinite sum 
of independent random variables, will still be Gaussian by the Central Limit Theorem for 
many well-behaved power spectra. It is straightforward to show that the power spectrum is 
simply cr^ times the three-dimensional Fourier transform of the correlation function, i. e. 



P{k) = AiTu'' j {^^^i{rydr. 



Note that P depends on k only through its magnitude k = |k|, because of the isotropy 
assumption. 

^Notc that this is true only for random fields that live on infinite spaces such as R", and does not hold 
for fields on compact manifolds such as the sphere . Thus the field of microwave background amisotropies 
(defined further on) is not ergodic, so that even if we could reduce our experimental errors to zero, we could 
still never measure any ensemble-averages with perfect accuracy. This phenomenon is known as "cosmic 
variance". It stems from the fact that in the spatial average above, one cannot average over an infinite 
volume {i.e. let R — > oo), since the volume of the compact manifold is finite. 
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Window functions 

If we smooth S by defining the weighted average 

5^(x) = I <5(r'Hr' - r)dV, 
w being some weight function, then 

al ^ (|5^(x)|2) = j^J \w{k)fmd'k. 

This is quite a typical pattern: the variance of some physical quantity is given by an in- 
tegral of the power spectrum against some function, referred to as a window function. In 
this particular example, where for an appropriately chosen w the physical quantity can be 
interpreted as the mass within a volume of some given shape, the window function is \w\'^. 
As we will see in Chapter 7, basically all cosmological experiments that measure a single 
number can be described by some window function, regardless of whether they involve mass 
distribution, bulk flows of galaxies or microwave background fluctuations. 

2.4.3 The Press-Schechter recipe 

In this section, we will describe a rather simple-minded model of structure formation that 
agrees remarkably well with observational data. 

The top hat solution 

Early on, while z S> ^' space is approximately flat and the Priedmann equation has the 
approximate solution 

a{t) cx t2/3 

regardless of the values of JIq and Aq. If an f2 = 1 universe has a completely uniform density 
p except for a "top hat" overdensity, a spherical region where the density at a given time is 
some constant p' > p, then this top hat region will gradually begin to expand slower than 
the rest of the universe, stop expanding and recoUapse to a point. By Birkhoff's theorem, 
the radius of this region will evolve according to the Priedmann equation, but with J7 > 1. 
It is easy to show that the overdensity 



evolves as 



9(a — sina)^ -, 3 o „, 
"+''-2 il-oos„)3 =^+20'' 
where the parameter a, the "development angle" , is related to the redshift through 

Here z^r is the redshift at which the top hat would collapse to a point. In reality, an over- 
dense region would of course not collapse to a point (and form a black hole). Since it would 
not be perfectly spherically symmetric, dark matter particles would mostly miss each other 



26 



CHAPTER 2. A COSMOLOGY PRIMER 



as they whizzed past the central region and out again on the other side, eventually settling 
down in some (quasi-) equilibrium configuration known as a virial state. For baryons, gas- 
dynamical processes become important, and pressure eventually slows the collapse. Strictly 
speaking, virial states are not stable over extremely long periods of time, and their density is 
certainly not uniform. For a virializcd lump, often referred to as a "halo", a typical density 
profile peaks around some constant value in its core and falls off like over some range 
of radii. Nonetheless, people often say that they have a "typical" density 

Pvir ~ 187r^po(l -I- Zyir)^, 

which is a useful rule of thumb. Thus in the top-hat collapse model, density in the perturbed 
region is assumed to evolve as in Figure 2.2. The density starts out decreasing almost as 
fast as the background density p, with 

early on, just as in linear theory, but gradually stops decreasing and increases radically as 
z approaches z^ir- It never increases past the virial value 187r^po(l + z-^ir)^, however, 
but stays at that density for all z < z^ir- Unfortunately, a cannot be eliminated from the 
equations that relate 5 and z, i.e. equations (2.4) and (2.5), by using elementary functions. 
For this reason, I made the following fit to the function 5{z), which is accurate to about 5% 
until z is within 10% of Zyir, at which the density is assumed to start approaching the virial 
value anyway: 

p{z) K, po(l + 2;)^exp 

where 

A{z) = — 

The Press-Schechter formula 

In the above top-hat solution, we have approximately 

5= |-(127r)2/3yl« 1.686yl 

early on, while A <^ 1. If we naively extrapolate this far beyond its region of validity, we 
thus obtain 

5 = 5c = ^(127r)2/3 ^ 1.69 

when A = 1, i.e. when z = Zmr- Press and Schechter proposed the following simplistic pre- 
scription for modeling structure formation: define a smoothed density field 5m by averaging 
5 over spherical regions of such a radius that they contain a mass M on the average. Define 
a{M) as the r.m.s. mass fluctuation in these spheres, i.e. 

a{M) ^ {5l)y\ 

Now make the interpretation that a point x is part of a virialized halo of mass exceeding M 
if 

(^m(x) > 5c ~ 1-69. 



1.9 A 



1 - 0.75^42 
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In linear theory, for = 1 and Aq = 0, we clearly have 

cj(M, 0) 



a{M, z) 



1 + z 



Since (5 is a Gaussian random field, so is 6m, which means that Sm{^) is a Gaussian random 
variable with zero mean and a standard deviation a{M,z). Thus with the Press- Schechter 
interpretation, the fraction of all mass that is in virialized lumps of mass exceeding M at a 
redshift z is 

' {l + z)Sc 



/g(M,^) = P[(5M(x)>4] = ^erfc 

where the complementary error function is defined by 

2 



V2a{M, 0) 



2 r°° 2 
erfcfx] = —= / dx. 



Unfortunately, the so defined fg{M, z) approaches ^ as t — oo, i.e. half of all mass remains 
unaccounted for. To remedy this. Press and Schechter introduced a "fudge factor" of 2 and 
predicted that 

/,(M,.)=erfc[4^±^' 

And lo and behold: this is quite a good fit both to the observed mass distribution of galaxies 
today and to the fg{M,z) that is observed in numerical n-body simulations. Although a 
score of papers have been written on the topic since this formula was proposed some twenty 
years ago, often attempting to justify the factor of two, not much has happened to the basic 
formula except that other values have been proposed for the constant b^, ranging from 1.33 
to 2.00. 



2.4.4 CDM and BDM 

After all this talk about the power spectrum P(A;), what does theory predict it to be? There 
is now a plethora of rival theories on the market predicting all sorts of power spectra that 
are tailor-made to fit various observational facts. In this extremely cursory discussion, we 
will give only the highlights of two models. The most popular model over the last decade 
has probably been the cold dark matter (CDM) model with adiabatic initial fluctuations 
and an n = 1 Harrison Zel'dovich initial power spectrum. We will refer to this mouthful as 
"standard CDM" for short. This model is simple and elegant, as it contains almost no free 
parameters, and roughly speaking fits all the observed power spectrum data to within about 
a factor of two when k ranges over many orders of magnitude. A recent rising star (albeit 
always surrounded by rumors of its imminent death) is the baryonic dark matter (BDM) 
model with isocurvature initial flTictuations and a power-law initial power spectrum. We 
will refer to this as just BDM for short. 

All theories of structure formation (theories that predict P{k)) split into two parts: one 
part that gives some "initial" or "primordial" power spectrum Pi{k) when the universe was 
very young, and a second part that describes how this power spectrum is altered by various 
physical processes to produce the P{k) we observe today. Neglecting speculative theoretical 
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constructs such as cosmic strings and other so cahed topological defects, the latter processes 
usually do not mix different Fourier modes as long as S <^ 1. Thus the second half of a 
theory can be specified by a single function, the transfer function T{k), defined so that 

P{k) = T{kfPi{k). 

The initial power spectrum 

For lack of any better assumptions, cosmologists usually assume that the primordial power 
spectrum was a simple power law 

P{k) oc k"" 

for some constant n termed the spectral index, at least over the range of A;-values that are 
accessible to our observations. In CDM, the primordial fluctuations are usually assumed to 
arise from quantum fluctuations during a speculative inflationary epoch when the universe 
was less than 10~^^ seconds old. Thirteen years after the theory of inflation was invented by 
Alan Guth, there are now scores of different versions of it. Most of them predict a spectral 
index between 0.7 and 1. These initial fluctuations are adiabatic, which means that the 
density varies slightly from point to point but that the ratio of photons to baryons and dark 
matter remains constant initially. For BDM, the primordial spectrum is an Achilles' heel, 
as the theory contains no physical mechanism for generating it. One simply postulates that 
it is a power-law. In addition, one postulates that the initial fluctuations are of isocurvature 
type, which means that the fluctuations in the photon density cancel the baryon fluctuations, 
so that the total density remains constant. One may question the value or calculating the 
BDM transfer function, since at the end of the day it is multiplied by a function that was 
put in arbitrarily by hand, but peoples' opinions differ on what is elegant. An aesthetic 
advantage of BDM is that it requires no new particles besides those we know exist, whereas 
CDM postulates that more than 90% of the matter is of a type that we have never detected 
in our laboratories. 

The CDM transfer function 

The standard way to compute transfer functions is to write down the Boltzman equation 
for the time-evolution of the phase-space distributions of all particle species in the model 
(baryons, photons, neutrinos, collisionless dark matter, etc.), linearize it by throwing away 
all quadratic and higher order terms, plug it into a fast computer and go for lunch. To 
model evolution of the power spectrum at late times, when the approximation 6 <^ 1 breaks 
down, so called A^-body simulations are employed: the time-evolution of some 10^ to 10^ 
point-particles in a box is studied by numerically integrating the Newtonian equations of 
motion, assuming periodic boundary conditions. The power spectra resulting from a few 
assorted models are plotted in Figure 2.3. 

Fortunately, the main features of the CDM spectrum can be understood without nu- 
merical calculations, with a bit of hand-waving, as wc will now sec. If the transfer function 
is to be anything else than a simple power law in k, some process must imprint a physical 
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length scale on it. Collisionless dark matter particles have no particular length scale at all 
associated with them that is relevant in this context. However, at redshifts z > 10^, the 
main contribution to the density and the self-gravity of the universe was CBR photons, not 
dark matter, so before this, the crucial issue was what the photons decided to do. Were it 
not for this fact, a power-law initial power spectrum would necessarily lead to a power-law 
spectrum today. In fact, the transfer function would be simply T{k) = 1, giving 

P{k) oc A;". 

As the universe expands, the horizon size keeps increasing, also in comoving coordinates. 
This means that the wavelength of any sinusoidal perturbation (a distance which remains 
constant in comoving coordinates) will sooner or later become smaller than the horizon 
scale. When this happens, the perturbation is said to "enter the horizon" . To avoid various 
painful GR gauge ambiguities associated with super-horizon sized modes, it is convenient 
to specify the amplitude of perturbations with a given wavenumber k when they enter the 
horizon. Since microphysical processes such as pressure support turn out to be irrelevant 
for super-horizon sized modes, the power spectrum at horizon entry must also have been 
some power law. 

Before recombination, the photons and baryons were tightly coupled and acted as a single 
fluid where the sound speed was of the order of the speed of light. This means that the 
above-mentioned Jeans length was enormous, indeed several times greater than the horizon 
size, so perturbations that had entered the horizon did not grow much at all early on, indeed 
logarithmically at best. As the redshifting of the photons caused the density of the universe 
to become matter-dominated around z ~ 10^, the perturbations in the dark matter began 
to grow. And at z ~ 10'^, recombination caused the Jeans length to drop drastically, down 
to a comoving scale much lower than the sizes of prcscnt-day galaxies, and perturbations in 
the baryons began to grow as well. Thus dark matter fluctuations that entered the horizon 
long after matter-radiation equality, fluctuations whose comoving wavelength exceeds 

Xeq ~ 13(J^o/!'^)"^MpC, 

were never affected by any of this microphysics that has a length scale associated with it, 
which means that T{k) 1 for A; ^ X~q. Fluctuations with k S> X~q, on the other hand, 
entered the horizon while dark matter was still only a minuscule contributor to the overall 
density, and did not start growing properly until after matter-radiation equality. Thus we 
know that for k S> X^q, T{k)'^ must be some different power law k'^, with m < 0. We will 
now sec that m = —4. Take Q = 1 and consider a perturbation of comoving wavelength A. 
Then its physical scale is 



whereas the horizon scale is 

dij octoc (l + z)-^/'^. 
Equating these two scales gives the redshift for horizon entry 

(1 + Zent) OC A"^ (X fc^. 
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Since this redshift, S has grown by a factor of (1 + Zent), so P{k) has grown by a factor of 
(1 + Zent)^ fc^- The modes with k S> A"^-*^, on the other hand, have all grown by the same 
factor since they entered the horizon, so in this part of the power spectrum, the slope is 
down by the sought for factor of k^. In conclusion, the net result of all this is that 



where X^q ~ 50 Mpc for O = 1 and h = 0.5. 

The transfer functions for baryons and photons 

All this concerned the power spectrum of the density of dark matter. The baryons and the 
photons are affected by two additional effects, known as Landau damping and Silk damping, 
which for the case of adiabatic fluctuations effectively wipes out their fluctuations on galactic 
scales. Since we know that galaxies do indeed exist, a successful theory of structure formation 
must find a way around this disaster. With CDM, the trick is that the cold dark matter is 
unaffected by this damping, and that gravitational instability eventually makes the baryon 
fluctuations "catch up" with the fluctuations in the dark matter. With BDM, the trick is to 
utilize isocurvature fluctuations instead, which do not suffer Silk damping. We will return 
to the transfer function for photons in the CBR section. 

The complaint department 

As to which model of structure formation is correct, the jury is still out. We will conclude 
this section by criticizing all the models, as democratically as possible. According to many 
workers, standard CDM is off by about a factor of two. This may not be much to fuss about 
in a field where one is often satisfied with getting the right answer to within an order of 
magnitude, but as data gets better, concern grows. This factor of two can be dealt with in 
at least six ways: 

• One can argue that what we see is a biased version of what is actually there, and 
introduce different fudge factors called bias factors for various types of measurements. 

• The spectral index n can be lowered to around 0.7 by tweaking certain parameters 
associated with inflation, which seems to underproduce high redshift quasars. So do 
several of the other models listed below. 

• Repeating what Einstein described as "the biggest blunder of my life" (Gamow 1970, 
p. 44), a non-zero cosmological constant A can be introduced. 

• Yet another free parameter can be introduced by adding just the right amount of 
massive neutrinos (about 30%) to the cosmological mixture. 

• Xeq can be increased by postulating that the r neutrino decays in an appropriate 
manner and that there arc new and unseen particles, for instance additional massless 
Nambu-Goldstone bosons. 




for k < Xgg, 
for k » Xgg, 



2.4. PERTURBING THE UNIVERSE 



31 



• Xeq can be increased by lowering the Hubble constant to about h = 0.25, causing 
raised eyebrows among those astronomers who claim that h = 0.8. 

In addition, all of these inflation-based theories are challenged by the fact that some CBR 
experiments have been argued to indicate n ~ 1.5 on large scales, whereas essentially all 
inflationary scenarios predict n < 1. 

Meanwhile, BDM is fighting a difficult battle, ensnarled by a number of tight constraints. 
The Compton y-parameter prevents reionization from being too complete, too hot and too 
early (see Chapter 5), degree-scale CBR-experiments prevent reionization from being too 
late, and the observed power-spectrum constrains the thermal history as well. As the author 
is writing this, one of his colleagues is running a Boltzmann code for BDM models which 
have nine free parameters, so it appears that in dodging the many constraints, the BDM 
model has lost much of its initial elegant simplicity, and with that perhaps most of its 
appeal. 
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2.5 The Microwave Background 

In 1964, Arno Pcnzias and Robert Wilson at Bell Laboratories discovered microwave radia- 
tion of an unknown source that they at one point thought might have been bird droppings in 
their antenna horn. Eventually, their discovery was given a more grandiose interpretation, 
which resulted in them both winning a free trip to Stockholm. This radiation is the topic 
of the present section. In the first few subsections, where its origin is described, we will also 
make a digression about the thermal history of the universe. 

2.5.1 The thermal history of the universe 

When electrons are accelerated, they emit photons. Thus if one heats a hydrogen plasma 
sufficiently, the fast-moving electrons will frequently scatter off of nuclei and other elec- 
trons, resulting in a significant number of photons being created and radiated away. For 
cosmological applications, the two most important of these scattering processes are thermal 
bremsstrahlung 

+ e~ 1-^ p'^ + e~ + J 

and double Compton scattering 

e~ + 7 I— > e~ + 7 + 7. 

These processes can also go in the reverse direction, so if such a hot plasma filled the 
entire universe, some of these photons would eventually get into collisions where they would 
become absorbed. One can show that if the plasma is held at a constant, uniform density 
long enough, then a thermal equilibrium situation will arise where the electrons and protons 
have Maxwell-Boltzman velocity distributions corresponding to some temperature Tg and 
the photons have a Planck blackbody distribution 

where Tj = Tg. Here the photon frequency u = |k|c/27r, and the wave vector k is not 
to be confused with Boltzmann's constant k. Note that since the number density n-y is a 
distribution over the six-dimensional photon phase space (x, k), it is dimensionless, as the 
units of m""^ arising from x are canceled by the m^ arising from k. 

In the beginning, there was thermal equilibrium 

This was not quite the situation in the early universe, since the expansion kept lowering 
the density (and hence the reaction rates). However, it is easy to show that there was 
nonetheless sufficient time for thermal equilibrium to be attained early on. Thus at high 
redshifts (say z S> 10^), the entire state of the hydrogen plasma and the photons could be 
described by a single number, the temperature T = Tg = T^. Indeed, if we go back to 
high enough redshifts, thermal equilibrium was attained not only for electron and proton 
velocities and for electromagnetic radiation, but indeed for all particles that we know of, as a 
plethora of reactions kept rapidly creating and destroying all sorts of particles left and right. 
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For instance, when the universe was much less than a second old, electrons and positrons 
where being created and destroyed so rapidly that their densities where accurately given by 
the standard Boltzman formula. And even earlier on, the same applied to heavier particles 
like protons and W bosons. 

Freezeout 

As the universe expanded, the reaction rates that affected the abundance and energy dis- 
tribution of any given particle species kept dropping until that species "decoupled" from 
the merry thermal mixture, or "froze out" as cosmologists say*. One by one, species after 
species froze out. If the rest mass of a particle was much greater than the thermal energy 
kT when its species froze out, then tough luck: the freezeout abundance would have been 
exponentially suppressed by a huge Boltzman factor, and almost no such particles would 
still be around today to bear witness. This was the grim fate of the W'^, W~ and Z bosons, 
which constitute only a minuscule fraction of the density of the universe today. This is also 
the assumed fate of many highly speculative particles^ whose existence today would be in 
conflict with observations. In the other extreme, where the particle is so weakly interacting 
that it freezes out while kT is still larger than its rest mass, there is still an ample abun- 
dance of such particles today. Such weakly interacting massive particles, known as WIMPS, 
are interesting candidates for the dark matter in the universe. Another classic example 
of freezeout survivors are the cosmic background neutrinos, which arc still with us today. 
A most interesting intermediate case is the neutron. Being about 0.2% heavier than the 
proton, the neutron to proton ratio in equilibrium is 

rin „-ioiok/t 

Up 

Thus the neutron and proton abundances would be almost identical while the temperature 
of the universe exceeded their energy difference, but eventually, as T ^ 0, the ratio would 
approach zero exponentially. Thus we might expect to find either virtually no neutrons today 
or a 50% abundance. In the latter case, it turns out that virtually all the hydrogen would 
have formed Helium during nucleosynthesis, and we would find things such as water quite 
difficult to come by. If we were around... Interestingly, as mentioned in the nucleosynthesis 
section above, nature chose to be right in the middle: neutron freezeout occurred just around 
this critical temperature, and the primordial abundance is about 12%. 

After this digression, we will limit ourselves entirely to electrons and photons. 



a more quantitative definition is to say that a given reaction is frozen out when the rate F at which a 
particle is subject to such reactions is so low, that less than one such event is expected to occur while the 
universe doubles its age. This happens roughly when T = H, the Hubble constant. 

^Particle physicists desperate for funding have postulated the existence of many dozens of particles never 
seen by man nor beast... 
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Recombination and the CBR freezeout 

The freezeout picture applies to hydrogen recombination as well. In thermal equilibrium, 
the fraction of hydrogen that is ionized is 

-1 



X 



1 + 2.8x10-6t-'/V-5«/^« 



where = T/IO^K, which means that if one slowly lowers the temperature of a hot hydrogen 
plasma, it will recombine when T Ri 15, OOOK. However, the recombination reaction 

p'^ + e~ 1-^ H 

has already begun freezing out by then. Quantum effects such as stimulated emission com- 
plicate the process further, and the end result is that recombination does not occur until 
much later, at a temperature of a few thousand K, ai z ^ 10^. At this point, something 
else of great importance happens as well: the photons freeze out. Until then, they were 
constantly Thomson scattering off of free electrons, but from this point on, a generic photon 
merely moves along in a straight line (or, more precisely, a geodesic) until today. Another 
way of phrasing this is that the universe became transparent at z 10^, so that in a 
sense we can see that far back. These frozen out photons gradually become redshifted as 
was described in the GR section, and it is easy to show that this redshifting preserves the 
blackbody nature of their spectrum, simply lowering the blackbody temperature as 

Xy = (1 + Z)T^Q 

as the the universe expands. It is these frozen out photons that are known as the cosmic 
microwave background radiation, the CBR. The best measurement to date to puts their 
current temperature at 

T^o ~ 2.726K ± 0.005K, 

and was made by the Far Infra-Red Absolute Spectrometer (FIRAS) detector on board the 
celebrated COBE satellite (Mather et al. 1994). 

2.5.2 CBR fluctuations 

If one corrects for our motion relative to the comoving CBR rest frame^'^, the observed CBR 

is almost perfectly thermal and isotropic. In other words, its spectrum is almost exactly a 

Planck spectrum, and the temperature Tjq is almost exactly the same in all directions in 

the sky. If it were not for the caveat "almost" , the CBR would not be so interesting, as it 

would contain no more information than the single number T^q. (In addition, a perfectly 

isotropic CBR would have implied that all our theories about structure formation would be 

dead wrong, as we shall see below.) It is the small deviations from thermality and isotropy 

that have the potential to give us information about the early universe. 

^°Note that whereas there are no preferred observers in classical physics or in generic GR metrics, which 
renders concepts such as "at rest" meaningless, this is in fact not the case in the FRW metric. Here comoving 

observers have a special status, and the CBR plays the role of the infamous "aether" which Michelson and 
Morley (1887) discredited with their famous speed-of-light experiment. 
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These deviations can be conveniently split into two categories: spectral fluctuations and 
spatial fluctuations. To study the former, we average the spectrum from all directions in 
the sky and investigate the extent to which it deviates from a Planck spectrum. To study 
the latter, we summarize the spectrum in each direction n by a single number, the best 
fit temperature T(n), and investigate how this temperature varies across the sky. In some 
instances, these two types of deviations are linked — we will return to this case when we 
discuss the so called SZ effect below. 

What sort of spectral distortions might we theoretically expect? One recurring scenario 
in this thesis is that the IGM was reheated and rcionizcd at some rcdshift z < 10^. As 
long as the IGM temperatures arc not too extreme, the thermal brcmsstrahlung and double 
Compton processes will remain frozen out, so that the number of CBR photons will not 
change. However, the photons will Thomson scatter off free electrons, which will modify 
their energy slightly and change their direction of motion radically. As we will see in the 
following two sections, the former effect tends to cause spectral distortions whereas the latter 
effect tends to suppress spatial fluctuations. 

2.5.3 Spectral distortions 

In this subsection, we will assume that the universe is completely uniform and isotropic. 
With this assumption, the time-evolution of the photon spectrum n^(r, k) in the presence 
of free electrons of temperature Tg is given by the Kompaneets equation (Kompaneets 1957) 



This equation describes the effect of Compton scattering only, and is valid after about 
z ^ when double Compton scattering and thermal Brcmsstrahlung have frozen out 
so that the number of photons is conserved. Here r and k are taken to be the comoving 
quantities, so that the redshift term disappears from the equation and we simply have ri-y = 
if He = 0. Direct substitution shows that n-y = also in the more interesting case where 



This means that if the photons start out with a Planck spectrum and if the electrons 
have the same temperature as the photons, then there will be no spectral distortions at 
all. Indeed, numerical integration of this nonlinear partial differential equation shows that 
for most physically reasonable initial spectra, the presence of the free electrons causes the 
photons to approach such a steady state. However, the resulting steady state typically has 
a "chemical potential" jJ- ^ 0. Thus if some process such as the decay of a mysterious relic 
particle creates non-thermal photons in the early universe, the result will be a spectrum 
with /i 7^ 0. The FIRAS experiment has placed the strong limit < 3.3 x 10~^ on the 
value of n (Mather et al. 1994). 




where 




1 



36 



CHAPTER 2. A COSMOLOGY PRIMER 



At the redshifts z < 10^ that we probe in this thesis, this convergence towards a steady 
state is no longer effective, being basically frozen out. Suppose the free electrons have a 
temperature Tg 7^ T-y. Substituting this into the right-hand-side of the Kompaneets equation 
and using the fact that the resulting spectral distortions will be quite small, we obtain the 
first order approximation that the spectral distortion today will be simply (Zel'dovich &; 
Sunyaev 1969) 



and the integral is to be taken from the early time when photon-creating processes froze out 
until today. Thus we see that to first order, the distortion will always have the same spectral 
profile, and that only the magnitude of the distortion depends on the thermal history of the 
electrons. This magnitude, referred to as the Compton y-parameter, thus summarizes the 
effect of free electrons by a single convenient number. Loosely speaking, the y-distortion 
indicates a relative excess of high-energy photons and a deficit of low-energy ones, which is 
caused by collisions in which hot electrons transfer energy to CBR photons. 

There are two types of y-distortions of current observational interest. The first type 
is the average y-distortion over all of the sky, typically caused by very early reionization. 
This is the topic of Chapter 5. To date, no such y-distortion has been observed, and the 
upper limit of y < 2.5 x 10~^ places interesting constraints on the thermal history of the 
IGM. The second type is the localized y-distortion caused by hot plasma in rich clusters 
of galaxies. This is usually referred to as the Sunyaev-Zel'dovich effect, or the SZ effect 
for short, after the Russians who invented the y-parameter (Zel'dovich &; Sunyaev 1969). 
The SZ-effect has indeed been observed in the direction of a number of large galaxy clusters 
— see for instance Birkinshaw Sz Hughes (1994). The experimental uncertainties are still 
large; around a factor of two. As soon as these shrink, the SZ-effect will become of great 
importance, since it provides an independent way of computing the distance to a galaxy 
cluster and thus can be used to determine the elusive Hubble constant. A recent review of 
the issues touched upon above is given by Hu &: Silk (1993). 

2.5.4 Spatial distortions 

In this subsection, we will take the opposite approach: we will assume that there are no 
spectral distortions, so that the CBR that we observe from the direction n in the sky has a 
Planck spectrum characterized by a temperature T^(n). It is the COBE observation of this 
function (Smoot et al. 1992) that has recently been plastered on front covers of magazines 
across the world and prompted various famous people to say silly things about holy grails 
and seeing the face of God. Our discussion of the ripples in T^(n), usually referred to 
simply as CBR fluctuations, will be organized as follows: first, we list some handy formulas 
for working with spherical harmonics. Then we review how the above-mentioned formulas 




for X ^ 1 
for X S> 1 



where 
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for random fields are altered when dealing with a sphere instead of R^. After this, we discuss 
how the CBR fluctuations are related to the fluctuations in the density field S. Finally, we 
show how early reionization can suppress CBR fluctuations, which is relevant to Chapters 
3, 4 and 5. A recent review of CBR fluctuations is given by White, Scott & Silk (1994). 

Working with spherical harmonics 

For functions that live on a sphere, the analogue of a Fourier expansion is an expansion in 
spherical harmonics. The spherical harmonics are defined as 

where Pl^ are the associated Legendre functions, and I and m are integers such that I > 
and \m\ < I. They have the symmetry property that 

where * denotes complex conjugation. These functions form a complete orthonormal set on 
the unit sphere. The orthonormality means that 

I Yii,{e,<p)Yir^,{e,<p)dn = 5u,5mm', 

and the completeness means that we can expand any function A as 

oo / 
1=0 m=-l 

where 

aim = j Yil{e,ip)A{e,ip)dn. 

Throughout this section, we use the differential solid angle notation 

dCl = sin 6 dOdip. 
Also, we will sometimes replace 6 and if by the unit vector 

fi = (sin 6 cos ip, sin 6 sin ip, cos 9) 

and write things like 

Yi^{e,f) = YUh). 

With this notation, the so called addition theorem for spherical harmonics states that 

' 27-1-1 

Y: Yr^ih)YUh') = ^Pi{h.h% 

m=—l 

where Pi = P^ are the Legendre polynomials. 

For the reader who wants more intuition about spherical harmonics, it is good to lump 
together all harmonics with the same I value. For / = 0, 1, 2, 3 and 4, these sets are 
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referred to as the monopole, the dipole, the quadrupole, the octupole and the hexadecapole, 
respectively. Geometrically, when one expands a function /(n) in spherical harmonics, 
I = picks up the constant part, I = 1 picks up the remaining linear part, I = 2 picks 
up the remaining quadratic part, I = 3 picks up the remaining cubic part, etc. In group 
theory jargon, the spherical harmonics corresponding to different Z-values are the irreducible 
representations of the group of rotations of the sphere. This means that if one expands a 
rotated version of the same function /, the new spherical harmonic coefficients aj^ will be 
some linear combination of the old ones ai^, 

I 

^'im' ~ Dijn'jn aim, 

m=—l 

such that different Z-values live separate lives and never mix. For instance, however one 
chooses to rotate a linear function (Z = 1), it will always remain linear and never say pick 
up quadratic terms. 

Random fields on the sphere 

Ty(n), the CBR temperature that we observe in the direction n in the sky, is modeled as a 
random field. It is more convenient to work with the dimensionless version 

(T^(n)) 

which is often denoted AT/T in the literature. This field is related to the random density 
field 6 that we discussed earlier on, and we will return to how they are relater further on. 
All the formulas and definitions we gave for 5 have spherical analogues for A, as we will now 
see. A is usually assumed to be Gaussian, which follows from the random phase assumption 

where the coefficients Ci constitute the spherical version of the power spectrum P{k). The 
angular correlation function is defined as 

c{e) = (A(n)A(n')), 

where n • n' = cos 9, and the right hand side is independent of the actual directions by 
the isotropy assumption. Just as the spatial correlation function was the Fourier transform 
of the spatial power spectrum, the angular correlation function is what might be called a 
"Legendre transform" of the angular power spectrum Q. Using the addition theorem, one 
readily obtains 

1 oo 

c{9) = —J2{21 + l)CiPiicose). 

1=0 

And just as the r.m.s. fluctuations of the smoothed density field was given by integrating 
the power spectrum against a window function w, the r.m.s. temperature fluctuations seen 
by a CBR experiment with some angular selection function (beam pattern) is given by a 
weighted average of the coefficients Ci, the weights being a sort of "Lcgcndrc transform" of 
the beam pattern. Explicit examples of this are given in Appendix C, and a few selected 
window functions are plotted Figure 2.4. 
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The monopole and the dipole 

The non-cosmologist reader may wonder why all sums over I in that chapter start with 
I = 2 rather than with 1 = 0. This is because we have no way of measuring the monopole 
coefficient cq or the dipole coefficient ci. 

If we knew the ensemble average (T-y) (or could measure it by repeating the COBE exper- 
iment in many different horizon volumes throughout the universe and invoking ergodicity), 
then we could calculate our observed monopole as the difference between the average Tj 
that we observe in our sky and {Tj). But we can't... In other words, our CBR monopole is 
the difference between the average temperature in our sky and the average temperature of 
the entire universe. But since we use the former as our estimate of the latter, this gets us 
nowhere with the monopole. Since direct estimates of the age of the universe arc still quite 
uncertain, but wc know the CBR temperature to three decimal places, the CBR in a sense 
merely tells us what time it is, i.e. how much our universe had expanded by the time this 
civilization of ours turned up on the scene. Thus it seems highly unlikely that wc will be 
able to use any kind of physics to predict (Tj) accurately (to one part in 10~^ or so) in the 
foreseeable future. 

The situation is similar with the dipole. If we knew that we were comoving observers, 
i.e. at rest with respect to the comoving FRW coordinates, then we could measure our local 
CBR dipole as the dipole we actually see in the sky. But the most accurate way we have 
to determine the comoving rest frame is to calculate the frame in which the CBR dipole 
vanishes. In other words, our CBR dipole is the difference between the average dipole in our 
sky and our velocity vector. Since we use the former as our estimate of the latter, this gets 
us nowhere with the dipole either. Thus the lowest multipole which gives us information 
about the power spectrum is Z = 2, the quadrupole. 

The relationship between 6 and A 

To date, the only really accurate way to compute the CBR power spectrum Q is to integrate 
the Boltzmann equation numerically. Today it takes a few hours of CPU time on a high-end 
workstation to compute the first few thousand Q-coefHcients. The results for a CDM model 
and a BDM model are plotted in Figure 2.5, multiplied by 1(1 + 1). To give some intuition 
about the general shape of these functions, we will very briefly mention the main physical 
effects at work. 

In the absence of rcionization, the CBR photons that reach us today are carrying infor- 
mation about the electron distribution on the last-scattering surface, the spherical region 
around us where the photons last Thomson scattered off of an electron some some 10^" 
lightyears away. There are three mechanisms through which adiabatic density fluctuations 
S on the last-scattering surface cause CBR anisotropics: 

1. Fluctuations in the gravitational potential cf) causing a gravitational redshift/blueshift 

2. Bulk velocities on the last-scattering surface, causing Doppler shifts 

3. Fluctuations in the photon density on the last-scattering surface 
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For adiabatic fluctuations, it turns out that the third eflFect cancels two thirds of the first 
effect. The net result in known as the Sachs- Wolfe effect, and for a matter power spectrum 
P{k) oc fc", the Sachs- Wolfe effect turns out to give 



which for the case n = 1 reduces to simply 



The Sachs- Wolfe effect is the dominant source of anisotropics on large angular scales 6^1°, 
corresponding to small I. This is why the n = 1 CDM curve in Figure 2.5 is flat for small I. 

The funny-looking bumps in the CDM curve, known as "Doppler peaks" are due to a 
combination of effects 2 and 3. Finally, the last scattering surface has a finite thickness, 
which means that when we look in a given direction, we are in fact seeing a mixture of 
photons emanating from different points along that line of sight. This averaging of the 
fiuctuations at many different points washes out fluctuations on very small angular scales, 
which is why all the curves approach zero for I 3> 10'^. 

2.5.5 How reionization suppresses fluctuations 

As mentioned in the introduction, reionization would affect the CBR in at least three ways: 

1. The hot electrons would cause a y-distortion, as discussed above. 

2. The fluctuations in the density and velocity of these hot electrons would generate new 
CBR anisotropics. For $7 = 1, due to various cancellations, it turns out that this so 
called Vishniac effect becomes important only on very small angular scales. We will 
not discuss it further here — see Hu et al. 1994 for a thorough discussion. 

3. Spatial fluctuations on angular scales below a few degrees could become suppressed, 
while fluctuations on larger scales such as those probed by CORE would remain fairly 



This third effect can be calculated numerically by integrating the linearized Boltzmann 
equation as mentioned above. Since hitting [RETURN] and waiting for hours does not 
greatly enhance ones intuitive understanding of the underlying physics, we will devote this 
section to deriving a simple analytic approximation which agrees fairly well with the nu- 
merical results for relatively late (z < 200) reionization. This is based on some ray tracing 
calculations made by the author. Examples of ray tracing work in the literature are Peebles 
(1987) and Anninos et al. (1991). 

The ray tracing approach 

Let nz{x,q) denote the six-dimcnsional phase-space density of photons at redshift z that 
have comoving position x and wave vector q = gq. If this function is known at our location 




unaffected. 
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(x = 0) today {z = 0), it is straightforward to predict the outcome of any microwave 
background experiment. The time evolution of this function between early times and today 
is described by the Boltzmann equation for radiative transfer. Since this equation is linear 
in n and invariant under spatial translations, we can write 



(x,q)= y"G(x-x',q,q',z,z')«^'(x',q')f^^a;'d'^9' (2-6) 



for some integral kernel G. One way to interpret this G is as the probability density that 
a photon that has position x and momentum q at redshift z was at x' with momentum 
q' at redshift z' . Thus an alternative to integrating the Boltzman equation numerically 
is to evaluate this probability distribution by Monte Carlo simulations of trajectories of 
large numbers of photons. As we will see, it is sometimes possible to obtain good analytic 
approximations for G as well. 

The geometrical smudging effect 

Reionization of the intergalactic medium (IGM) causes CBR photons to Thomson scatter 
off of free electrons, which affects CBR fluctuations in the three ways listed above. The 
third of these effects, to which we are limiting our attention, is purely geometrical in nature, 
i. e. independent of the fluctuations in the electron temperature and velocity. Thus ignoring 
the frequency change that photons experience when they scatter, the magnitude q simply 
stays constant for all photons (there redshift is taken care of by the fact that we are using 
comoving position and momentum coordinates). Hence we can write equation (2.6) as 

nzi^, ^,q) = j G{x - x', q, q', z, z')nz' (x, q', q) d^x'dO!, (2.7) 

where the angular integration is over all directions q'. Assuming Planck spectra and using 
the notation introduced in Section 2.5.4, this implies that 

A^(x, (i) = j G(x - x', q, q', z, /)Ay (x', q') d^x'd^'. (2.8) 
The isotropy approximation 

A very useful approximation is that the radiation is isotropic early on, at the last scattering 
epoch z = Zrec- In this approximation, we can write 

^2rec(x,q) = A^ec(x), 

since the left hand side is independent of q. For a thin last-scattering surface, it is easy to 
see that this is the case for the contribution from both the Sachs- Wolfe and the intrinsic 
density fluctuations, whereas the Doppler term will have a dipole anisotropy. With the 
isotropy approximation, we have 

A^(X, q) = y G(X - x', q, Z, Zrec)Kec{x')d^x', (2.9) 

where we have defined the angularly averaged kernel 

G(Ax, q, z, z') = J G(Ax, q, q', z, z')dQ! . 
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Prom here on, we will suppress z and z' , taking z = and z' = Zrec- Fourier transforming 
equation (2.9) with respect to x and using the convolution theorem leaves us with 

Ao(k,q) = G(k,q)A,ec(k). (2.10) 

Thus there is no mode coupling at all, and for fixed q, the transfer function is simply the 
Fourier transform of the probability distribution G( Ax) = G( Ax, q) . 

The transfer function G 

Before we turn to the problem of explicitly evaluating G, we will briefly discuss its interpre- 
tation. Assuming that A is a Gaussian random field, it is easy to show that 

(Ao(k,q)Ao(k',q)) = |G(k, q)|2(A(k),ecA.ec(k')), 

i.e. that the quantity |G(k, q)| plays the role of a transfer function. Thus the radiation 
power spectrum on scales A = 27r/A; is suppressed by a factor |G(k, q)p. Let us make a few 
elementary observations that are valid for an arbitrary probability distribution G. 

Observation (I): |G(k)| < 1, with equality if k = or G(Ax) = (5(Ax - Axq). 

Observation (II): G(k) — >^ as |k| — >^ oo if G is an integrable function (as opposed to say 
a tempered distribution like (5(Ax — Axq)). 

Observation (III): If we define a to be the mean and S to be the covariance matrix of 
the probability distribution G, then 

- . 1 3 

G(k) = 1 — iajnkjn — -^{O'mO'n + Smn)krakn + 0(|k| ), 

SO 

|G(k)|2 = 1 - S^nkmkn + 0(|k|3). (2.11) 

(Repeated indices are to be summed over, from 1 to 3.) 

Observation (I) tells us that Thomson scattering causes a low-pass filtering of the Fourier 
modes, i.e. leaves very long wavelengthes unaffected. G(Ax) = 5(Ax — Axq) corresponds 
to no reionization, in which case there is no smudging at all and no modes are suppressed. 
Observation (II), which is known as Riemann-Lebesgue's Lemma, tells us that unless there 
is a finite probability for no scattering at all, very high frequency modes get almost entirely 
damped out. 

Observation (III) gives us a Taylor expansion of G(k) around the origin, which will prove 
useful below. 

All quantities of physical interest (G^, C{9), etc.) can be computed from the initial power 
spectrum at Zrec once the transfer function G(k, q) is known, and explicit formulas for this 
are given in most standard texts. A crTidc rule of thumb simply identifies the comoving length 
scale A with the angle that this distance subtends on the last scattering surface, a spherical 
region whose radius is our roughly our horizon radius uoTq k, 6000/i~^Mpc. Another rule 
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of thumb is to identify a multipole moment £ with the angular scale 9 180° /Z. Thus for 
small angles, we have the rough correspondence 

X (X 6 (X - (X —, 
I k 

where the proportionality constants are given by the correspondence 

6* « 1° ~ Z « 200 ~ A « 100/i~^Mpc. 
Let us define the smudging scale 

Xc= (det5)^/^ 

i.e. as the geometric mean of the three eigenvalues of S^^^. For realistic scenarios, the 
smudging tends to be fairly isotropic. Thus all three eigenvalues of the covariance matrix S 
tend to be of the same order of magnitude, and the Taylor expansion (2.11) gives 

|G(k)|2 « 1 - (Ac/fc)2 (2.12) 

fov k <^ . When we detect a CBR photon arriving from some direction q in the sky and 
ask where it was at 2; = Zrec, the smudging scale Ac is roughly the standard deviation of the 
answer (see Figure 2.7). 

Approximating the transfer function G 

For fairly late (say z < 200) reionization, there is a substantial probability that a given 
CBR photon manages to dodge all the reionized electrons without being scattered a single 
time. Writing this probability as e~'^* , where the quantity Tt is called the optical depth and 
will recur in subsequent chapters, the probability that a photon is scattered n times will be 
e^'^'-p^, a Poisson distribution, (rt is not to be confused with the conformal time r.) Let us 
expand the propagator G in a type of Born expansion 

G=Y.e-^^^Gn, (2.13) 

where G„ is the probability distribution corresponding to the case where the photon scatters 
exactly n times. Omitting the trivial dependence on q and q', we have simply 

G'o(Ax) = 5(Ax-Axo), 

where Axq is some constant, since if the photon never scattered, we know exactly where it 
came from. Gi can be readily calculated analytically, which we will do in the next section, 
after which the higher order terms can be calculated by repeated integrations if desired. 
A nice feature is that for late reionization, Tt is often so small that a good approximation 
is obtained by dropping all but a few of the first terms in equation (2.13). For instance, 
reionization around z = 50 with = 1 and hVii, = 0.03 gives e^"^* ~ 0.70, e^'^^Tt ~ 0.25 and 

2 

e~'^*-^ ~ 0.044, so if we keep only the Gq and Gi terms, we will be off by only 5%, and if we 
include the G2 term as well, we are off by less than 1%. Physically, this simply means that 
if it is not that likely that a photon scattered even a single time, then it is quite unlikely 
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that it scattered more than a few times, so that we can safely neglect that possibility. If one 
nonetheless wants to include all terms, it is readily done by Monte Carlo simulation that 
follow individual photons as they propagate upwards in Figure 2.1, occasionally changing 
direction as they Thomson scatter off of free electrons. Some examples of this are shown in 
Figures 2.6 and 2.7. Let us write the infinite sum (2.13) as 

G = e-^'Go + (1 - e-"*) 

where Gi+ contains the contributions from all terms with n> 1. We Fourier transform this 
with as before: 

G(k,q)=e--^ + (l-e--=)Gi+(k,q) 

Since Gi+ contains no delta-functions, its Fourier transform approaches zero as |k| — ^ oo 
according to Observation I. Thus 

G(k, q) e""^" as k oo. 

Combining this with the estimate (2.11), we thus know how the transfer function behaves 
in both limits: the approximation 

|G(k,q)| (1 - e"^*) e-^^'^^^ + e"^* (2.14) 

becomes exact both as — > and as k ^ oo. Using equation (2.12) yields the isotropic 
approximation 

|G(k,q)| « (1 - e"^*) e-(^=^)'/2 + e-^*. (2.15) 

This function is plotted in Figure 2.8 for a few reionization scenarios. The values of Ac used 
in this plot where computed through Monte Carlo simulation in conjunction with Wayne 
Hu. 

The backward light cone 

We will conclude our discussion of how reionization suppresses fluctuations by computing 
the function Gi analytically. Let us assume that J7 = 1. To simplify the calculations that 
follow, we will use comoving coordinates x and conformal time r, so that all light rays (null 
geodesies) make 45 degree angles with the r axis in (r, x)-space, just as in Penrose diagrams. 
Let us choose the scale factor ao such that r = 1 corresponds to today, i.e. so that 

1 

The situation is illustrated in Figure 2.1, where we have suppressed one of the three spatial 
coordinates for the sake of the plot. Wc arc at the apex of the cone, our backward light cone, 
and in the absence of reionization, the CBR photons that reach us today have been moving 
on the surface of this cone ever since recombination. 
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Where where the photons at recombination? 



If we detect a CBR photon arriving from a given direction in the sky, where was is at a 
redshift z7 The presence of free elections makes it possible for a photon to change direction a 
number of times by inelastic Thomson scattering. We neglect other ways by which photons 
can change direction, such as gravitational lensing. 

To avoid double use of the letter z, we will write comoving coordinates as x = (n, v, w). 
Suppose a photon arrives from the lu-direction to the point x = at some conformal time 
r. Where was it at the earlier time ti = r — Ar? If it has not been scattered during the 
interval Ar, it has followed a null geodesic, and its position at ri is given by 

x(°)(ti) = (0,0,At). (2.16) 

If it has been scattered exactly once during the interval, say at the conformal time ti + Tg, 
and previously came from the direction (9, (j)) in spherical coordinates, then 



{u, v,w) = [ts sin 9 cos 0, Tg sin 9 sin cp, Tg cos 9 + (Ar — Tg 
For |x| < Ar, this is readily inverted, yielding 

's ~ 2(Ar-w)) 

? = 2arctan^^-^' 



(2.17) 



(2.18) 



b = arctan - 

u 

To answer the question "where did the photon come from if it scattered exactly once" , 
let us take Tg 9 and cj) to be independent random variables with probability distributions 
fr, fe and f(j), and calculate the probability distribution for the random variable x(ri) as 
defined by equation (2.17). This probability distribution is 



Gi(x) = 

where the volume element is 

d{Tg,9,cj)) 
d{u, V, w) 



/.[r.(x)]/,[0(x)]/^[<^(x)]^ 




d{Ts,e,<t>) 

to) 



if |x| < Ar, 
if |x| > Ar, 



(2.19) 
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For Thomson scattering, the angular dependence is = (1 + cos^^), so 



_ 3 



(l + cos^ 9) sine. 



The function depends on the ionization history of the universe, so we will leave it arbitrary 
for the time being. Figure 2.1 shows a situation where the photon was scattered at 2; = 10: 
we then know that it must have arrived on the smaller backward light cone emanating from 
the scattering event. This corresponds to the probability distribution /r(T) being a delta 
function. If we know that the photon scattered sometime around z ~ 10, and exactly once, 
but do not know exactly when, we obtain a probability distribution Gi(x) that is smeared 
out around the perimeter of this z = 10 backward light cone. For more physically realistic 
cases, where scattering occurs predominantly at higher redshifts, the resulting probability 
distribution tends to be unimodal. 
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2.6 Cosmological Chemistry 

Almost all of the following chapters will involve the ehemistry of the IGM, the intergalaetic 
medium, and frequent references will be made to constraints from nucleosynthesis. The 
purpose of this section is to convince the reader that it is quite straightforward to understand 
chemical equations, even if one has never been anywhere near a test tube. 

2.6.1 Evolution of the IGM 

As an example, let us take the evolution of the ionization state of the IGM in the simplest 
case, neglecting helium and ionizing photons. Let x denote the ionization fraction and n the 
total number density of protons, free and bound. Then the densities of hydrogen atoms, free 
electrons and free protons are nn = (1 ~ x)n, Ue = xn and Up = xn, respectively. Thermal 
ionization occurs through the reaction 

e~ + H e~ + e~ + p'^ , 

i.e. the electron is a catalyst that collides with the hydrogen atom and uses its kinetic 
energy to ionize it. Thus the rate per unit volume for this reaction is proportional to 
nef^H = n'^x{l — x). Similarly, the rate for the recombination reaction 

e~ + p"*" I— > + photon 
is proportional to rierip = r?x^ . Thus the time-evolution of x is 

X = n\Ax{\ - x) - Bx\ (2.20) 

where the proportionality constants A and B are functions of temperature (since the particles 
come near each other more often if they are moving faster, and since the reaction cross 
sections depend on the kinetic energies of the particles). The cosmologist generally takes 
such functions from the atomic physics literature without worrying to much about how 
they were calculated. In addition to the equation specifying i, one also needs an equation 
for T, to be able to compute the evolution of temperature. Such an equation is quite 
straightforward to interpret, the various terms simply corresponding to the various cooling 
and heating effects that one wishes to take into account. 

Before integrating the above-mentioned equations numerically, it is often convenient to 
change the independent variable from time to redshift. Then the IGM density is simply 
given by n = no(l + 2;)^, no denoting the density today. Using equation (2.1) for dz/dt., 
equation (2.20) becomes 

-— = HQnQ^====[Ax{l -x)- Bx% 

2.6.2 Nucleosynthesis 

For a more complicated problem like nucleosynthesis the formation of light elements 

during the first three minutes after the Big Bang, the mathematical structure of the problem 

^^This calculation would usually be referred to as nuclear physics rather than chemistry, but the resulting 
equations arc of the same type the only difference being that nucleons (protons and neutrons) rather than 
atoms are the fundamental entities. 
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is the same, merely more complicated. Instead of one single fraction x to keep track of, there 
will be a whole bunch of fractions Xj, i = 1, 2, n: the fraction of protons in hydrogen, the 
fraction of protons in helium-3, the fraction of protons in helium-4, the fraction of protons 
in lithium, etc. The time derivative of each one of them is given by an equation like the 
one above, where the left hand side is Xj and the right hand side contains one term for 
each process that destroys or creates species i. The resulting system of n coupled first-order 
ordinary differential equations is fed into a computer and solved numerically. The numerical 
part is somewhat complicated by the fact that chemical equations tend to be stiff, but there 
is a plethora of standard packages that handle this sort of problems, for instance in the 
NAG library. 

When this calculation is carried out — see Malaney & Mathews (1993) for a recent re- 
view, results are obtained that are in fairly good agreement with the primordial abundances 
of helium-4, helium-3, deuterium and lithium that are inferred from observational data. 
Pioneered by George Gamow and others in the forties (Alphcr, Bethe Sz Gamow 1948^^), 
nucleosynthesis predictions are one of the most striking successes of the Big Bang model. 
Apart from being a powerful weapon against rival cosmological theories, nucleosynthesis 
calculations also place an important constraint on the standard cosmological model, since 
the results turn out to depend strongly on the ratio between the densities of baryons and 
microwave background photons. This translates into a strong dependence on h^Qi,. The 
most recent calculations to date (Smith et al. 1993; Walker et al. 1991) give the 95% 
confidence interval 

0.010 < h^^lb < 0.015. 

(Levels of confidence must be taken with a grain of salt in cosmology, where systematic errors 
often dominate.) This interval keeps shifting around and getting narrower as observational 
abundance data gets better and the the knowledge of the nuclear cross sections that go into 
the theoretical calculation improves. 

2.7 Astronomy Jargon 

Since the scales involved in cosmology are many orders of magnitude larger than those 
usually encountered in physics, cosmologist tend to use larger units of measurement. A 
convenient unit of time is the Giga-year: 

IGyr = 10^ years « 3.15 x lO^S « IO^^tt s. 

An astronomical unit, or AU for short, is the average distance between the earth and the 
sun. 

lAU ^ 1.4960 X lO^^cm. 

For historical reasons, related to the parallax triangulation method by which the distance 
to nearby stars is calculated, the standard astronomical distance unit is the parsec. This is 

Gamow added his friend Hans Bethe to the author hst of this so called a/37 paper as a joke... 
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the distance at which an AU makes an angle of one arcsecond (1/3600 of a degree) in the 
sky. Thus 

Ipc = ^ ^^^^ AU K 3.086 X 10^^ cm ^ 3.2615 hght-years. 

This is still a tiny distance by cosmology standards. The standard length unit in cosmology 
is the megaparsec: 

IMpc 3.086 X lO^^cm. 

The standard astrophysical unit of mass is solar mass, denoted AIq. The mass of the sun is 

Mq 1.989 X lO^^g. 

Apart from involving many powers of ten, some cosmological parameters arc quite un- 
certain. It is convenient to write such parameters as the product of a large round number 
that is in the right ballpark and a small dimensionless "fudge factor" . For instance, the di- 
mensionless quantity h is ubiquitous in the cosmology literature. This is because the Hubble 
constant is factored as 

V Mpc J 

Whereas Hubble originally estimated that ~ 5, most cosmologists today assume that 
0.5 < h < 0.8. The jury is still out. If some form of matter x has a uniform density px, then 
this density is usually factored as 

Px ~ Pcrit^xi 

where ^ 

Peru = 1^ ~ 1-9 X IQ-^^h^ ^ 



SttG ■ cm3 
is the critical density mentioned above. Thus is a measure of the baryon density, ^igm 
is a measure of the IGM density, etc. 

The only other unexplained astronomy jargon that I have been able to spot in the 
subsequent chapters is some reference to different types of stars. For historical reasons (ever 
heard that excuse before...?), stars are classified in categories called O, B, A, F, G, K, and 
M, where 0-stars arc the hottest and M-stars the coolest. Within each category, there are 
ten sub-categories specified by numbers, where for instance an 03 star is hotter than an 09 
star. 
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2.8 An Abstract for Non-Cosmologists 

So in plain English, what is this thesis all about? 

It is now widely accepted that something heated and ionized the hydrogen that fills 
intergalactic space, but we don't know quite when this happened. 

If this is to have had interesting efi'ects on the microwave background fluctuations, then 
it must have occurred quite long ago, at redshifts greater than 50 or so. This is the topic of 
Chapters 3 and 4. 

Some theories predict that the hydrogen has indeed been ionized pretty much all along. 
In Chapter 5 we show to what extent such scenarios are ruled out by recent measurements 
of the microwave background spectrum. 

To affect the microwave background, 100% ionization is not required — some reasonable 
fraction like 30% is quite sufficient, providing that the ionization happens early enough. By 
studying the spectra from quasars, however, one gets more extreme constraints: in "recent" 
times, for redshifts less than about five, it appears that at least 99.999% of the hydrogen 
atoms must have been ionized. It is still not entirely clear whether photoionization alone is 
able to produce the high temperatures required to achieve this. In Chapter 6 we provide an 
alternative model for how this may have happened. 

Finally, there is increasing evidence that the primordial power spectrum may be a more 
complicated function than a simple power law, an assumption that has been used in much of 
the literature, including Chapters 3, 4 and 6. Motivated by this evidence, the final chapter 
is devoted to developing a formalism that allows cosmological constraints to be placed that 
are completely independent of the shape of power spectrum. Some supplementary material 
for Chapters 3 and 7 is given in the appendices. 
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Figure 2.1: Our backward light cone. 

Our backward light cone in a fiat CI = 1 universe is shown in comoving 
spatial coordinates, with conformal time on the vertical axis (one 
of the three spatial dimensions has been suppressed). As can be 
seen, in these coordinates, light rays propagate in 45° lines just as 
in Euclidean space. The horizontal circles are labeled with their 
corresponding redshifts, the dotted circle corresponding to the Big 
Bang. If we (at the apex of the cone) detect a CBR photon arriving 
from the right in the picture, it could have been anywhere in the 
shaded region aX z = 10^, the recombination epoch. If we know that 
it has not been scattered since z = 10, it must have come from the 
double-hatched region. 
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Figure 2.2: Evolution of a top hat overdensity. 

The time-evolution of the density in a spherical ( "top hat" ) overdense 
region in a flat universe is plotted as a function of redshift z. Prom top 
to bottom, the three solid curves correspond to virialization redshifts 
Zmr of 200, 100 and 50, respectively. The dashed curve gives the 
density of the IGM. The number density is given in nucleons per 
unit volume, i.e. n = p/nip. 
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Figure 2.3: Assorted power spectra. 

The power spectrum P{k) is plotted for a CDM model and two vari- 
ations on the theme. The solid curve is a standard CDM power 
spectrum with h = 0.5, = 1, and n = 1 (Bond & Efstathiou 
1984). The dotted curve has a lower "shape parameter" hQo = 0.25, 
which can be attained in for instance the MDM and rCDM scenar- 
ios discussed in subsequent chapters, and corresponds to sliding the 
CDM curve to the left. The dashed curve has a lower spectral index 
n = 0.7, the so called tilted model, which approximately corresponds 
to rotating the CDM curve clockwise. All three power spectra have 
been normalized to COBE. 
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Figure 2.4: Assorted angular windows functions. 

The angular window functions Wi are plotted for a number of re- 
cent CBR experiments. The data is courtesy of White, Scott & Silk 
(1994). 
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Figure 2.5: Assorted angular power spectra. 



The angular power spectrum Q is plotted for a CDM and a BDM 
models. The CDM model has ^^o = 1, = 0.06, h = 0.5 and n = 1. 
The BDM model has ^0 = 9,^ = 0.15, h = 0.8 and n = -0.5. Both 
are normalized to COBE, i.e. give the same value when integrated 
against the window function for the COBE DMR 10° pixel variance. 
The models are courtesy of Naoshi Sugiyama and Wayne Hu. 
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Figure 2.6: Monte Carlo photons, top view. 



The trajectories of seven photons in a Monte Carlo simulation of a 
fully ionized universe are shown in comoving coordinates {u,v,w), 
with the ii-coordinate perpendicular to the page. (In terms of Fig- 
ure 2.1, this is a view from "above".) The question being asked is 
where photons arriving from the w-direction were in the past, so one 
can just as well model the photons as emanating from the apex of the 
light cone and traveling backwards in time. With this terminology, 
they all enter this figure from below, and part ways as they Thomson 
scatter. Notice how their mean free path gets shorter and shorter, 
which is due to the increase in electron density at higher redshifts. 
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Figure 2.7: Monte Carlo photons, side view. 

The trajectories of ten photons in the Monte Carlo simulation de- 
scribed in Figure 2.6 are plotted "from the side" with reference to 
the light cone in Figure 2.1. The horizontal axis is {u^ + v'^ + itj^)^/^, 
the comoving distance from "here" (the apex of the light cone). The 
curves all end at r « 0.3, corresponding to z ^ 10^, i.e. standard 
recombination. The horizontal spread of these endpoints roughly 
corresponds to the smudging scale Ac- 
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Figure 2.8: The transfer function |G|. 

The transfer function i.e. the factor by which the r.m.s. CBR 
fluctuations are suppressed by reionization, is plotted as a function 
of comoving wavenumber for scenarios where the universe becomes 
reionized at some redshift Zjon- As A; ^ oo, the suppression factor 
simply approaches e""^*. For these scenarios, hQigm = 0.03, so Tt ~ 
(zjon/QS)^/^. The cutoff scales Ac where determined by Monte Carlo 
simulation, and turn out to be depend only quite weakly on the 
ionization redshift Zion- In all cases, degree-scale anisotropies are 
seen to be suppressed whereas fluctuations on the scale probed by 
COBE are almost unaffected. 



CHAPTER 2. A COSMOLOGY PRIMER 



Chapter 3 

Early Reionization and CBR 
Fluctuations 

In this chapter, early photoionization of the intergalactic medium is discussed in a nearly 
model-independent way, in order to investigate whether early structures corresponding to 

rare Gaussian peaks in a CDM model can photoionize the intergalactic medium sufficiently 
early to appreciably smooth out the microwave background fluctuations. We conclude that 
this is indeed possible for a broad range of CDM normalizations and is almost inevitable 
for unbiased CDM, provided that the bulk of these early structures are quite small, no 
more massive than about 10*^ Mq. Typical parameter values predict that reionization occurs 
around z = 50, thereby suppressing fluctuations on degree scales while leaving the larger 
angular scales probed by COBE relatively unaffected. However, for non-standard CDM, 
incorporating mixed dark matter, vacuum density or a tilted primordial power spectrum, 
early reionization plays no significant role. 

3.1 Introduction 

The first quantitative predictions of cosmic microwave background anisotropies in cold dark 
matter (CDM)-dominated cosmological models recognized that reionization by rare, early- 
forming objects could play a role in suppressing temperature fluctuations on small angular 
scales (Bond & Efstathiou 1984; Vittorio & Silk 1984). Now that the COBE DMR ex- 
periment has detected fluctuations on large angular scales (Smoot et al. 1992) at a level 
(within a factor of two) comparable to that predicted by CDM models, it is especially rele- 
vant to examine whether reionization can affect the degree scale searches that are currently 
underway. 

Cold dark matter models are generally characterized by a late epoch of galaxy formation. 
However, the smallest and oldest objects first go nonlinear at relatively large redshift. In 
this chapter we investigate, for a wide range of CDM normalizations, power spectra and 
efficiency parameters, whether reionization associated with energy injection by early forming 
dwarf galaxies can reionize the universe sufficiently early to smooth out primordial CBR 
temperature fluctuations. 

Although we go into some detail in Appendix A to make estimates of a certain effi- 
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ciency parameter, our overall treatment is fairly model-independent, and can be used as a 
framework within which to compare various photoionization scenarios. Our basic picture 
is roughly the following: An ever larger fraction of the baryons in the universe falls into 
nonlinear structures and forms galaxies. A certain fraction of these baryons form stars or 
quasars which emit ultraviolet radiation, and some of this radiation escapes into the ambi- 
ent intergalactic medium (IGM) and ends up photoionizing and heating it. Due to cooling 
losses and recombinations, the net number of ionizations per UV photon is generally less 
than unity. 

Apart from photoionization, early galaxies can also ionize the IGM through supernova 
driven winds, an ionization mechanism that will not be treated in this chapter. Although 
such winds can ionize the IGM by z = 5, early enough to explain the absence of a Gunn- 
Peterson effect (Tegmark et al. 1993), the relatively low velocities of such winds makes 
them unable to distribute the released energy throughout space at redshifts early enough 
(by z K, 50) to measurably affect the CBR. 

Our approach will be to first write the ionization fraction of the IGM as a product of a 
number of factors, and then discuss the value of each of these factors in more detail. Let us 
write 

X ~ fsfuvppfionj (^•-'-) 

where 

X = fraction of IGM that is ionized, 
fs = fraction of baryons in nonlinear structures, 

fuvpp = UV photons emitted into IGM per proton in nonlinear structures, 
fion = net ionizations per emitted UV photon. 

Let us first consider the case where the UV photons are produced by stars, and return to 
the quasar case later. Using the fact that a fraction 0.0073 of the rest mass is released in 
stellar burning of hydrogen to helium, we obtain 

13 6eV / l^'^) 

where 

fn = mass fraction hydrogen in IGM, 
fburn = mass fraction of hydrogen burnt, 
fuv = fraction of energy released as UV photons, 
fesc = fraction of UV photons that escape from galaxy. 

We will take the primordial mass fraction of helium to be 24%, i.e. fn = 76%. Now define 

the net efficiency 

fnet ~ fburn fuv fesc fion-i 

and equation (3.1) becomes 

X ~ 3.8 X 10^ fnet fs. (3.3) 

The key feature to note about this expression is that since 3.8 x 10^ is such a large number, 
quite modest efficiencies fnet still allow x to become of order unity as soon as a very small 
fraction of the baryons are in galaxies. As will be seen in the next section, this means that 
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reionization is possible even at redshifts far out in the Gaussian tail of the distribution of 
formation redshifts, at epochs long before those when the bulk of the baryons go nonlinear. 
This appears to have been first pointed out by Couchman and Rees (1986). 

3.2 The Mass Fraction in Galaxies 

In this section, we will discuss the parameter fg. Assuming the standard PS theory of 
structure formation (Press &; Schechter 1974), the fraction of all mass that has formed 
gravitationally bound objects of total (baryonic and non-baryonic) mass greater than M at 
redshift z is the integral of the Gaussian tail, 



fs = erfc 



Sc 



(3.4) 



y2a{M, z) 

where the complementary error function erfc(a;) = 2tt~^^^ e~^^du and a{M,z) is the 
r.m.s. mass fluctuation in a sphere containing an expected mass M at redshift z. is given 
by top-hat filtering of the power spectrum as 

"sin/cro cos /erg "'^ 



/■oo 

a{M,zf oc / P{k) 
Jo 



dk, (3.5) 



where P{k) is the power spectrum at redshift z and ro is given by ^Trrgp = M, p = 
being the density of the universe at redshift z. Although this approach has been criticized 
as too simplistic, numerical simulations (Efstathiou et al. 1988; Efstathiou & Rees 1988; 
Carlberg &: Couchman 1989) have shown that it describes the mass distribution of newly 
formed structures remarkably well. Making the standard assumption of a Gaussian density 
field, Blanchard et al. (1992) have argued that it is an accurate description at least in the 
low mass limit. Since wc arc mainly interested in extremely low masses such as 10^ Mq, it 
appears to suffice for our purposes. 

For our middle-of-the road estimate, we choose Sc = 1.69, which is the linearly ex- 
trapolated over density at which a spherically symmetric perturbation has collapsed into a 
virialized object (Gott &; Rees 1975; Efstathiou et al. 1988; Brainerd &; Villumsen 1992). 
We take Sc = 1.44 (Carlberg &; Couchman 1989) for the optimistic estimate, although 
the even lower value Sc = 1.33 has been discussed (Efstathiou & Rees 1988), and Sc = 2.00 
(Gelb &; Bertschinger 1992) for the pessimistic estimate. (Here and throughout this chapter, 
parameter choices are referred to as optimistic if they permit earlier reionization.) 

The fact that cr(M, z) — oo as M — implies that if we consider arbitrarily small scales, 
then all dark matter is in non-linear structures. Thus if no forces other than gravity were at 
work, so that the baryons always followed the dark matter, we would simply have = 1 at 
all z. However, it is commonly believed that galaxies correspond only to objects that are able 
to cool (and fragment into stars) in a dynamical time or a Hubble time (Binney 1977; Rees & 
Ostriker 1977; Silk 1977; White & Rees 1978). The former apphes to ellipticals and bulges, 
the latter to disks. Let us define the virialization redshift (1 + z^r) = {V^/Sc)o'{Mc,0), 
where Mc is some characteristic cutoff mass which is the total mass (baryonic and dark) of 
the first galaxies to form. Zyir is roughly the redshift at which the bulk of all baryons goes 
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non-linear. Using equation (3.4) and the fact that a{M,z) = (t(M, 0)/(l + z) in the Unear 
regime of CDM, we thus have 



fs = erfc 



l + z 



(3.6) 



1 -\- Zui^ _ 

A common assumption is that Mc ~ 10® M0, roughly the Jeans mass at recombination. 
Blanchard et al. (1992) examine the interplay between cooling and gravitational collapse in 
considerable detail, and conclude that the first galaxies to form have masses in the range 
1O^M0 to 10^ Mq, their redshift distribution still being given by equation (3.6), whereas 
Couchman & Rees (1986) argue that the first galaxies to form may have had masses as low 
as IO^Mq. 

As our CDM power spectrum today, we will use that given by BBKS (Bardeen et al. 
1986) and an n = 1 Harrison-Zel'dovich primordial spectrum: 

Pik)o.( .rMn(l + 2.34<,) 



V [1 + 3.89g + (16.1^)2 + (5.46g)3 + {Q.llqY]^/^ , 

where q = k/[h?Q.QM.\)C~^]. Throughout this chapter, we will take ilo = 1- 
Evaluating the cr^-integral in equation (3.5) numerically yields 

(t(1O^M0,O) « 33.76-^ 

for h = 0.8 and 

(t(10^Mo,0) «i 13.66"^ 

for h = 0.5, where the so called bias factor b = (7(8/t~^Mpc, 0)~^ has been estimated to lie 
between 0.8 (Smoot et al. 1992) and 2.5 (Bardeen et al. 1986). Our pessimistic, middle-of- 
the-road and optimistic CDM estimates of Zyir are given in Table 3.1, and the dependence 
of Zyir on Mc is plotted in Figure 3.1. This figure also shows three alternative models of 
structure formation: CDM with cosmological constant (Efstathiou et al. 1992); tilted CDM 
(Cen et al. 1992) and MDM, mixed hot and dark matter (Shafi & Stecker 1984; Schaefer 
&; Shafi 1992; Davis et al. 1992; Klypin et al. 1993). For the model with cosmological 
constant, we have taken a flat universe with h = 0.5, = 0.4 and Aq = 0.6. For the 
tilted model, the power spectrum P{k) is simply multiplied by a factor A;"^^, where wc have 
taken n = 0.7. For the tilted case, equation (3.6) still applies. For the MDM case, however, 
perturbations in the cold component grow slower than linearly with the scale factor (l + z)~^ 
and equation (3.6) is not valid. For the low masses we are considering, we have (Bond &: 
Szalay 1983) 

(^mdm[Mc,z)^ — — — — — , where 
(1 + z)°' 

1 



""4 



^25 - 24nHDM - 1 



Using the parameters from Davis et al. (1992), who take Q = 1 and ^hdm = 0.3, the MDM 
version of equation (3.6) becomes 
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Mixed 


Tilted 


Lambda 


Pess. 


Mid. 


Opt. 


Me 


75 


71 

10'' Mo 


71 

IO'^Mq 


n 

1O«M0 


71 

lO^'Mo 


10^ Mo 


Model 


MDM 


Tilted 


Lambda 


CDM 


CDM 


CDM 


h 


0.5 


0.5 


0.5 


0.5 


0.5 


0.8 


b 


1 


1 


1 


2 


1 


0.8 




1.69 


1.69 


1.69 


2.00 


1.69 


1.44 




2.9 


10.1 


8.4 


4.8 


17.2 


41.4 



Table 3.1: Galaxy formation assumptions 





Pess. 


Mid. 


Opt. 


fburn 


0.2% 


1% 


25% 


fesc 


10% 


20% 


50% 


fuv 


5% 


25% 


50% 


fion 


10% 


40% 


95% 


fnet 


1 X 10-*^ 


2 X 10"^ 


6 X 10"^ 


fuvpp 


4 


190 


24,000 



Table 3.2: Efficiency parameters used 



where a ^ 0.8 and we redefine z^r by 

(1 + ^iiir) 



V2a(Me,0) 



1/a 



with a{Mc, 0) referring to a pure CDM power spectrum. 

For the A case, perturbations grow approximately linearly until the universe becomes 
vacuum dominated at z ~ ilg ^ — 1 = 1.5, after which their growth slowly grinds to a halt. A 
numerical integration of the Friedmann equation and the equation for perturbation growth 
using h = 0.5, fio = 0.4 and Aq = 0.6 gives 

for z » 3. 

Since our A-modcl yields a value of Zyir very similar to our tilted model, we will omit 
the former from future plots. 



3.3 Efficiency Parameters 

In this section, we will discuss the various parameters that give fnet when multiplied to- 
gether. The conclusions are summarized in Table 3.2. 

/burn) the fraction of galactic hydrogen that is burnt into helium during the early life of 
the galaxy (within a small fraction of a Hubble after formation), is essentially the galactic 
metallicity after the first wave of star formation. Thus it is the product of the fraction 
of the hydrogen that forms stars and the average metallicity per star (weighted by mass). 
This depends on the stellar mass function, the galactic star formation rate and the final 
mctallicitics of the high-mass stars. For our middle-of-the-road estimate, we follow Miralda- 
Escude &; Ostriker (1990) in taking fbum = 1%; half the solar value. An upper limit to fbum 
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UV source 



Spectrum ^(1^) 



fuv 




03 star 
06 star 
09 star 



T = 50, OOOK Planck 
T = 40, OOOK Planck 
T = 30, OOOK Planck 
T = 50, OOOK Vacca 

a = I power law 



0.57 
0.41 
0.21 
0.56 



17.3 eV 
16.6 eV 
15.9 eV 

18.4 eV 
18.4 eV 
17.2 eV 
20.9 eV 
19.9 eV 



28,300K 
23,400K 
18,000K 
36,900K 
37,400K 
27,800K 
56,300K 
49,000K 



2.9 
3.4 
3.9 
2.2 
1.7 
2.7 
0.6 
1.6 



Pop. Ill star 
Black hole, QSO 



? 



a = 2 power law 
cc = power law 
T = 100, OOOK Planck 



? 



? 



0.89 



Table 3.3: Spectral parameters 



is obtained from the extreme scenario where all the baryons in the galaxy form very massive 
and short-lived stars with M ^ 3OM0 , whose metallicity could get as high as 25% (Woosley 
&; Weaver 1986). Although perhaps unrealistic, this is not ruled out by the apparent absence 
of stars with such metallicities today, since stars that massive would be expected to collapse 
into black holes. 

In estimating fesc, the fraction of the UV photons that despite gas and dust manage to 
escape from the galaxy where they are created, we follow Miralda-Escude Sz Ostriker (1990). 

For fuy, the fraction of the released energy that is radiated above the Lyman limit, we 
also follow Miralda-Escude & Ostriker (1990). The upper limit refers to the extreme SOM© 
scenario mentioned above. For reference, the values of /„„ for stars with various spectra 
arc given in Tabic 3.3, together with some other spectral parameters that will be defined 
and used in Appendix A. All these parameters involve spectral integrals, and have been 
computed numerically. 

The parameter fion is estimated in Appendix A. 

An altogether different mechanism for converting the baryons in nonlinear structures 
into ultraviolet photons is black hole accretion. If this mechanism is the dominant one, 
equation (3.2) should be replaced by 



' fbh = mass fraction of nonlinear structures that end up as black holes, 

face = fraction of rest energy radiated away during accretion process, 

fuy = fraction of energy released as UV photons, 

^ fesc = fraction of UV photons that escape from host galaxy. 

There is obviously a huge uncertainty in the factor /j,^. However, the absence of the factor 
0.0073 X 0.76 compared to equation (3.2) means that the conversion of matter into radiation 
is so much more efficient that the black hole contribution might be important even if fjjh 
is quite small. For instance, face = 10% and fesc = 100% gives fuvpp ~ ^^^fbhfuv, which 
could easily exceed the optimistic value fuvpp ~ 24, 000 for the stellar burning mechanism 
in Table 3.2. 

In Figure 3.2, the ionization fraction xi^) is plotted for various parameter values using 
equations (3.1) and (3.6). It is seen that the ionization grows quite abruptly, so that we 
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may speak of a fairly well-defined ionization redshift. Let us define the ionization redstiift 
Zion as the redshift when x becomes 0.5, i.e. 

1 + Zion = (1 + Zyir)eiic~^ i ""^ — ) . (3.7) 

This dependence of Zjo„ of the efficiency is shown in Figure 3.3 for our various galaxy 
formation scenarios. It is seen that the ionization redshift is fairly insensitive to the net 
efficiency, with the dependence being roughly logarithmic for fnet > 0.0001. 

3.4 Scattering History 

For a given ionization history x{^)i the Thomson opacity out to a redshift z, the probability 
that a CBR photon is Thomson scattered at least once after z, is 

P,(^) = l-e-^W, 

where the optical depth for Thomson scattering is given by 

1 - (l - ^) fHe] ^ ^ 0.057m 



Stt 



igmi 



where we have taken the mass fraction of helium to be fHe ~ 24% and assumed XHe ~ x, i.e. 
that helium never becomes doubly ionized and that the fraction that is singly ionized equals 
the fraction of hydrogen that is ionized. The latter is a very crude approximation, but makes 
a difference of only 6%. We assume that = 1 throughout this chapter, ^igm denotes the 
density of the intergalactic medium divided by the critical density, and is usually assumed 
to equal O;,, the corresponding density of baryons. the probability that a CBR photon is 
Thomson scattered at least once after the standard recombination epoch at ^ 10^. 
The profile of the last scattering surface is given by the so called visibility function 

/.(.) . f (.). 

which is the probability distribution for the redshift at which a photon last scattered. An 
illuminating special case is that of complete ionization at all times, i.e. x{z) = 1) which 
yields 



Ps{z) = 1 - exp (^-^r* [(1 + zf'^ - l]) «i 1 - exp - 



i) J 

for z S> 1 and hfligm = 0.03. Hence we see that in order for any significant fraction of the 
CBR to have been rescattered by reionization, the reionization must have occurred quite 
early. Figures 3.4 and 3.5 show the opacity and last-scattering profile for three different 
choices of h^igm- In the optimistic case Mligm = 0.1, it is seen that even as low an ionization 
redshift as Zjon = 30 would give a total opacity Pg ~ 50%. In Figures 3.6 and 3.7, we have 
replaced z by the angle subtended by the horizon radius at that redshift. 
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which is the largest angular scale on which Thomson scattering at z would affect the mi- 
crowave background radiation. In Figure 3.7, we have plotted the angular visibility function 
dPs/d{—6) instead of dPg/dz, so that the curves are probability distributions over angle 
instead of redshift. 

In the sudden approximation, the ionization history is a step function 

X{z) = 9{Zion - Z) 

for some constant zion, and as was discussed in Section 3.3, this models the actual ion- 
ization history fairly well. In this approximation, the visibility functions are identical to 
those in Figures 3.5 and 3.7 for z < Zion, but vanish between Zion and the recombination 
epoch at z « 10'^. Figure 3.8, which is in a sense the most important plot in this chapter, 
shows the total opacity Ps{ function of fnet for a variety of parameter values, as 

obtained by substituting equation (3.7) into (3.8). As can be seen, the resulting opacity 
is relatively insensitive to the poorly known parameter fnet, and depends mainly on the 
structure formation model {i.e. z^r) and the cosmological parameter hQ,igm- 

Thomson scattering between CBR photons and free electrons affects not only the spatial 
but also the spectral properties of the CBR. It has long been known that hot ionized IGM 
causes spectral distortions to the CBR, known as the Sunyaev-Zel'dovich effect. A useful 
measure of this distortion is the Comptonization y-parameter (Kompaneets 1957; Zel'dovich 
& Sunyaev 1969; Stebbins & Silk 1986; Bartlett & Stebbins 1991) 

yc= I (^^^rieatcdt, 

where the integral is to be taken from the reionization epoch to today. Let us estimate this 
integral by making the approximation that the IGM is cold and neutral until a redshift zicm, 
at which it suddenly becomes ionized, and after which it is remains completely ionized with 
a constant temperature T. Then for $7 = 1, Zjon S> 1, we obtain 

= (^) (^) r " ^ ^""*"-»-^^ 

where T4 = T/IO^K and Ueo, the electron density today, has been computed as before 
assuming a helium mass fraction of 24% that is singly ionized. Substituting the most recent 
observational constraint from the COBE FIRAS experiment, yc < 2.5 x 10~^ (Mather et al. 
1994), into this expression yields 

so all our scenarios are consistent with this spectral constraint. 

3.5 Discussion 

A detailed discussion of how reionization affects the microwave background anisotropics 
would be beyond the scope of this chapter, so we will merely review the main features. If 
the microwave background photons are rescattered at a redshift z, then the fluctuations 
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we observe today will be suppressed on angular scales smaller than the angle subtended by 
the horizon at that redshift. This effect is seen in numerical integrations of the linearized 
Boltzmann equation {e.g. Bond & Efstathiou 1984; Vittorio & Silk 1984), and can be simply 
understood in purely geometrical terms. Suppose we detect a microwave photon arriving 
from some direction in space. Where was it just after recombination? In the absence of 
reionization, it would have been precisely where it appears to be coming from, say 3000 Mpc 
away. If the IGM was reionized, however, the photon might have originated somewhere else, 
scattered off of a free electron and then started propagating towards us, so at recombination 
it might even have been right here. Thus to obtain the observed anisotropy, we have to 
convolve the anisotropics at last scattering with a window function that incorporates this 
smoothing effect. Typical widths for the window function appropriate to the last scattering 
surface range from a few arc-minutes with standard recombination to the value of a few 
degrees that we have derived here for early reionization models. 

In addition to this suppression on sub-horizon scales, new fluctuations will be generated 
by the first order Dopplcr effect and by the Vishniac effect. The latter dominates on small 
angular scales and is not included in the linearized Boltzmann treatment because it is a 
second order effect. The current upper limit on CBR fluctuations on the 1 arcminute 
scale of AT/T < 9 x 10~^ (Subrahmanyan et al. 1993) provides an interesting constraint 
on reionization histories through the Vishniac effect. In fact, according to the original 
calculations (Vishniac 1987), this would rule out most of the reionization histories in this 
chapter. However, a more careful treatment (Hu et al. 1994) predicts a Vishniac effect a 
factor of five smaller on this angular scale, so all reionization histories in this chapter are 
still permitted. 

The COBE DMR detection of AT/T has provided a normalization for predicting CBR 
anisotropics on degree scales. Several experiments are underway to measure such anisotropics, 
and early results that report possible detections have recently become available from exper- 
iments at the South Pole (Meinhold k, Lubin 1991; Shuster et al. 1993) and at balloon 
altitudes (Devlin et al. 1992; Meinhold et al. 1993; Shuster et al. 1993). There is some rea- 
son to believe that these detected signals are contaminated by galactic emission. Were this 
the case, the inferred CBR upper limits to fluctuations on degree scales might be lower than 
those predicted from COBE extrapolations that adopt the scale-invariant power spectrum 
that is consistent with the DMR result and is generally believed to be the most appropriate 
choice on large scales from theoretical considerations {e.g. Gorski et al. 1993; Kashlinsky 
1992). In the absence of such contaminations, the detected fluctuations in at least some 
degree-scale experiments are, however, consistent with the COBE extrapolation {e.g. Jubas 
&; Dodelson 1993). The variation from field to field, repeated on degree scales, also may 
argue either for galactic contamination or else for unknown experimental systematics, or 
even non-Gaussian fluctuations. The results of other recent experiments such as ARGO (de 
Bernardis et al. 1993), PYTHON (Dragovan et al. 1993) and MSAM (Cheng et al. 1993) 
have reinforced the impression that the experimental data is not entirely self-consistent, and 
that some form of systematic errors may be important. 
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The controversy over the interpretation of the degree-scale CBR fluctuations makes our 
reanalysis of the last scattering surface particularly timely. We have found that canonical 
dark matter, tailored to provide the 10 degree CBR fluctuations detected by the COBE 
DMR experiment, results in sufliciently early reionization (before z ^ 50) over a fairly wide 
range of parameter space, to smooth out primordial degree-scale fluctuations. Our middle- 
of-the-road model produces suppression by roughly a factor of two; it is difficult, although 
not impossible, to obtain a much larger suppression. This smoothing, because it is of order 
unity in scattering optical depth, is necessarily inhomogcncous. We predict the presence of 
regions with large fluctuations and many "hot spots" and "cold spots", corresponding to 
"holes" in the last-scattering surface, as well as regions with little small-scale power where 
the last scattering is more efficient. The detailed structure of the CBR sky in models with 
reionization will be left for future studies. Here we simply conclude by emphasizing that 
anomalously low values of AT/T over degree scales are a natural corollary of reionization 
at high redshift. 
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Figure 3.1: Virialization redshifts for objects of various masses. 

The virialization redshift, the redshift at which the bulk of the objects 
of mass Mc form, is plotted for a number of cosmological models. In 
all cases shown, + = 1 and 5c = 1-69. 
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Figure 3.2: Volume fraction ionized for various scenarios. 

The volume fraction of the universe that is in ionized Stromgren 
bubbles is plotted as a function of redshift for various parameter 
choices, corresponding to n = 1 CDM (optimistic, less optimistic, 
middle-of-the-road and pessimistic cases) and the tilted power spec- 
trum (n=0.7) variant of CDM. 
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Figure 3.3: Ionization redshift for various scenarios. 

The redshift at which x = 0.5 plotted as a function of the net effi- 
ciency. The four curves correspond to four of the choices of z^ir in 
Table 3.1: 41.4, 17.2, 8.4 and 4.8 from top to bottom. 
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Redshift 

Figure 3.4: Opacity for completely ionized IGM. 

The Thomson opacity Ps{z), the probability that a CBR photon has 
been scattered at least once after the redshift z, is plotted for four 
different choices of hQigm for the case where the IGM is completely 
ionized at all times. For more realistic scenarios where ionization 
occurs around some redshift Zion, the opacity curves simply level out 
and stay constant for z S> Zion- 
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Figure 3.5: Last-scattering surface for completely ionized IGM. 

The probability distribution for the redshift at which a CBR photon 
was last scattered, the so called visibility function, is plotted for four 
different choices of hi^igm for the case where the IGM is completely 
ionized at all times. For more realistic scenarios where ionization 
occurs around some redshift Zion, the curves are unaffected for z <C 
Zion, vanish for Zi(m <C ^; <C 10^ and have a second bump around 
z « 10=^. 
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Degrees subtended by horizon at reionization epoch 



Figure 3.6: Opacity for completely ionized IGM as function of angle. 

The total Thomson opacity Ps, the probability that a CBR photon 
has been scattered at least once since the recombination epoch, is 
plotted as a function of the angle in the sky that the horizon sub- 
tended at the reionization epoch. This is the largest angular scale on 
which fluctuations can be suppressed. 
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Figure 3.7: Last-scattering surface for completely ionized IGM as function of ang 

The probability distribution for the angle subtended by the horizon 
when a CBR photon was last scattered, the so angular visibility func- 
tion, is plotted for four different choices of hQig„i for the case where 
the IGM is completely ionized at all times. For more realistic scenar- 
ios where ionization occurs around some redshift Zion, corresponding 
to an angle 6ion, the curves are unaffected for 6 S> 9ion, vanish for 
2° ^ 6* ^ 9ion and have a second bump around 9 2°, the horizon 
angle at recombination. 
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Figure 3.8: Total opacity for various models. 



The total opacity, the probability that a CBR photon has been scat- 
tered at least once since the recombination epoch, is plotted for a 
variety of models. 
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Figure 3.9: Temperature evolution in intergalactic Stromgren bubbles. 

The temperature evolution is plotted for IGM exposed to a UV flux 
strong enough to keep it completely photoionized. In this example, 
h = 0.5, ^igm = 0.06, and T* = 36, 900. The upper dashed line is T* , 
the temperature corresponding to the average energy of the released 
photoelectrons, towards which the plasma is driven by recombina- 
tions followed by new photoionizations. The lower dashed line is the 
temperature of the CBR photons, towards which the plasma is driven 
Compton cooling. The three solid curves from left to right correspond 
to three different redshifts for becoming part of a Stromgren bubble. 
The first time the hydrogen becomes ionized, its temperature rises 
impulsively to T*/2. After this, Compton cooling rapidly pushes the 
temperature down to a quasi-equilibrium level, where the Compton 
cooling rate equals the recombination heating rate. 
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Figure 3.10: Ionization efficiencies for various scenarios 

The ionization efficiency, the fraction of the UV photons that pro- 
duce a net ionization, is plotted for three different parameter com- 
binations. In all cases, T* = 36, 900K, the value appropriate for the 
radiation from the population 3 star in Tabic 3.3. The solid lines are 
the exact results from numerical integration of equation (A. 5). The 
dotted lines are the analytic fits, which are seen to agree well in the 
redshift range of interest, which is typically z twice or three times 




Chapter 4 

Reionization in an Open Universe 



In this chapter, we generahze the results of the previous chapter to CDM models with J7 < 1. 
Such models have received recent interest because the excess power in the large-scale galaxy 
distribution is phenomenologically fit if the "shape parameter" T = hClo 0.25. It has been 
argued that essentially all CDM models may require early reionization to suppress degree- 
scale anisotropics in order to be consistent with experimental data, if the lowest degree- 
scale measurements are indeed characteristic of the primordial temperature anisotropics. It 
is found that if the cosmological constant A = 0, the extent of this suppression is quite 
insensitive to Qq, as opposing effects partially cancel. Given a ug-normalization today, the 
loss of small-scale power associated with a lower Qq is partially canceled by higher optical 
depth from longer lookback times and by structures forming at higher redshifts as the 
universe becomes curvature-dominated at z ~ ^- "^^^ maximum angular scale on which 
fluctuations arc suppressed decreases when Qq is lowered, but this effect is also rather weak 
and unlikely to be measurable in the near future. For flat models, on the other hand, where 
Aq = 1 — Oq, the negative effects of lowering dominate, and early reionization is not likely 
to play a significant role if <C 1. The same applies to for CDM models where the F is 
lowered by increasing the number of relativistic particle species. 

4.1 Introduction 

The inference from some, but not all, cosmic microwave background anisotropy experiments, 
is that many models, including the standard cold dark matter (CDM) model, may produce 
excessive temperature fluctuations on degree and sub-degree angular scales (Gorski et al. 
1993; Vittorio Sz Silk 1992). Early reionization at redshift z > 30 produces an optical 
depth to scattering r > 20%, and suffices, with the known baryon density from primordial 
nucleosynthesis constraints, to reconcile CDM with all observational limits (Sugiyama et al. 
1993). In Chapter 3, early photoionization of the intcrgalactic medium was discussed in a 
fairly model-independent way, in order to investigate whether early structures corresponding 
to rare Gaussian peaks in a CDM model indeed could photoionize the intcrgalactic medium 
sufficiently early to appreciably smooth out the microwave background fluctuations. In 
this chapter, the results of Chapter 3 will be generalized to < 1 models with non-zero 
cosmological constant A. Essentially all the notation used in this chapter was defined in 
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Chapter 3 and, in the interest of brevity, some of the definitions will not be repeated here. 

Just as in Chapter 3, our basic picture is the following: An ever larger fraction fg of 
the baryons in the universe falls into nonlinear structures and forms galaxies. A certain 
fraction of these baryons form stars or quasars which emit ultraviolet radiation. Some of 
this radiation escapes into the ambient intergalactic medium (IGM), which is consequently 
photoionized and heated. Due to cooling losses and recombinations, the net number of 
ionizations per UV photon, /jon, is generally less than unity. 

The results that we present here generalize the previous work to the case of an open 
universe. Lowering Qq has four distinct effects: 

1. Density fluctuations gradually stop growing once z < JIq ^. Thus given the observed 
power spectrum today, a lower implies that the first structures formed earlier. 

2. Matter-radiation equality occurs later, which shifts the turning-point of the CDM 

power spectrum toward larger scales. This means less power on very small scales 
(such as ~ 10^ Mq) relative to the scales at which we normalize the power spectrum 
(namely galaxy cluster scales, ~ 8/i~^Mpc or even the much larger COBE scale). One 
consequence is that the first structures form later. 

3. The lookback time to a given ionization redshift becomes larger, resulting in a higher 
optical depth. 

4. The horizon at a given ionization redshift subtends a smaller angle on the sky, thus 
lowering the angular scale below which CBR fluctuations are suppressed. 

Thus in terms of the virialization redshift Zyir defined in Chapter 3, the redshift at which 
typical structures go nonlinear, effect 1 increases z^r whereas effect 2 decreases Zmr- So 
these two effects influence fg, the fraction of baryons in nonlinear structures, in opposite 
directions. As to effect 2, it should be noted that this applies not only to CDM, but 
to any spectrum that "turns over" somewhere between the very smallest nonlinear scales 
(~ 10^ Mq) and the very largest (~ IO^^Mq) scales at which we COBE-normalize. Yet 
another effect of lowering is that the ionization efficiency fion drops slightly, at most by 

1/2 

a factor . This is completely negligible compared to the above-mentioned effects, as 
Zion depends only logarithmically on the efficiency. 

In the following sections, we will discuss each of these four effects in greater detail, and 
then compute their combined modification of the ionization history in a few scenarios. 

4.2 The Boost Factor 

When curvature and vacuum density are negligible, sub-horizon-sized density fluctuations 
simply grow as the scale factor aoc (l + z)^^. Thus at early times z 3> ^q^, we can write 

= —— do 

1 + z 
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for some function B independent of z that we will refer to as the boost factor. Clearly 
B{1,0) = 1. Thus if certain structures are assumed to form when S equals some fixed value, 
then given the observed power spectrum today, the boost factor tells us how much earlier 
these structures would form than they would in a standard flat universe. The boost factor 
is simply the inverse of the so called growth factor, and can be computed analytically for a 
number of special cases (see e.g. Peebles 1980). For the most general case, the fit 



2 



is accurate to within a few percent for all parameter values of cosmological interest (Carroll 
et al. 1992). The exact results are plotted in Figure 4.1 for the case Aq = and the flat 
case Ao = 1 — rio- Since we will limit ourselves to these two cases, the simple power-law fits 

rS(r!o,0) 

\B{n^,l-VL^) ^%'^-^\ 
which are accurate to within 1% for 0.2 < < 1, will suffice for our purposes. Note that 
the standard rule of thumb that perturbations stop growing at 1 + z k, VLq^, indicating 
B oc ^' particularly accurate in this context. 

4.3 The Power Spectrum Shift 

As mentioned above, lowering hn^ causes the first structures go nonlinear at a later redshift. 
This is quantified in the present section. 

The standard CDM model with power-law initial fluctuations proportional to /c" predicts 
a power spectrum that is well fitted by (Bond & Efstathiou 1984; Efstathiou et al. 1992) 

n 

P{k) oc ^ 



(l + [ag+(6g)i-5 + (cg)2]i-i3) 

where a = 6.4, b = 3.0, c = 1.7, q = (l/i~^Mpc)A;/r and the "shape parameter" F will be 
discussed further on. Although this fit breaks down for scales comparable to the curvature 
scale Tcurv = ^curv — ^o^\^ ~ ^o|~"^^^) it is quite accurate for the much smaller scales that 
will be considered in the present chapter. Rather, its main limitation is that it breaks down 
if 0,Q is so low that the baryon density becomes comparable to the density of cold dark 
matter. Thus for ^ligm ~ 0.05, the results cannot be taken too seriously for f^o < 0.2. We 
will limit ourselves to the standard n = 1 model here, as the tilted (n < 1) case was treated 
in Chapter 3 and was seen to be be essentially unable to reionize the universe early enough 
to be relevant to CBR anisotropics. The same applies to models with mixed hot and cold 
dark matter. 

Let us define the amplitude ratio 

i?(F,ri,r2)^^, 

where cr(ri) and a{r2) are the r.m.s. mass fluctuation amplitudes in spheres of radii ri and 
r-2, i.e. 

\2 /'°° n/i\ \sinkr cos/crl^ 
""^'^ io ^^^^ [ W " W 
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As in Chapter 3, we normalize the power spectrum so that a{8h~^Mpc) equals some constant 
denoted erg, and Mf. will denote the characteristic mass of the first galaxies to form. The 
corresponding comoving length scale Vc is given by Mc = ^irr^p, where p is the mean density 
of the universe. Thus given (jg, what is relevant for determining when the first galaxies form 
is the amplitude ratio 

R{r,rc,8h~^Mpc). 

This ratio is computed numerically, and the results are plotted as a function of F in Figure 4.2 
for a few different values of the cutoff mass Mc- It is easy to see why the amplitude ratio 
increases with F, since on a logarithmic scale, a decrease in F simply shifts the entire power 
spectrum towards lower k, thus decreasing the amount of power on very small scales relative 
to that on large scales. The fit 

R{T, Tc, 8/i-^Mpc) « 3 + 7.1 ln(l/i-^Mpc/rc)F 

is accurate to within 10% for 0.05 < F < 2 and lOOpc < Tc < lOOkpc. 



4.4 The Optical Depth 



Since a lower Qq implies a larger \dt/dz\ and an older universe, the optical depth out to a 
given ionization redshift Zion is greater for small ilo- For a given ionization history x(2;), the 
optical depth for Thomson scattering is given by 



=x{z')dz', 
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where we have taken the mass fraction of helium to be fne ~ 24% and assumed XHe ~ Xi 
i.e. that helium never becomes doubly ionized and that the fraction that is singly ionized 
equals the fraction of hydrogen that is ionized. The latter is a very crude approximation, 
but has the advantage that the error can never exceed 6%. If the universe is fully ionized 
for all redshifts below z, the integral can be done analytically for Aq = 0: 



t{z) 



2t* 
3ng 



1/2 



(4.1) 



for z S> ■ As is evident from the asymptotic behavior of the integrand, r is independent 
of Ao in the high redshift limit. Thus optical depth of unity is attained if reionization occurs 
at 



g2 I h^igm 
0.03 



-2/3 



a 



1/3 



4.5 The Angular Scale 



It is well known that reionization suppresses CBR fluctuations only on angular scales below 
the horizon scale at last scattering. Combining the standard expressions for horizon radius 
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{e.g. Kolb &; Turner 1990) and angular size, this angle is given by 



9 = 2 tan 



-1 



(4.2) 



noz - (2 - J^o)(Vl + ^^o-2 - 1) 



For z » fi. 



1-1 

'0 ' 



this reduces to 



e 




(4.3) 



but as is evident from Figure 4.3, this is quite a bad approximation except for z ^ 100. If 
we substitute it into equation (4.1) nonetheless, to get a rough estimate, we conclude that 
optical depth unity is obtained at an epoch whose horizon scale subtends the angle 



i.e., the dependence on all three of these cosmological parameters is relatively weak. As 
discussed in Chapter 3, fluctuations on angular scales much smaller than this are suppressed 
by a factor 



the opacity, which is the probability that a photon was Thompson scattered after redshift 
z. Its derivative, the visibility function fz = dP/dz, is the probability distribution for the 
redshift at which last scattering occurred, the profile of the last scattering surface. The 
angular visibility function 



is plotted in Figure 4.4 for the case where the universe never recombines (the curves for the 
more general case with reionization at some redshift ^jon can be read off from Figure 4.4 as 
described in the previous chapter). These functions give a good idea of the range of angular 
scales on which suppression starts to become important. In plotting these curves, the exact 
expression (4.2) has been used, rather than the approximation (4.3). It is seen that the 
qualitative behavior indicated by equation (4.4) is correct: as is lowered, the peak shifts 
down toward smaller angular scales, but the J7o-dependence is quite weak. 

4.6 Cosmological Consequences 

We will now compare the effect of lowering F in three cosmological models. The first model, 
which will be referred to as "open CDM" for short, has Aq = 0. The second model, referred 
to as "ACDM", has Aq = 1 — r^o- The shape parameter essentially tells us how early the 
epoch of matter-radiation equality occurred, and is given by 



where g^, = 3.36 corresponds to the standard model with no other relativistic degrees of 
freedom than photons and three massless neutrino species. In open CDM and ACDM, we 
have 7* = 3.36, so that VLq = T/h. The third model, referred to as rCDM (Dodelson, Gyuk 




(4.4) 



P{z) = 1 - e-^(^) 
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Pess. 


Mid. 


Opt. 


Very opt. 


0-8 


0.5 


1 


1.1 


1.2 




2.00 


1.69 


1.44 


1.33 


h 


0.5 


0.5 


0.8 


0.8 


Mc[Mq] 


10^ 


10^ 


10^ 


10^ 


f fionfuvpp 


1 


120 


23,000 


10^ 


erfc-i[l/2/] 


0.48 


2.03 


3.00 


3.55 




0.010 


0.013 


0.015 


0.020 



Table 4.1: Parameters used 



&; Turner 1994b), has Aq = and f2 = 1, and achieves a lower value of V by increasing 
instead. 

Including the effect of the boost factor, equation (3.7) in Chapter 3 becomes 

1 + = R{hnQ, rc, 8/i-^Mpc) B{VLo, Aq) erfc"^ 

The ionization redshift Zjon is plotted as a function of the shape parameter in Figure 4.5 
for the various scenarios specified in Table 4.1. It is seen that for the open model, the 
dependence otiVLq = T /h\s typically much weaker than the dependence on other parameters. 
One reason for this is that changes in the boost factor and the amplitude ratio partially 
cancel each other. For ACDM, the Jlo-dcpcndcnce is stronger, since the boost factor is 
weaker. In the rCDM model, the dependence on T /h is even stronger, as there is no boost 
factor whatsoever to offset the change in the amplitude ratio. 

The scenarios in Table 4.1 are similar to those in Chapter 3. In the one labeled "very 
optimistic", the high value for fuvpp-, the net number of produced UV photons per proton, 
is obtained by assuming that the main source of ionizing radiation is black hole accretion 
rather than conventional stars. Note that this speculative assumption still only increases 
Zion by 3.55/3.00 — 1 18%, the efficiency dependence being merely logarithmic. 

Figure 4.6, in a sense the most important plot in this chapter, shows the opacity as a 
function oiV/h for the various scenarios. Because of the increase in optical depth due to 
larger lookback times, the open model now gives slightly larger opacities for lower = F/Zi. 
However, this dependence is seen to be quite week. For ACDM, where the boost factor 
contributes less, the net result is seen to be the opposite; a slight decrease in the opacity 
for lower VLq = T/h. For the rCDM model, where there is neither a boost factor nor an 
increase in the lookback time, this drop in opacity is seen to be much sharper. Note that the 
dependence on other uncertain parameters, summarized by the four scenarios in Table 4.1, 
is quite strong. Indeed, this dependence is stronger than the effect of moderate changes in 
T/h, so in the near future, it appears unlikely that opacity limits will be able to constrain 
the shape parameter except perhaps in the rCDM model. 

The rCDM situation is summarized in Figure 4.7. To attain at least 50% opacity, hQigm 
must lie above the heavy curve corresponding to the scenario in question. On the other hand, 
nucleosynthesis (Malaney & Mathews 1993) places a strict upper bound on this quantity 



fuvpp f ton 
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if we assume that h > 0.5 It is seen that a shape parameter as low as F ~ 0.25, which 
would match large-scale structure observations (Peacock Sz Dodds 1994), is quite difficult 
to reconcile with these two constraints. 

4.7 Discussion 

Lowering F is an attractive resolution of the problem that arises in reconciling the observed 
structure in the universe on large scales with observations on megaparsec scales. The em- 
pirical power spectrum is well fit by F 0.25 (Peacock & Dodds 1994). Kamionkowski &; 
Spergel (1994) have found that primordial adiabatic fluctuations in an open universe with 
f2 0.3 are reconcilable with large-scale CBR anisotropy. On degree scales Kamionkowski, 
Spergel &; Sugiyama (1994) require reionization with optical depth r ~ 1 in order to reconcile 
the low density open model with recent experimental limits, if the lowest limits are adopted. 
In a low J7o ACDM model, the situation is not so critical, but reionization is required if the 
lowest limits (SP91) are adopted on degree scales (Gaier et al. 1992); r ~ 0.5 suffices 
however. A similar but slightly more favorable situation occurs in a rCDM model, where 
F = /iQo(fl'*/3.36)~^/^ is reduced by increasing by a factor of ~ 4, but some reionization 
is still required to match SP91. 

We have found that reionization giving r in the range 0.5 to 1 is readily produced and 
even natural in open models. This is because of the early formation of structure in combi- 
nation with the increased age of the universe, effects which compensate for the flattening 
of the power spectrum due to the delay in matter domination. However, the ACDM and 
rCDM models with low fare less well in this regard, since the loss of small-scale power 
is not balanced by significantly earlier structure formation. With Qf, in the range given by 
standard nucleosynthesis, a significant optical depth r > 0.5 is difficult to attain in either 
the ACDM or rCDM scenarios. 



^The rCDM model also alters the nucleosynthesis process (Dodelson, Gyuk & Turner 1994a), but this 
can only marginally relax the bounds unless the r neutrino is in a mass range incompatible with rCDM 
(Gyuk & Turner 1994). 
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Figure 4.1: The boost factor. 

The boost factor B{^q,Xq) is plotted as a function of for two 
classes of cosmologies. The upper curve corresponds to a standard 
open universe, i.e. Xq = 0, whereas the lower curve corresponds to 
fiat universes with Aq = 1 — i^o- 
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Figure 4.2: The amplitude ratio. 



The ratio of the fluctuation amplitude on the small scale rc to that 
at 8^~^Mpc is plotted as a function of the shape parameter F. From 
top to bottom, the four curves correspond to scales of 3.5/i~^kpc, 
7.5/i~^kpc, 16/i~^kpc and 35/i~^kpc, respectively. For h = 0.5 and 
fio = 1) these four length scales correspond to the masses IO^Mq, 
IO^Mq, 10''' Mq and lO^M©. The weak additional dependence on 
Qo/h that would result from holding Mc rather than rc fixed is clearly 
negligible, since as can be seen, Mc must vary by an entire order of 
magnitude to offset a mere 20% change in F. 
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Figure 4.3: The horizon angle. 

The angle in the sky subtended by a horizon volume at redshift z is 
plotted as a function of CIq for the case with no cosmological constant. 
The solid lines are the exact results for the four redshifts indicated, 
and the dashed lines are the corresponding fits using the simplistic 
approximation 6 ^ (Oq/^)^''^- 
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Degrees subtended by horizon at last scattering 



Figure 4.4: Visibility functions. 

The angular visibility function for a fully ionized universe is plotted 
for different values of fio and diffuse baryon content hQ,igm- The 
left group of curves corresponds to hQigm = 0.03 and the right to 
hO,igm = 0.1. Within each group, from left to right starting at the 
lowest peak, =0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9 and 1.0. 
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Figure 4.5: The ionization redshift. 



The ionization redshift is plotted for the four scenarios described in 
Table 4.1. The solid lines correspond to the open case where Aq = 0. 
The dashed lines correspond to the flat case where Aq = 1 — ^o- 
The dotted lines correspond to the rCDM model, where Aq = and 
Oq = 1- Note that the combination T/h is really an /i- independent 
quantity: for the open and flat cases, it is simply equal to Qq, and 
for the rCDM case it depends only on g^, the number of relativistic 
particle species. 
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Figure 4.6: The opacity. 

The opacity, the probabihty that a CBR photon is Thomson scat- 
tered at least once since since the standard recombination epoch, 
is plotted for the four scenarios described in Table 4.1. The solid 
lines correspond to the open case where Ao = 0. The dashed lines 
correspond to the flat case where Aq = 1 — Oq. The dotted lines cor- 
respond to the rCDM model, where Aq = and J7o = 1- Just as in 
Figure 4.5, note that the combination T/h is really an /i-independent 
quantity: for the open and flat cases, it is simply equal to 0,q, and 
for the rCDM case it depends only on , the number of relativistic 
particle species. 
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Figure 4.7: Reionization in rCDM. 

The two curves show the baryon density required for 50% opacity 
in rCDM, i.e. for reionization to rescatter 50% of the CBR pho- 
tons. The upper and lower heavy curves correspond to the middle- 
of-the road and optimistic scenarios, respectively. Thus even in the 
optimistic scenario, 50% opacity cannot be obtained outside of the 
fine-hatched region. The horizontal shaded region corresponds to the 
values of hO,b allowed by standard nucleosynthesis (0.01 < /i^Jlf, < 
0.015) in conjunction with the constraint 0.5 < h < 0.8. The verti- 
cal shaded region corresponds to values of the "shape parameter" T 
preferred by power spectrum measurements. Note that these three 
regions do not intersect. 



Chapter 5 

Did the Universe Recombine? 



In this chapter we will show that one still cannot conclusively assert that the universe 
underwent a neutral phase, despite the new COBE FIRAS limit y < 2.5 x 10"^ on Compton 
y-distortions of the cosmic microwave background. Although scenarios where the very early 
{z ~ 1000) ionization is thermal (caused by IGM temperatures exceeding lO^K) are clearly 
ruled out, there is a significant loophole for cosmologies with typical CDM parameters if the 
dominant ionization mechanism is photoionization. If the ionizing radiation has a typical 
quasar spectrum, then the y-constraint implies roughly h^^^Uig^QQ^'^^ < 0.06 for fully 
ionized models. This means that BDM models with J7o ~ 0.15 and reionization at z ~ 1000 
are strongly constrained even in this very conservative case, and can survive the y test only 
if most of the baryons form BDM around the reionization epoch. 

5.1 Introduction 

Recombination of the primeval plasma is commonly assumed but was by no means inevitable. 
Theories exist that predict early reionization are as diverse as those invoking primordial seed 
fluctuations that underwent early collapse and generated sources of ionizing radiation, and 
models involving decaying or annihilating particles. The former class includes cosmic strings 
and textures, as well as primordial isocurvature baryon fluctuations. The latter category 
includes baryon symmetric cosmologies as well as decaying exotic particles or neutrinos. 

The Compton y-distortion of the cosmic microwave background (CBR) provides a unique 
constraint on the epoch of reionization. In view of the extremely sensitive recent FIRAS 
limit of y < 2.5 X 10~^, we have reinvestigated constraints on the early ionization history of 
the intergalactic medium (IGM), and have chosen to focus on what we regard as the most 
important of the non-standard recombination history models, namely the primordial isocur- 
vature baryon scenario involving a universe dominated by baryonic dark matter (BDM), as 
advocated by Peebles (1987); Gnedin & Ostriker (1992) (hereafter "GO"); Gen, Ostrikcr & 
Peebles (1993) and others. This class of models takes the simplest matter content for the 
universe, namely baryons, to constitute dark matter in an amount that is directly observed 
and is even within the bounds of primordial nucleosynthesis, if interpreted liberally, and 
can reconstruct essentially all of the observed phenomena that constrain large-scale struc- 
ture. The BDM model is a non-starter unless the IGM underwent very early reionization. 
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in order to avoid producing excessive CBR fluctuations on degree scales. Fortunately, early 
nonlinearity is inevitable with BDM initial conditions, Sp/p oc M"^/^^, corresponding to a 
power-spectrum oc A;""^/^ for the observationally preferred choice of spectral index (Cen, 
Ostriker & Peebles 1993). 

Is it possible that the IGM has been highly ionized since close to the standard recombi- 
nation epoch at z ~ 1100? Perhaps the most carefully studied BDM scenario in which this 
happens is that by GO. In their scenario, 0,q = ri^o ~ 0.15. Shortly after recombination, a 
large fraction of the mass condenses into faint stars or massive black holes, releasing energy 
that reionizes the universe and heats it to T > 10, OOOK by z = 800, so Compton scattering 
off of hot electrons causes strong spectral distortions in the cosmic microwave background. 
The models in GO give a Compton y-parameter between 0.96 x 10"^ and 3.1 x 10~^, and 
are thus all ruled out by the most recent observational constraint from the COBE FIRAS 
experiment, y < 2.5 x 10""^ (Mather et al. 1994). 

There are essentially four mechanisms that can heat the IGM sufficiently to produce 
Compton y-distortions: 

• Photoionization heating from UV photons (Shapiro &; Giroux 1987; Donahue &; Shull 
1991) 

• Compton heating from UV photons 

• Mechanical heating from supernova-driven winds (Schwartz et al. 1975; Ikeuchi 1981; 
Ostriker & Cowie 1981) 

• Cosmic ray heating (Ginzburg &; Ozernoi 1965) 

The second effect tends to drive the IGM temperature towards two-thirds of the temperature 
of the ionizing radiation, whereas the first effect tends to drive the temperature towards a 
lower value T* that will be defined below. The third and fourth effect can produce much 
higher temperatures, often in the millions of degrees. The higher the temperature, the 

greater the y-distortion. 

In the GO models, the second effect dominates, which is why they fail so badly. In this 
chapter, we wish to place limits that are virtually impossible to evade. Thus we will use the 
most cautions assumptions possible, and assume that the latter three heating mechanisms 
are negligible. 

5.2 The Compton ^/-Parameter 

Thomson scattering between CBR photons and hot electrons affects the spectrum of the 
CBR. It has long been known that hot ionized IGM causes spectral distortions to the 
CBR, known as the Sunyaev-Zel'dovich effect. A useful measure of this distortion is the 
Comptonization y-parameter (Kompaneets 1957; Zel'dovich &: Sunyaev 1969; Stebbins &; 
Silk 1986; Bartlett k Stebbins 1991) 



5.3. IGM EVOLUTION IN THE STRONG UV FLUX LIMIT 



95 



where 

^) >■] (^) (^) 

Here Tg is the electron temperature, T^ is the CBR temperature, AT4 = (Te — r-y)/10^K, 
^igm is the fraction of critical density in intergalactic medium, and x{z) is the fraction of 
the hydrogen that is ionized at redshift z. Note that we may have ^igm ^ ^^6, i-e. all 
baryons may not be in diffuse form. The integral is to be taken from the reionization epoch 
to today. In estimating the electron density rig, we have taken the mass fraction of helium 
to be y 24% and assumed XHe ~ x, i.e. that helium never becomes doubly ionized and 
that the fraction that is singly ionized equals the fraction of hydrogen that is ionized. The 
latter is a very crude approximation, but makes a difference of only 6%. 

Let us estimate this integral by making the approximation that the IGM is cold and 
neutral until a redshift Zion, at which it suddenly becomes ionized, and after which it remains 
completely ionized with a constant temperature Tg. Then for Zion 3> 1 and Tg 3> Zion x 2.7K 
we obtain 

y « 6.4 X 10-^ hnigm^o '/^n 4n' 

where T4 = Te/IO^K. Substituting the most recent observational constraint from the COBE 
FIRAS experiment, y < 2.5 x 10~^ (Mather et al. 1994), into this expression yields 

..„.< 554r-^/^Qy^(^)"''. (5.2) 

Thus the only way to have Zion as high as 1100 is to have temperatures considerably below 
lO^K. In the following section, we will see to what extent this is possible. 



5.3 IGM Evolution in the Strong UV Flux Limit 

In this section, we will calculate the thermal evolution of IGM for which 

• the IGM remains almost completely ionized at all times, 

• the Compton y-distortion is minimized given this constraint. 

5.3.1 The ionization fraction 

In a homogeneous IGM at temperature T exposed to a density of C, UV photons of energy 
hv > 13.6 eV per proton, the ionization fraction x evolves as follows: 

= 7^=^ - ^) + A«x(l -x)- Kecx'] , (5.3) 

where Hq^{1 + z)~^ times the rates per baryon for photoionization, collisional ionization 
and recombination are given by 

' Xpi « 1.04 X 10^2 [hnigmais] C, 
Aci « 2.03 X lO^hntgmTl^^e-^^-^/^^' (5.4) 



^rec ~ O.TlT/lfijg^T^ 



1/2 



1.808 - 0.51nT4 + O.WT^ 



1/3 
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and T4 = Tg/IO^K. Here ais is the spectrally-averaged photoionization cross section in units 
of 10~^^cm^. The differential cross section is given by (Osterbrock 1974) 

dais fO ifi/<13.6eV, 

^(^)-{6.30$^i^ if^>13.6eV, (5-5) 

where 



hu ^ 



V 13.6 eV 

The recombination rate is the total to all hydrogenic levels (Scaton 1959; Spitzcr 1968). 
Recombinations directly to the ground state should be included here, since as will become 
evident below, the resulting UV photons are outnumbered by the UV photons that keep 
the IGM photoionized in the first place, and thus can be neglected when determining the 
equilibrium temperature. 

At high redshifts, the ionization and recombination rates greatly exceed the expansion 
rate of the universe, and the ionization level quickly adjusts to a quasi-static equilibrium 
value for which the expression in square brackets in equation (5.3) vanishes. In the absence of 
photoionization, an ionization fraction x close to unity requires Tg > 15, OOOK. Substituting 
this into equation (5.2) gives consistency with Zion > 1000 only if hQ.igm < 0.008, a value 
clearly inconsistent with the standard nucleosynthesis constraints (Smith et al. 1993). Thus 
any reheating scenario that relies on collisional ionization to keep the IGM ionized at all 
times may be considered ruled out by the COBE FIRAS data. 

However, this does not rule out all ionized universe scenarios, since photoionization can 
achieve the same ionization history while causing a much smaller y-distortion. The lowest 
temperatures (and hence the smallest y-distortions) compatible with high ionization will 
be obtained when the ionizing flux is so strong that \pi 3> Xd- In this limit, to a good 
approximation, equation (5.3) can be replaced by the following simple model for the IGM: 

• It is completely ionized {x = 1). 

• When a neutral hydrogen atom is formed through recombination, it is instantly pho- 
toionized again. 

Thus the only unknown parameter is the IGM temperature Tg, which determines the re- 
combination rate, which in turn equals the photoionization rate and thus determines the 
rate of heating. 

5.3.2 The spectral parameter T* 

The net effect of a recombination and subsequent photoionization is to remove the kinetic 
energy |fcT from the plasma and replace it with the kinetic energy ^kT*, where T* is defined 
by |/cT* = (Euv) - 13.6 eV and (Euv) is the average energy of the ionizing photons. Thus 
the higher the recombination rate, the faster this effect will tend to push the temperature 
towards T*. 
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UV source 


Spectrum P(i^) 






03 star 


T = 


50, OOOK Planck 


17.3 eV 


28,300K 


06 star 


T = 


40, OOOK Planck 


16.6 eV 


23,400K 


09 star 


T = 


30, OOOK Planck 


15.9 eV 


18,000K 


Pop. Ill star 


T = 


50, OOOK Vacca 


18.4 eV 


36,900K 


Black hole, QSO 


a = 


1 power law 


18.4 cV 


37,400K 


? 


a = 


2 power law 


17.2 eV 


27,800K 


? 


a = 


power law 


20.9 eV 


56,300K 


? 


T = 


100, OOOK Planck 


19.9 eV 


49,000K 



Table 5.1: Spectral parameters 



The average energy of the ionizing photons is given by the spectrum -P(i^) as {Euv) = 
h{i'), where 

^ J,^ P{u)a{u)du 

Here a is given by equation (5.5). Note that, in contrast to certain nebula calculations 
where all photons get absorbed sooner or later, the spectrum should be weighted by the 
photoionization cross section. This is because most photons never get absorbed, and all 
that is relevant is the energy distribution of those photons that do. Also note that -P(i^) 
is the energy distribution (W/Hz), not the number distribution which is proportional to 

The spectral parameters {Euv) and T* are given in Table 5.1 for some selected spectra 
(this is merely a subset of Table 3.3). A power law spectrum P{i^) oc z>^" with a = 1 fits 
observed QSO spectra rather well in the vicinity of the Lyman limit (Cheney &: Rowan- 
Robinson 1981; O'Brien et al. 1988), and is also consistent with the standard model for 
black hole accretion. A Planck spectrum P{i') oc v"' / ^g'^^/^'^ — 1^ gives a decent prediction 
of T* for stars with surface temperatures below 30, OOOK. For very hot stars, more realistic 
spectra (Vacca 1993) fall off much slower above the Lyman limit, thus giving higher values 
of T*. As seen in Table 5.1, an extremely metal poor star of surface temperature 50, OOOK 
gives roughly the same T* as QSO radiation. The only stars that are likely to be relevant to 
early photoionization scenarios are extremely hot and short-lived ones, since the universe is 
less than a million years old at z = 1000, and fainter stars would be unable to inject enough 
energy in so short a time. Conceivably, less massive stars could play a the dominant role 
later on, thus lowering T*. However, since they radiate such a small fraction of their energy 
above the Lyman limit, very large numbers would be needed, which could be difficult to 
reconcile with the absence of observations of Population HI stars today. 

5.3.3 The thermal evolution 

At the low temperatures involved, the two dominant cooling effects^ are Compton drag 

against the microwave background photons and cooling due to the adiabatic expansion of 

^Another cooling mechanism is coUisional excitation of atomic hydrogen followed by radiative de- 
excitation, which cools the IGM at a rate of (Dalgarno & McCray 1972) 

rtce ~ 7.5 X 10 e ' *n (1 — x)x ergcm s . 
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Model 


^0 




h 






y/0.000025 


Verdict 


QSO BDM I 


0.15 


0.15 


0.8 


37,400K 


1100 


6.30 


Ruled out 


QSO BDM II 


0.15 


0.15 


0.8 


37,400K 


200 


1.43 


Ruled out 


09 BDM 


0.15 


0.15 


0.8 


18,000K 


800 


2.91 


Ruled out 


QSO BDM III 


0.15 


0.04 


0.8 


37,400K 


1100 


0.67 


OK 


QSO CDM I 


1 


0.06 


0.5 


37,400K 


1100 


0.17 


OK 


QSO CDM II 


1 


0.03 


0.8 


37,400K 


1100 


0.16 


OK 



Table 5.2: Compton y-parameters for various scenarios 



the universe. Combining these effects, we obtain the following equation for the thermal 
evolution of the IGM: 



dT 2 _ 1 + z 



d{-z) 1 + z VI + ^oz 



KompiTy — T) + -Xrec{T)(T* — T) 



(5.6) 



where 



is (1 + z)^^ times the Compton cooling rate per Hubble time and = T^'q(1 + z). The 
factor of ^ in front of A^ec is due to the fact that the photoelectrons end up sharing their 
energy with the protons. We have taken T^^ « 2.726K (Mather et al. 1994). Numerical 
solutions to this equation are shown in Figure 5.1, and the resulting y-parameters are given 
in Table 5.2. 

The temperature evolution separates into two distinct phases. In the first phase, which 
is almost instantaneous due to the high recombination rates at low temperatures, T rises 
very rapidly, up to a quasi-equilibrium temperature slightly above the temperature of the 
microwave background photons. After this, in the second phase, T changes only slowly, 
and is approximately given by setting the expression in square brackets in equation (5.6) 
equal to zero. This quasi-equilibrium temperature is typically much lower than T*, since 
Compton cooling is so efficient at the high redshifts involved, and is given by 

AT = Te - « -^(T* - Te) cx -^g{T,)h^nigm{T* - T), (5.7) 

independent of J7o, where g{Te) 9S T~^ '^ encompasses the temperature dependence of A^ec- 
We typically have T <^ T* . Using this, making the crude approximation of neglecting 
the temperature dependence of Xreci ^^nd substituting equation (5.7) into equation (5.1) 
indicates that 

y^h'^lJl-^'^Tlz]!^. 



The ratio of this coohng rate to the Compton coohng rate is 

Hce _ exp[9.2-11.8/T4](l-z) 

He omp 

(1 + Z)T4 

a quantity which is much smaller than unity for any reasonable parameter values when T < lO^K. As will 
be seen, the temperatures at z ^ 1000 are typically a few thousand K, which with h^Qigm < 0.1 and 
X > 0.9 renders collisional excitation cooling more than nine orders of magnitude weater than Compton 
cooling. Hence we can safely neglect collisional excitations when computing the IGM temperature, the 
reason essentially being that the temperatures are so low that this process is suppressed by a huge Boltzman 
factor. 



5.4. CONCLUSIONS 



99 



Numerically selecting the best power-law fit, we find that this is indeed not too far off: the 
approximation 

y ^ (}.mi2h^-^n}gln-'/\T:f'{zion/noof' (5.8) 

is accurate to about 10% within the parameter range of cosmological interest. We have 
used equation (5.8) in Figure 5.3 by setting y = 2.5 x 10~^ and Zicm = 1100. The shaded 
region of parameter space is thus ruled out by the COBE FIRAS experiment for fully ionized 
scenarios. 

5.4 Conclusions 

A reanalysis of the Compton y-distortion arising from early reionization shows that despite 
the radical sharpening of the FIRAS limit on y, one still cannot conclusively assert that 
the universe underwent a neutral phase. Non-recombining scenarios where the ionization is 
thermal, caused by IGM temperatures exceeding lO^K, are clearly ruled out. Rather, the 
loophole is for the dominant ionization mechanism to be photoionization. We have shown 
that for spectra characteristic of both QSO radiation and massive metal-poor stars, the 
resulting IGM temperatures are so low that typical CDM models with no recombination 
can still survive the FIRAS test by a factor of six. This conclusion is valid if the flux of 
ionizing radiation is not so extreme that Compton heating becomes important. This is 
not difficult to arrange, as the cross section for Thomson scattering is some six orders of 
magnitude smaller than that for photoionization. 

For BDM models, the constraints are sharper. Non-recombining "classical" EDM models 
with ^ligm = r^o ~ 0.15 are ruled out even with the extremely cautious reheating assump- 
tions used in this chapter, the earliest ionization rcdshift allowed being z ~ 130. Such 
models involving early non-linear seeds that on energetic grounds can very plausibly pro- 
vide a photoionization source capable of reionizing the universe soon after the period of 
first recombination inevitably generate Compton distortions of order 10~^. These include 
texture as well as BDM models, both of which postulate, and indeed require, early reion- 
ization [z > 100) to avoid the generation of excessive anisotropy in the cosmic microwave 
background on degree angular scales. 

Thus BDM models with reionization at 2; 1000 can survive the y test only if most of 
the baryons form BDM when reionization occurs, and are thereby removed as a source of 
y-distortion, at least in the diffuse phase. This may be difficult to arrange at 2; > 100, since 
once the matter is reionized at this high a redshift, Compton drag is extremely effective in 
inhibiting any further gas collapse until z < 100. Since it takes only a small fraction of the 
baryons in the universe to provide a source of photons sufficient to maintain a fully ionized 
IGM even at z ~ 1000, we suspect that most of the baryons remain diffuse until Compton 
drag eventually becomes ineffective. Moreover, the possibility that the IGM is only partially 
reionized at z 1000 {e.g. GO), a situation which allows a lower value of the y-parameter, 
seems to us to be implausible as a delicate adjustment of ionization and recombination 
time-scales over a considerable range in z would be required. A complementary argument 
that greatly restricts the parameter space allowable for fully ionized BDM models appeals 
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to temperature fluctuations induced on the secondary last scattering surface, both by first 
order Doppler terms on degree scales and by second order terms on subarcminute scales (Hu 
et al. 1994). Thus, BDM models would seem to be in some difficulty because of the low 
limit on a possible y-distortion. 
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Figure 5.1: Thermal histories for various models. 

The temperature of the photoionized IGM is plotted for four of the 
cosmological models and spectra of ionizing radiation listed in Ta- 
ble 5.2. The lowermost curve gives the temperature of the CBR 
photons. 
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Figure 5.2: dy/dz iov various models 



The contribution to the y-parameter from different redshifts is plot- 
ted four of the cosmological models and spectra of ionizing radiation 
listed in Table 5.2. Thus for each model, the area under the curve 
is the predicted y-parameter. The area under the horizontal dashed 
line is 2.5 x 10-^ i.e. the COBE FIRAS limit. 
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Figure 5.3: Predicted and ruled out regions of parameter space. 

The hatched regions of parameter space are ruled out by the the 
COBE FIRAS limit y < 2.5 x 10"^ for Zion = HOG. J^o = 1 in the 
CDM plot and J7o = 0-15 in the BDM plot. The rectangular regions 
are the assumed parameter values for the CDM and BDM models, 
respectively. For CDM, the range 0.012 < h'^/^Qigm < 0.024 IS given 
by the nucleosynthesis constraint 0.010 < h^ri}, < 0.015 and the 
assumption that 0.5 < h < 0.8. (If 0,iym < ^b, the rectangle shifts to 
the left.) For the BDM models, h = 0.8 and 0.03 < ^igm < ^o- The 
vertical range corresponds to feasible values of the spectral parameter 
T*. The upper limit corresponds to highly speculative star with 
surface temperature 100, OOOK and T* = 49, OOOK. The lower line 
corresponds to an 09 star. The dotted horizontal line corresponds 
to the spectrum expected from quasars/accreting black holes. 
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Chapter 6 



Late Reionization by 
Supernova-Driven Winds 

In this chapter, a model is presented in which supernova-driven winds from early galaxies 
reionize the intergalactic medium by 2; = 5. This scenario can explain the observed absence 
of a Gunn-Pctcrson trough in the spectra of high-redshift quasars providing that the bulk 
of these early galaxies are quite small, no more massive than about 10® Mq. It also predicts 
that most of the IGM was enriched to at least 10% of current metal content by z = 5 and 
perhaps as early as 2: = 15. The existence of such early mini-galaxies violates no spectral 
constraints and is consistent with a pure CDM model with b < 2. Since the final radius 
of a typical ionized bubble is only around 100 kpc, the induced modification of the galaxy 
autocorrelation function is negligible, as is the induced angular smoothing of the CBR. Some 
of the gas swept up by shells may be observable as pressure-supported Lyman-alpha forest 
clouds. 

6.1 Introduction 

The absence of a Gunn-Peterson trough in the spectra of high-redshift quasars has provided 
strong evidence for the intergalactic medium (IGM) being highly ionized as early as 2; = 4 
(Gunn & Peterson 1965; Steidel and Sargent 1987; Webb et al. 1992). The hypothesis that 
photoionization of the IGM by quasars could account for this ionization (Arons &; McCray 
1969; Bergeron Sz Salpeter 1970; Sherman 1980) has been challenged (Shapiro 1986; Shapiro 
Sz Giroux 1987; Miralda-Escude & Ostriker 1990). Other studies have maintained that pho- 
toionization by quasars (Donahue &; ShuU 1987) or active galactic nuclei (Teresawa 1992) 
may nonetheless be sufficient. However, in view of the large uncertainties in crucial param- 
eters such as ionizing fluxes, the issue of what reionized the IGM must still be considered 
open. 

In comparing the Gunn-Peterson constraints with our work in Chapters 3 and 4, the 
crucial difference is the degree of ionization required. To affect the CBR, it does not really 
matter whether the ionization fraction x is 90% or 99.999%, as this makes only a 10% 
difference in the optical depth r. The Gunn-Pctcrson limits constrain not x but (1 — x), the 
neutral fraction. Thus in this context, the difference between 90% and 99.999% ionization 
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is four orders of magnitude. In Chapters 3 and 4, we found that photoionization by early 
galaxies could easily reionize the IGM by a redshift z = 5, but the issue here is whether 
photoionization alone can provide the extremely high ionization fraction required to pass 
the Gunn- Peterson test. 

In this chapter, we investigate an alternative rcionization scenario, which produces con- 
siderably higher IGM temperatures than those attained by the photoionization models in 
previous chapters. Supernova driven winds from luminous galaxies have long been conjec- 
tured to be an important ionization source for the IGM (Schwartz et al. 1975; Ikeuchi & 
Ostriker 1986; Carlberg & Couchman 1989). Cold dark matter (CDM)-based models of 
structure formation (Blumenthal et al. 1984; Efstathiou et al. 1985) predict the forma- 
tion of gravitationally bound objects of mass as small as 10^ Mq in large numbers before 
z = h. Recent work (Blanchard et al. 1992) indicates that such objects can cool rapidly 
and presumably fragment into stars. These early mini-galaxies would be expected to release 
great amounts of kinetic energy into the surrounding IGM, thereby creating large, fairly 
spherical voids filled with thin, hot, ionized plasma. We analyze the effect of expanding 
bubbles driven by supernova winds from early mini-galaxies, and show that this mechanism 
of distributing energy can indeed provide the required ionization without violating any of 
the current spectral constraints. 

In Section 6.2, we will treat the expansion of a shell in a uniform, cold and neutral IGM. 
As these bubbles become larger and more numerous and fill most of space, this obviously 
becomes a very poor model of the IGM. In Section 6.3 we estimate bulk properties of this 
new processed IGM such as temperature, density and ionization. 

6.2 The Explosion Model 

Since the pioneering work on spherically symmetric explosions by Sedov (1959), a profusion 
of analytic solutions have been given by numerous authors for models of ever-increasing 
complexity (Cox & Smith 1974; McKee & Ostriker 1977; Weaver et al. 1977; McCray & 
Snow 1979; Bruhweiler et al. 1980; Tomisaka et al. 1980; McCray & Kafatos 1987; Ostriker 
& McKee 1988). Most of these models pertain to bubbles in the interstellar medium of a 
galaxy, where the expansion of the universe can be ignored. Ostriker & McKee have given 
asymptotic self-similarity solutions that incorporate this latter complication, but unfortu- 
nately they are not sufficiently accurate for our needs. The reason is that since neither 
energy nor momentum is conserved in the regime before the shell becomes self-similar, there 
is no accurate way to normalize the self-similar solution using the initial data. 

Let and denote the average densities of baryonic and non-baryonic matter in the 
universe. We will assume that all baryons are in diffuse for early on, so that is also the 
density of the IGM. We will write p;, = ri^pc and p^ = ^dPc-, where the critical density 
Pc = SH^/SttG. We will use a three-phase model for the expanding bubbles: 

• a dense, fairly cool spherical shell of outer radius R and thickness R5, containing a 
fraction (1 — /^) of the total baryonic mass enclosed. 
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• uniform neutral ambient intergalactic medium (IGM) of density + pd and zero 
pressure outside, 

• a hot, thin, isothermal plasma of pressure p and temperature T inside the shell. 

The shell is driven outwards by the pressure of the hot interior plasma but slowed by the 
IGM and by gravity. The plasma is heated by kinetic energy from supernova explosions and 
collisions with IGM and cooled by bremsstrahlung and Compton drag against the cosmic 
background radiation. 

6.2.1 The expanding shell 

We assume that the expanding shell sweeps up almost all baryonic IGM that it encounters 
and loses only a small fraction of it through evaporation into the interior, so that its total 
mass is given by m(t) = ^TrR(t)^{l — fjn)pb, where the constant ^ 1. Since Pb/Pb = — 3-H" 
for any cosmological model, we get 

— = (R^pb) — (R^Pb) = 3\ —- h] if — >i7, zero otherwise. 
m ^ at ^ ) \R yit 

(The shell will acquire new mass when it is expanding faster than the Hubble flow, and will 
never lose mass.) It turns out that the Hubble flow catches up with the shell only as t — ^ oo, 
so we will always have R > HR and m > 0. 

When new mass is swept up, it must be accelerated from the velocity HR to R, so the 
shell experiences a net braking force {R — HR)m. The interior pressure p drives the shell 
outward with a force pA = AttR'^p = 3mp/pbR in the thin shell approximation 6, fm ^ 1- 
Finally there is a gravitational braking force, which in the thin-shell approximation (Ostriker 
Sz McKee 1988) gives the deceleration ^7rGR{pd + ^Pb)- Adding these three force terms, the 
radial equation of motion becomes 

6.2.2 The interior plasma 

The equation of state for the plasma in the hot interior gives the thermal energy 

Et = ^pV = 2'KpR\ (6.2) 

and energy conservation for the interior yields 

Et = L- pdV/dt = L- iirpR'^R, (6.3) 

where the luminosity L incorporates all sources of heating and cooling of the plasma. We 
will consider five contributions to L and write 

~ -^sn -l-'comp -l-'brems -^ion ~l" -^diss! 
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where Lgn is the energy injection from supernova explosions, Lcomp the cooUng by Compton 
drag against the CBR, Li^e^^s the coohng by bremsstrahlung, Ljo„ the coohng by ionization 
of neutral IGM and L^iss the heating from collisions between the shell and the IGM. 

In stellar burning from zero to solar metallicity, the mass fraction 0.02 x 0.007 is released, 
mostly as radiation. Dtic to low cross-sections, only a negligible fraction of this radiation 
will contribute towards heating the gas, so we will only be interested in the energy that 
is released in kinetic form. From empirical observations of active galactic winds (Hcckman 
1990) about 2% of the total luminosity from a galaxy is mechanical. Another empirical 
observation is that for a solar neighborhood initial stellar mass function, one has roughly 
one supernova for every 150Mq of baryons that form stars, with a typical kinetic energy 
output of 10^1 ergs per explosion. Both of these observations lead to the same estimate 

Lsn = — 7 ^{tburn — i) ~ ^■'^LQ——9{thurn ~ 0> 

t^burn -'"O 

where the efficiency fsn ~ 4 x 10~^ and where we have assumed that the energy is released 
at a constant rate during a period tburn ~ 5 x lO'^ years. 

Now let us examine cooling. The interior baryon density is pi = Pbfm/{^ — <^)'^ whereas 
the shell density is ps = pb{^ — fm)/{^ ~ (1 ~ ^T) ~ if fm-^ ^ 1- Compton drag 

against the microwave background radiation causes energy loss at a rate (Kompaneets 1957) 



Lcomp = ^ (^.cne) (^) ncV, (6.4) 



where ut is the Thomson cross section, V = |7ri?^, and T^ = T, the temperature of the 
interior plasma, which is given by 

E, = [l + l) kT^V. 



2 2/ nip 

(We will assume almost complete ionization and low metallicity, so that He ~ fmnb-) Using 
equation (6.2), we see that Compton drag causes cooling on a timescale 

Et 45 ( 'hc\ rUe „ ,„i9 ,^ a 

—J— = —Lr^2x 10^2 years x (1 + z)"^, 

that is, it becomes important only at high redshifts. It turns out that Lbrems ^ Lcomp in our 
regime of interest, so we will simply make the approximation Lbrems ~ 0- Assuming that 
the ambient IGM is completely neutral, the power required to ionize the hydrogen entering 
the interior is simply 



Lion = fmnbEo X AttR^ 



R-HR 



where Eo 13.6 eF. 

The equation of motion (6.1) assumes that the collisions between the expanding shell 
and the ambient IGM are perfectly inelastic. The kinetic energy dissipated has one of three 
fates: It may 



(a) radiate away in shock cooling processes. 
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(b) ionize the swept up IGM, or 

(c) heat the shell and by conduction heat the interior plasma. 

Let fd denote the fraction that is reinjected into the interior plasma through processes 
(b) and (c). This is one of the major uncertainties of the model. Now a straightforward 
kinematic calculation of the kinetic energy loss per unit time gives 

Making accurate estimates of fa is difficult, so we simply use the two extreme cases fd = 
and = 1 in the simulations. Perhaps surprisingly, the results will be seen to be relatively 
independent of the choice of fd. 

6.2.3 Solutions to the equations 

Combining (6.2) and (6.3) leaves us with 

L R 

The system of equations (6.1) and (6.5) reduces to that derived by Weaver et al. (1977) in 
the special case where L{t) is constant and the expansion of the universe is ignored. 
Let us define dimensionless variables as follows: 

I = L/L,, L, = fsnMf,c\^^^ ^ 1.2 X 10^ Lq X M; 
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e= E/E^, E^= L^hurn « 7.2 X lO^'^ erg X M; 

r = R/R^, R^ = lI'^G^/^U ~ 0.13 Mpc x h-^{l + z^)~^''^mI'^ 

q = p/p^, p^ = lI'^G-^/H-^ « 1.4 X 10-16 Pa x h'^{l + z^ful'^ 

Here we have taken /i = 0.5 and defined M5 = M^/IO^ Mq. If f2 = fi;, + ^2^=1, then is 
the age of the universe at the rcdshift z^, when the shell begins its expansion, i.e. the Big 
Bang occurred at r = — 1 and the shell starts expanding at r = 0. For this simple case, we 
have r/ = (1 + t)-^ = (1 + z*)-3/2(l + zfl^. 
Equations (6.1) and (6.5) now become 

r"{r) = iff .(r)-$l - 3 (l - f^)' ^ - {p, + 1^,) ,(r)V(r) 

H6r6 ^ ^sn ^comp ^ion ~l~ ^diss-) wh.6r6 

icomp ~ 0.0m-l(l + Z*)-3/2(l + z)Vg, 

iion ~ 2.2fAM~'^/^ X rfr"^ [r' - |r?r) , and 
hiss = ^d^b X {wf ir' - Iw) ■ 
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In computing icomp, we have taken T-yo = 2.74K. The interior temperature, the thermal 
energy and the kinetic energy are given by 

T ^ 4.5 X 10' K X 



£t = r q, 

Tburn 

mE?/2 2 3 /2 

£k = = Q ^ ■ 

The solution to the system (6.6) evolves through three qualitatively different regimes: 
T ^ 1, r 1 and r S> 1. 

a) In the limit of small times r <C 1 , gravity and Hubble flow are negligible and we obtain 
the asymptotic power law solution 

r-(r) = ar^/^ ?(r) = hT'^'^ 

as may be verified by direct substitution. This solution reduces to that found by 
Weaver et al. in the special case fd = 0- Since the total energy injected is simply 
£m = r/Tburn, t^is gives 

£t 35 £k 15 

and — 



Sin 77 - 27 fd e^n 77 - 27 fd 

for small r. Hence even though et+£k = Em only for the most optimistic case /d = 1, 

5+1. 
77 



we see that no more than 1 - « 25% of the injected energy is lost as radiation 



even in the worst case /d = 0. 

b) The behavior in the intermediate regime is a complicated interplay between several 
different effects: 

1. After approximately 5 x 10^ years, the supernova explosions cease, which slows 
the expansion. In this pressure-driven snowplow phase, we would asymptotically 
have R oc t^/^, t^^^^, t^/'^, t^^^^ or t^^' if there were no gravity, no Hubble flow 
and no cooling with fd = 0, |, |, | or 1, respectively. 

2. Cooling (and pdV) work reduces the pressure and the thermal energy to virtually 
zero, which slows the expansion. With zero pressure, we would approach the 
momentum-conserving snowplow solution R oc t^^'^ if there were no gravity and 
no Hubble flow. 

3. The density of the IGM drops and the IGM already has an outward Hubble 
velocity before it gets swept up, which boosts the expansion and adds kinetic 
energy to the shell. 

4. Gravity slows the expansion. 
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5. Dark matter that has been accelerated outward by the shell catches up with it 
again and speeds up the expansion. (This last effect has been neglected in the 
equations above, since it generally happens too late to be of importance for our 
purposes.) 

c) As t —>■ oo, the shell gets frozen into the Hubble flow, i.e. R oc t^/^ if $7 = 1. An 
approximate analytic solution for r S> 1 is given by Ostriker&; McKee (1988), but 
since neither energy nor momentum is conserved in the intermediate regime, there is 
no simple way to connect this solution with the short-time solution above. 

Numerical solutions for the comoving radius (1 + z)R are plotted in Figure 6.1 for different 
values of and fd. The asymptotic solution (6.7) has been used to generate initial data 
at r = 0.01 for the numerical integration. In this Figure, we have truncated R when the 
interior temperature drops below 15,000K, after which newly swept up IGM fails to become 
ionized. Figures 6.2, 6.3 and 6.4 show what becomes of the injected energy for different 
parameter values. Note that the relative fractions are approximately constant early on, 
while the supernovae inject energy, in accordance with the asymptotic solution (6.7). The 
reason that the total energy exceeds 100% of the input is that the shell gobbles up kinetic 
energy from swept-up IGM that already has an outward Hubble velocity. 

6.3 Cosmological Consequences 

Once the expanding bubbles discussed in the previous section have penetrated most of space, 
the IGM will presumably have a frothy character on scales of a few 100 kpc, containing thick 
and fairly cool shell fragments separated by large, hot, thin and ionized regions that used 
to be bubble interiors. 

In Section 6.3.1, we calculate at what point the IGM becomes frothy, more specifically 
what fraction of space is covered by expanding shells at each z. In 6.3.2 we discuss the 
resulting enrichment of the IGM with heavy elements. In 6.3.3 the thermal history of the 
IGM after this epoch is treated. Finally, in 6.3.4 the residual ionization is computed, given 
this thermal history, and we discuss the circumstances under which the Gunn-Peterson 
constraint is satisfied. 

6.3.1 IGM porosity 

Assuming the standard PS theory of structure formation (Press & Schechter 1974), the 
fraction of all mass that has formed gravitationally bound objects of total (baryonic and 
non-baryonic) mass greater than M at redshift z is 

1 A ' 

1 — erf —= 

V2a{M) 

where erf (z) = 2tt~^^'^ Jq e~^^ du and is the linearly extrapolated r.m.s. mass fluctuation 
in a sphere of radius tq. The latter is given by top-hat filtering of the power spectrum as 



o- = — — oc 



2 1 ,.00 r g^j^ 1 2 



1 + zJ {l + zyJo 



P{k)dk, (6.8) 
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where ro is given by ^ttvqP = M and where P{k) is the power spectrum. Although this 
approach has been criticized as too simphstic, numerical simulations (Efstathiou et al. 1988; 
Efstathiou Sz Rees 1988; Carlberg & Couchman 1989) have shown that it describes the mass 
distribution of newly formed structures remarkably well. Making the standard assumption 

of a Gaussian density field, Blanchard et al. (1992) have argued that it is an accurate 
description at least in the low mass limit. Since we are interested only in extremely low 
masses such as W^Mq, it appears to suffice for our purposes. 

Wc choose 6c = 1.69, which is the linearly extrapolated overdensity at which a spherically 
symmetric perturbation has collapsed into a virializcd object (Gott &: Rees 1975). Letting 
fg denote the fraction of all baryons in galaxies of mass greater than M at z, this would 
imply that 

^1.69(1 + z)' 



/p 1 - erf 



(6.9) 



if no other forces than gravity were at work. However, it is commonly believed that galaxies 
correspond only to such objects that are able to cool (and fragment into stars) in a dynamical 
time or a Hubble time (Binney 1977; Rees & Ostriker 1977; Silk 1977; White & Rees 1978). 
Hence the above value of fg should be interpreted only as an upper limit. 

A common assumption is that the first galaxies to form have a total (baryonic and dark) 
mass Mc ~ 10^ M©, roughly the Jeans mass at recombination. Blanchard et al. (1992) 
examine the interplay between cooling and gravitational collapse in considerable detail, 
and conclude that the first galaxies to form have masses in the range IO^Mq to 10^ Mo, 
their redshift distribution still being given by equation (6.9). To keep things simple we 
will assume that all early galaxies have the same mass Mc and compare the results for 
Mc = 2 X IO^Mq, IO^Mq and 10^^ Mq. 

Let R{z;z^:) denote the radius of a shell at z that was created at by a galaxy of 
baryonic mass Mf, = ^IbMc as in Section 6.2. Then the naive filling factor, the total bubble 
volume per unit volume of the universe, is 



oo 



where 



3'''^^^'^*^ M6d(-z.) ^ Jz (i+z^r/^ di-z,)"^'*' ^^-'^^ 



(f>^ fii 1600/t-^M5 ^/^(n;,/0.06). 



Clearly nothing prohibits (p from exceeding unity. This means that nearby shells have 
encountered each other and that certain volumes are being counted more than once. If the 
locations of the bubbles are uncorrelated, then the fraction of the universe that will be in a 
bubble, the porosity, is given by 

P = 1 - e-"^. 

If the early galaxies arc clustered rather than Poisson-distributed, this value is an over- 
estimate. For an extreme (and very unrealistic) example, if they would always come in 
clusters of size n and the clusters would be much smaller than the typical bubble size of 
100 kpc, then it is easy to see that P m 1 — e~^/^. For more realistic cases, simple analytic 
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(^o(Me) 


bcdmfi 


bcdm,8 


hiltedfi 


btilted,8 


1% 


3.94 


4.8 


3.5 


2.5 


2.0 


10% 


6.18 


3.1 


2.2 


1.6 


1.3 


20% 


7.92 


2.4 


1.7 


1.2 


1.0 


50% 


15.06 


1.3 


0.9 


0.6 


0.5 



Table 6.1: Correspondence between various ways of normalizing the power spectrum 

expressions for P are generally out of reach. Since we expect the clustering to be quite 
weak, we will use the Poisson assumption for simplicity. 

The uppermost panels of Figures 6.5 and 6.6 contain P{z) for various parameter values, 
calculated numerically from equation (6.10) using the numerical solutions for r{z-,z^..). It is 
seen that the lower mass in Figure 6.5 (2 x 10^ M0 versus 10^ M©) gives higher filUng factors, 
so that the expanding shells fill almost 100% of space by z = 5 for three of the four choices 
of /p(5). In 6.6, we see that almost 20% of the baryons must be in galaxies by 2; = 5 to 
achieve this. The greater efficiency of small galaxies is to be expected, since 4>^ OC M-2/5. 
Although some parameters still yield the desired P 100% by z = 5 in Figure 6.6, using 
present-day masses like Mc = 10^^ M© fails dismally (not plotted) for all choices of the other 
parameters. Roughly, the largest Mc that works is 10^ M© 

As can be seen, the dependence on (dashed versus solid lines) is rather weak. 

In order to calculate ao{Mc) from the fluctuations observed on larger scales today, we 
need detailed knowledge of the power spectrum down to very small scales, something which 
is fraught with considerable uncertainty. For this reason, we have chosen to label the curves 
by the more physical parameter /g(5), the fraction of all baryons that have formed galaxies 
by z = 5. The four sets of curves correspond to fractions of 50%, 20%, 10% and 1%. 
These percentages should be compared with observational estimates of metallicity, as will 
be discussed in Section 6.3.2. 

The second column of Table 6.1 contains the values (Jo(Mc) necessary to obtain various 
values of /g(5), calculated by inverting the error function in equation (6.9). The last four 
columns contain the bias factors necessary to yield this value of ao{Mc) for two choices 
of power spectra (CDM and n=0.7 tilted CDM) and two choices of cutoff mass (Mc = 
2 X 10^ M0 and Mc = 10^ Mq). Thus b = j/ao{Mc), where we define 7 to be the ratio 
between a at Mc and a at 8h~^ Mpc = b~^. Performing the integral (6.8) numerically with 
the CDM transfer function given by Bardeen et al. 1986 (BBKS), h = 0.5, ^ = 1, Qb <^ 1 
and an n = 1 Harrison-Zel'dovich initial spectrum gives 7 ^ 19.0 for Mc = 2 x 10^ M© 
and 7 13.6 for Mc = 10^ M©. Using the CDM transfer function of Bond and Efstathiou 
(1984) instead gives 7 Ri 18.1 and 7 ^ 13.7, respectively. The BBKS transfer function is 
more applicable here since it includes the logarithmic dependence that becomes important 
for very low masses. The BBKS transfer function with a tilted (n=0.7) primordial spectrum 
yields the significantly lower values 7 « 9.71 and 7 ~ 7.93, respectively. 
Basically, Table 6.1 shows that any of our values of fg{5) become consistent with a feasible 
bias factor for some choice of power spectrum and cutoff mass. 
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6.3.2 IGM enrichment 

These values of fg{5) should be compared with observational estimates of metallicity, since 
if the stars in these early mini-galaxies produce the same fractions of heavy elements as 
do conventional stars, then these percentages are directly linked to the fraction of currently 
observed metals that were made before z = 5. Some of the enriched shells may be observable 
as quasar absorption line systems, as intracluster gas, and, indirectly, as in the metallicities 
of old disk and halo stars. 

Observations of iron abundances in intracluster gas by HEAO-1, Exosat and Ginga {e.g. 
Mushotzky 1984; Hughes et al. 1988; Edge 1989; Hatsukade 1989) have shown that most 
clusters have abundances between 25% and 50% of the solar value. Einstein observations 
have showed the presence of a large variety of other heavy elements in the intracluster gas 
(Lea et al. 1982; Rothenfiug et al. 1984). Most of this gas and some of these metals are 
believed to be "primordial" , since the gas mass in clusters is typically several times greater 
than the observed stellar mass in the cluster galaxies (Blumenthal et al. 1984; David et al. 
1990; Arnaud et al. 1991). 

There are indications that the most of these heavy elements may have been produced as 
recently as around z = 2 — 3, and that the metallicity in the halo gas of some 2; ~ 3 galaxies 
inferred from QSO absorption line studies are as low as 0.1% of the solar value (Steidel 
1990). However, this and other observations of extremely metal-poor objects (Pettini et al. 
1990) does not necessarily rule out our scenario, since it is highly uncertain whether all the 
hydrogen in the swept-up IGM would get thoroughly mixed with the metal-rich supernova 
ejecta. 

6.3.3 IGM temperature 

Let T{z; z*) denote the temperature of the interior of a bubble at z that was created at 
as in Section 6.2. Then the volume-averaged temperature of the IGM is 



When (j) becomes of order unity, the IGM swept up by the expanding shells is no longer 
cold, neutral and homogeneous, so the treatment in Section 6.2 breaks down. The resulting 
temperatures will be underestimated, since less thermal energy needs to be expended on 
heating and ionization. As can be seen in Figures 6.5 and 6.6, the transition from cf) <^ 1, 
where the treatment in Section 6.2 is valid, to ^ 3> 1, where the IGM becomes fairly uniform, 
is quite rapid. Since Tjgm defined above is proportional to the thermal energy per unit 
volume, energy conservation leads us to assume that Tigm remains fairly constant during 
this transition and therefore is a good estimate of the bulk IGM temperature immediately 
afterwards. From this time on, we will approximate the IGM outside the scattered dense 
and cold shell remnants by a uniform isothermal plasma. Applying equations (6.2), (6.3) 




(6.11) 
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and (6.4) to the IGM yields the following equation for its thermal evolution: 

n+iinj, (6.12) 



—T 
dz ^ 



2 +^(1 + ^)3/2' 



1 + z 

where the first term encompasses cooling from adiabatic expansion and the second term 
Compton cooling. T5 = T/gm/IO^ K, A = l.5{to/tcomp) ~ 0.0042/i~^ and £inj = t^Linj/k x 
lO^K, where is the power injected into the IGM per proton from all heat sources com- 
bined. The Compton cooling term is seen to increase with redshift, equaling the adiabatic 
term at z « VlhP-^ - 1. 

In the most pessimistic case of no reheating whatsoever, i.e. for l^nj = Oj equation (6.12) 
has the solution 

roc(i + z)V-4^(^+^)'''. 

A more optimistic assumption is that some fraction finj of the total energy released from 
stellar burning in newly formed galaxies continues to heat the IGM, i.e. 

/0.O2 X 0.007mpc2\ dfg /2.1xl0^\ / 1 [1.69(1 + z)1 2\ 

- J^nj y ^ X lO^if j [ ao ) ''''^ [~2 [ ao(Me) J ) ' 

which would also incorporate other modes of energy injection such as radiation. The middle 
panels in Figures 6.5 and 6.6 show the temperature resulting from /j„j = 0,0.001 and 
0.02 and different initial values. For P < 0.8, Tjgm has been calculated numerically from 
equation (6.11) by using the numerical solutions for r{z; 2;*). Then the value at the redshift 
for which P = 0.8 has been used as initial data for equation (6.12). For comparison, two 
horizontal lines have been added showing what temperatures would be required to obtain 
neutral fractions of 10^^ and 10^''^ in equilibrium, using equation (6.13). Since the plasma 
is in fact out of equilibrium, these highly ionized states can be maintained at much lower 
temperatures, as is seen in the bottom plots. 

6.3.4 IGM ionization and the Gunn-Peterson effect 

Assuming that any neutral hydrogen in the remains of the shells will have insufficient time 
to difi"usc far into the hot ionized regions that used to be shell interiors, we can treat the 
latter as an isolated mixture of gas and plasma where the ionization fraction x evolves as 

X = nx [(1 - X)^ioniT) - xKeciT)] , 

and where the rates for collisional ionization and recombination are given by (Stebbins &; 
Silk 1986) 

\meJ 



^rec — {(^recV) 



647r 4 2 fkT\~'^/^ 
3v37r V^-o, 



where oq is the Bohr radius and a is the fine structure constant. Changing the independent 
variable to redshift, the i7 = 1 case leaves us with 



116 



CHAPTER 6: LATE REIONIZATION.. 



XioniT) = n^toAioniT) « 5.7 x lO^hUb 



2.1/2 gl.58/T5 



Xrec{T) = ncoioArec(T) « [OAGhflh] T5 

For large enough z, the ionization fraction will adjust rapidly enough to remain in a qua- 
sistatic equilibrium and hence be given by x = 0, i.e. 



X 



1 + 2.8 X lO-^Ts'^/^^e^-^^/^^l ^ . (6.13) 



The observed absence of a Gunn-Petcrson trough in the spectra of high-redshift quasars 
strongly constrains the density of neutral hydrogen in the IGM. The most thorough study 
to date, involving eight quasars (Steidel &: Sargent 1987), concluded that 

nHj{z = 2.64) < (1.2 ± 3.1) X W^h^Q 

if Q = 1. In our model this corresponds to (1 - x) < (1-2 ± 3.1) x 10'^ /{f^Qb) ~ 2 x 10"^ 
for = 0.06 and /m,=0.1. Thus we are helped not only by the IGM being ionized, but also 
by it being diffuse. In a recent study of a single quasar, Webb et al. (1992) find the data 
consistent with cither VLhj{z = 4.1) = or VtHjiz = 4.1) = 1.5 x lO^^/igQ^, depending on 
model assumptions. We will use the latter value as an upper limit. Finally, recent Hubble 
Space Telescope spectroscopy of 3C 273 has been used to infer that ^Hj{z = 0.158) < 
1.4 X lQ~'^h'^Q. The constraints from these three studies are plotted in Figures 6.5 and 6.6 
together with the ionization levels predicted by our scenario. 

To achieve x = 10~^, 10~^ and 10~^ in equilibrium would by equation (6.13) require 
r > 5.5 X 10^ T > 1.1 X IQ^K and T > 3.6 x l{pK, respectively. As can be seen 
from the numerical solutions in the bottom panels of Figures 6.5 and 6.6, the recombination 
rate is generally too slow for equilibrium to be established, and the IGM remains almost 
completely ionized even when T <C 15, OOOif and equilibrium would have yielded x 0. 
In both 3a and 3b, a very moderate reheating (/j„j = 0.001, heavy lines) is seen to suffice 
to satisfy the three observational constraints. In the absence of any reheating whatsoever, 
the only models that satisfy the constraints are those with very low density (fi^ = 0.01 or 

/m = 0.01). 

In summary, the only parameters that are strongly constrained by the Gunn-Peterson 
test are Mc and c7o(Mc). 

6.3.5 Other spectral constraints 

Let us estimate to what extent Compton cooling of the hot plasma will distort the cosmic 
microwave background radiation (GBR). Since for Tg S> T^ the Gomptonization y-parameter 
(Stebbins & Silk 1986) 

kTe 
It rUeC^ 



kTe 

yC = / lyUeCFtCdt 



is linear in the plasma energy density ^| + |j kTeUe at each fixed time (Tg = Tjgm), all 
that counts is the spatially averaged thermal energy density at each redshift. Since the 
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former is simply et{z;Z:^) times the density of injected energy fgfsrS^bPC^-, the calculation 
reduces to mere energetics and we obtain 

d{—z*) 

where 



yc = y* / Vl + z et{z;z^)dzdz^, 

Jz d{-z^) Jo 



Tt = —fin^b — 7T ~ o.ombh. 

Svr m„G 



y* = Ifsn^l^ « 9 X w-^hnl 

The current observational upper limit on y is 2.5 X 10-5 (Mather et al. 1994), so even if we 
take fg{0) as high as 100% and make a gross overestimate of the integral by making all our 
galaxies as early as at = 30 and by replacing et{z; z^) by its upper bound 60% for all z, 
z*, our y is below the observational limit by three orders of magnitude for J7 = 1. 

Now let us estimate the optical depth of the IGM. It has long been known that reioniza- 
tion can cause a spatial smoothing of the microwave background as CBR photons Thomson 
scatter off of free electrons. Since Ue = XiGMnb, the optical depth for Thomson scattering, 
i.e. the number of mean free paths that a CBR photon has traversed when it reaches our 
detectors, is 

'^t= (TtXiGM-nbcdt = Tt / \/l + zxiGMdz, 

J tree ''0 

where 

3 ^ ^ Hocat 
ipC 

Let us evaluate the integral by making the approximation that XIGM increases abruptly from 
to 1 at some redshift Zion- Then even for Zjon as high as 30, Tt ~ 7.9hQbfm ^ 0.02 <§; 1 
for our fiducial parameter values h = 0.5, ilfe = 0.06 and fm = 0.1, so the probability that 
a given CBR photon is never scattered at all is e~^'^'^ ~ 98%. Hence this scenario for late 
reionization will have only a very marginal smoothing effect on the CBR. If the shells are 
totally ionized as well, then the factor fm disappears from the expressions above which helps 
only slightly. Then Zion = 15 would imply that 8% of the CBR would be spatially smoothed 
on scales of a few degrees. 

6.4 Discussion 

We have calculated the effects of supernova driven winds from early galaxies assuming a 
Press-Schechter model of galaxy formation and a CDM power spectrum. The calculations 
have shown that reionization by such winds can indeed explain the observed absence of a 
Gunn-Peterson effect if a number of conditions are satisfied: 

1. The masses of the first generation of galaxies must be very small, not greater than 
about 10^ Mq. 

2. There is enough power on these small scales to get at least 10% of the baryons in 
galaxies by z = 5. 

3. Except for the case where Qh is as low as 0.01, there must be some reheating of the 
IGM after z = 5 to prevent the IGM from recombining beyond allowed levels. 
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4. The commonly used thin-shell approximation for expanding bubbles must remain valid 
over cosmological timescales, with the mass fraction in the interior remaining much 
less than unity. 

Whether 1) is satisfied or not depends crucially on the model for structure formation. 
This scenario is consistent with a pure CDM model and some low-bias tilted CDM models, 
but not with top-down models like pure HDM. 

Observations of nearly solar abundances of heavy elements in intracluster gas have given 
some support for 2), which is roughly equivalent to requiring that at least 10% of the heavy 
elements in the universe be made before z = 5 (or whenever ^ 1). As discussed in 
Section 6.3.2, the observations of some extremely metal-poor objects in QSO absorption 
line studies do not necessarily rule out our scenario, since it is highly uncertain whether 
all the hydrogen in the swept-up IGM would get thoroughly mixed with the metal-rich 
supernova ejecta. The fact that large numbers of mini-galaxies are not seen today need not 
be a problem either. Possible explanations for this range from mechanisms for physically 
destroying them (Dekel & Silk 1986, for instance) to the fact that the faint end of the 
luminosity function is still so poorly known that old dwarf galaxies in the field may be too 
faint to see by the present epoch (Binggeli et al. 1988). 

To violate 3), the actual reheating would have to be extremely small. A current IGM 
temperature between lO'^K and lO^K suffices, depending on other parameter values, since 
the low density IGM never has time to reach its equilibrium ionization. 

The thin-shell approximation 4) is obviously a weak point in the analysis, because of 
the simplistic treatment of the dense shell and its interface with the interior bubble. For 
instance, could the shell cool and fragment due to gravitational instability before it collides 
with other shells? An approximate analytic model for such instability has been provided by 
Ostriker Sz Cowie (1981). Their criterion is that instability sets in when S > 1, where 

^ _ ^GpshellR^ 

Rvs 

and the sound speed Vg = ^bhT/Snip. In terms of our dimensionless variables, this becomes 

^ ^ V0.06; VO.i; \l + zj r'' 

which indicates that with our standard parameter values, gravitational instability does not 
pose problems even with fairly low shell temperatures. The reason that the shell density is 
not limited to four times the ambient IGM density is that the jump condition is not adiabatic, 
due mainly to effective Compton cooling at the high redshifts under consideration. 

After the critical z (typically between 20 and 5) at which the expanding shells have 
collided with neighbors and occupied most of space, the IGM is "frothy" on scales around 100 
kpc, with dense cool shell remnants scattered in a hot thin and fairly uniform plasma. Since 
the dark matter distribution is left almost unaffected by the expanding bubbles, formation 
of larger structures such as the galaxies we observe today should remain fairly unaffected as 
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far as concerns gravitational instability. There is indirect influence, however: the ubiquitous 
metals created by the early mini-galaxies would enhance the ability of the IGM to cool, which 
as mentioned in Section 6.3.1 is commonly believed to be crucial for galaxy formation. 

Blanchard et al. (1992) argue that if the IGM has a temperature higher than the virial 
temperature of a dark halo, pressure support will prevent it from falling into the potential 
well and thus stop it from forming a luminous galaxy. The virial temperature they estimate 
for an object of mass M formed at a redshift z is approximately 

T„,. 5.7X10' A- (j^) + 

for h = 0.5, so requiring Tyir > Tjgm for say Tjgm = at z = 5 would give a minimum 

galaxy mass of about 10^^ M©. Such arguments indicate that the IGM reheating of our 
scenario might produce a "mass desert" between the earliest mini-galaxies and the galaxies 
we see today: The first generation of galaxies, mini-galaxies with masses of perhaps 10^ or 
10^ Mq, would keep forming until their expanding bubbles had occupied most of space and 
altered the bulk properties of the IGM. After that, formation of galaxies much smaller than 
than those of today would be suppressed, since the IGM would be too hot. Eventually, 
as the IGM cools by adiabatic expansion, a progressively larger fraction of the IGM can 
be accreted by dark matter potential wells. Indeed, even with the volume averaged IGM 
temperature remaining hot due to some form of reheating, cooling flows in deep potential 
wells, in particular galaxy clusters, would not be suppressed. Late formation of galaxies is 
therefore possible. 

Pressure balance between the shell and the interior during the expansion would give the 
ration TghM/Tinterior = Pb,interior/pb,shell ~ « 0.03 for 6 = fm = 0.1, SO the sheU 

fragments would expected to contain non-negligible fractions of neutral hydrogen and thus 
absorb some Lyman-alpha. A typical shell radius is about 100 kpc for Mc = 10^ M©, a size 
comparable to that of the clouds of the Lyman-alpha forest. As to the number density of 
Lyman-alpha clouds, the observed velocity separations greatly exceed those we would expect 
if all shell fragments were to be identifled with Lyman alpha clouds. Thus the majority of 
these fragments must have been destroyed by some other process. There are a number of 
ways in which this could occur, for instance through photoionization by UV flux from the 
parent galaxy or by collapse to form other dwarf galaxies. The resulting numbers resemble 
the abundance of minihalos in an alternative explanation of the Lyman-alpha forest (Rees 
1986). Strong evolution, in the sense of an increasing cloud abundance with decreasing 
redshift, is expected to occur as cooling becomes effective. 
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Redshift 



Figure 6.1: Comoving radius of expanding shell. 

The comoving shell radius (1 + z)R is plotted for galaxies of total 
mass 2 x 10^ Mq, forming at integer redshifts from 1 to 29. Here 
O = 1, Jlft = 0.06, h = 0.5, and f„i = 0.1. = 1 for the upper set 
of lines and fd = for the lower set. R has been truncated when 
T drops below 15,000 K, after which newly swept up IGM fails to 
become ionized. 
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Figure 6.2: Energetics of expanding shell, example 1. 

This and the two following figures show the energy contents of an 
expanding bubble as a function of redshift, for different choices of 
fd and Zif. Q = 1, = 0.06, h = 0.5 and fm = 0.1 for all three 
plots. Figures 6.2 and 6.3 illustrate the difference between fd = and 
fd = i (there is no shock cooling loss in the second case). Figure 6.4 
has fd = and illustrates that the Compton cooling loss is larger at 
higher redshift. 
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Figure 6.3: Energetics of expanding shell, example 2. 
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Figure 6.4: Energetics of expanding shell, example 3. 
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Figures 6.5 and 6.6: IGM evolution. 

Three different properties of the IGM (fiUing factor, temperature and 
ionization) are plotted as a function of redshift for different choices 
of Mc, fg{5) and f^. n = I, Qf, = 0.06, h = 0.5 and /„ = 0.1 in the 
two previous figures, 6.5 and 6.6. In all panels, the different families 
of curves correspond to different values of /g(5); 50%, 20%, 10% and 
1% from left to right, with the rightmost cases being omitted where 
they fail dismally. In the porosity and temperature plots (the upper 
two panels of 6.5 and 6.6), dashed lines correspond to = 1 and 
solid ones to fd = 0, whereas only the pessimistic fd = case is 
plotted in the ionization plots (the lower third). In the temperature 
and ionization plots (the lower two panels), the three branches of 
each curve correspond to the three reheating scenarios: finj = 0, 
finj = 0.001 and finj = 0.02. 



Chapter 7 

Power Spectrum Independent 
Constraints 

In this chapter, a formahsm is presented that allows cosmological experiments to be tested for 
consistency, and allows a simple frequentist interpretation of the resulting significance levels. 
As an example of an application, this formalism is Tiscd to place constraints on bulk flows of 
galaxies using the results of the microwave background anisotropy experiments COBE and 
SP91, and a few simplifying approximations about the experimental window functions. It is 
found that if taken at face value, with the quoted errors, the recent detection by Lauer and 
Postman of a bulk flow of 689 km/s on scales of 150ft.^^Mpc is inconsistent with SP91 at a 
95% confidence level within the framework of a Cold Dark Matter (CDM) model. The same 
consistency test is also used to place constraints that arc completely model-independent, in 
the sense that they hold for any power spectrum whatsoever — the only assumption being 
that the random fields are Gaussian. It is shown that the resulting infinite-dimensional 
optimization problem reduces to a set of coupled non-linear equations that can readily be 
solved numerically. Applying this technique to the above-mentioned example, we find that 
the Lauer and Postman result is inconsistent with SP91 even if no assumptions whatsoever 
are made about the power spectrum. 

7.1 Introduction 

Together with the classical cosmological parameters h, CI, etc., the power spectrum P{k) of 
cosmological density fluctuations is one of the most sought-for quantities in modern cosmol- 
ogy, vital for understanding both the formation of large-scale structure and the fluctuations 
in the cosmic microwave background radiation (CBR). 

The traditional approach has been to assume some functional form for P{k) (like that 
predicted by the cold dark matter (CDM) scenario, for instance), and then investigate 
whether the predictions of the model are consistent with experimental data or not. The large 
amounts of data currently being produced by new CBR experiments and galaxy surveys, all 
probing different parts of the power spectrum, allow a new and more attractive approach. 
We can now begin to probe exact shape of the function P{k), without making any prior 
assumptions about P{k). More speciflcally, we measure different weighted averages of the 
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function, the weights being the experimental window functions. 

This new approach is quite timely (Juszkiewicz 1993), as there are now many indications 
that the primordial power spectrum may have been more complicated than an n = 1 power 
law. There are several sources of concern about the standard CDM cosmology, with inflation 
leading to ~ 1 and a primordial n ~ 1 Harrison-Zel'dovich power spectrum. Compared to 
COBE-normalized CDM, observational data shows unexpected large-scale bulk flows (Lauer 
& Postman 1994), too weak density correlations on small scales (Maddox et al. 1990), a 
rather quiet local velocity field (Schlcgcl et al. 1993) and a deficit of hot x-ray clusters 
(Oukbir & Blanchard 1992). The combined data from the COBE DMR (Smoot et al. 
1992) and the Tcnerife anisotropy experiment (Hancock et al. 1994) point to a spectral 
index exceeding unity (Watson &; Gutierrez de la Cruz 1993) which, if correct, cannot be 
explained by any of the standard inflationary models. The recent possible detections of 
halo gravitational microlensing events (Alcock et al. 1993) give increased credibility to the 
possibility that the dark matter in our galactic halo may be baryonic. If this is indeed 
the case, models with J7 < 1 and nothing but baryonic dark matter (BDM) (Peebles 1987; 
Gnedin & Ostriker 1992, Cen, Ostriker & Peebles 1993) become rather appealing. However, 
in contrast to CDM with inflation, BDM models do not include a physical mechanism that 
makes a unique prediction for what the primeval power spectrum should be. Rather, the 
commonly assumed P{k) oc k~^/'^ is chosen ad hoc to fit observational data. Moreover, for 
fluctuations near the curvature scale in open universes, where the Q = 1 Fourier modes are 
replaced by hyperspherical Bessel functions with the curvature radius as a built-in length 
scale, the whole notion of scale-invariance loses its meaning (Kamionkowski &: Spergel 1993). 

In summary, it may be advisable to avoid theoretical prejudice as to the shape of the 
primordial power spectrum. In this spirit, we will develop a consistency test that requires 
no such assumptions whatsoever about the form of the power spectrum. This approach was 
pioneered by Juszkiewicz; Gorski and Silk (1987), who developed a formalism for comparing 
two experiments in a power-spectrum independent manner. We generalize this method to 
the case of more than two experiments, and then use the formalism to assess the consistency 
of three recent observational results: the CBR anisotropy measurements made by the COBE 
Differential Microwave Radiometer (Smoot et al. 1992), the South Pole anisotropy experi- 
ment (SP91, Gaier et al. 1992), and the measurement of bulk velocity of Abell clusters in 
a 150 /i~^Mpc sphere (Lauer &; Postman 1994, hereafter LP). 

In Section 7.2, we develop a formalism for testing cosmological models for consistency. 
In Section 7.3, we apply this formalism to the special case of cold dark matter (CDM) and 
the LP, SP91 and COBE experiments. In Section 7.4, we solve the variational problem 
that arises in consistency tests of models where we allow arbitrary power-spectra, and apply 
these results to the LP, SP91 and COBE experiments. Section 7.5 contains a discussion of 
our results. Finally, two different goodncss-of-fit parameters arc compared in Appendix A, 
and the relevant window functions are derived in Appendix B. 
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7.2 Consistency Tests for Cosmological Models 

In cosmology, a field where error bars tend to be large, conclusions can depend crucially 
on the probabilistic interpretation of confidence limits. Confusion has sometimes arisen 
from the fact that large-scale measurements of microwave background anisotropics and bulk 
flows are fraught with two quite distinct sources of statistical uncertainty, usually termed 
experimental noise and cosmic variance. In this section, we present a detailed prescription for 
testing any model for consistency with experiments, and discuss the appropriate probabilistic 
interpretation of this test. By model we will mean not merely a model for the underlying 
physics, which predicts the physical quantities that we wish to measure, but also a model for 
the various experiments. Such a model is allowed to contain any number of free parameters. 
In subsequent sections, we give examples of both a very narrow class of models (standard 
CDM where the only free parameter is the overall normalization of the power spectrum), 
and a wider class of models (gravitational instability with Gaussian adiabatic fluctuations 
in a flat universe with the standard recombination history, the power spectrum being an 
arbitrary function). 

Suppose that we are interested in N physical quantities ci, ...,cn, and have N exper- 
iments Ei,...,Ei\f devised such that the experiment Ei measures the quantity Cj. Let Sj 
denote the number actually obtained by the experiment Ei. Because of experimental noise, 
cosmic variance, etc., we do not expect Sj to exactly equal q. Rather, Sj is a random vari- 
able that will yield different values each time the experiment is repeated. By repeating 
the experiment M times on this planet and averaging the results, the uncertainty due to 
experimental noise can be reduced by a factor v^M. However, if the same experiment were 
carried out in a number of different horizon volumes throughout the universe (or, if we have 
ergodicity, in an ensemble of universes with different realizations of the underlying random 
fleld), the results would also be expected to differ. This second source of uncertainty is 
known as cosmic variance. We will treat both of these uncertainties together by simply 
requiring the model to specify the probability distribution for the random variables Sj. 

Let us assume that the random variables Sj are all independent, so that the joint prob- 
ability distribution is simply the product of the individual probability distributions, which 
we will denote fi{s). This is an excellent approximation for the microwave background and 
bulk flow experiments we will consider. Finally, let si, ...,sjv denote the numbers actually 
obtained in one realization of the experiments. 

The general procedure for statistical testing will be as follows: 

• First, define a parameter r] that is some sort of measure of how well the observed data 
Si agree with the probability distributions fi, with higher r] corresponding to a better 
fit. 

• Then compute the probability distribution fr){r]) of this parameter, either analytically 
or by employing Monte-Carlo techniques. 

• Compute the observed value of rj, which we will denote fj. 
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• Finally, compute the probability P{ri < fj), i.e. the probability of getting as bad 
agreement as we do or worse. 

We will now discuss these four steps in more detail. 

7.2.1 Choosing a goodness-of-fit parameter 

Obviously, the ability to reject models at a high level of significance depends crucially on 
making a good choice of goodness-of-fit parameter 77. In the literature, a common choice is 
the likelihood product, i. e. 

N 

1=1 

In this chapter, we will instead use the probability product, i.e. the product of the probabil- 
ities Pi that each of the experiments yield results at least as extreme as observed. Thus if 
the observed Si is smaller than the median of the distribution /j, we have Pi = 2P{si < Si), 
whereas Si larger than the median would give Pi = 2P{si > Si). The factor of two is present 
because we want a two-sided test. Thus Pj = 1 if Sj equals the median. Pi = 2% if Sj is at 
the high 99th percentile, etc. 

These two goodness-of-fit parameters are compared in Appendix A for a few explicit 
examples, and the conclusions can be summarized as follows: 

1. Rather unphysical probability distributions can be concocted that "fool" either one of 
these two parameters but not the other. Thus neither the likelihood product nor the 
probability product can be hailed as fundamentally better than the other. 

2. For probability distributions encountered in the type of cosmology applications dis- 
cussed in this chapter, always smooth and unimodal functions, the likelihood product 
and the probability product yield very similar results. 

3. As described in the following sections, the probability distribution of the probability 
product can always be computed analytically. The probability distribution for the 
likelihood product, however, depends on the distributions of the underlying random 
variables, and except in a few fortuitous simple cases, it must be computed numerically 
through either repeated convolutions or Monte Carlo simulations. Thus the probability 
product is considerably more convenient to use. 

Since 1) and 2) indicate that neither of the two goodness-of-fit parameters is preferable over 
the other on scientific grounds, the authors feel that 3) tips the balance in favor of the 
probability product. 

7.2.2 Its probability distribution 

Apart from the simple interpretation of the probability product rj, it has the advantage that 
its probability distribution can be calculated analytically, and is completely independent 
of the physics of the model — in fact, it depends only on A'^. We will now give the exact 
distributions. 
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Confidence level 


A^ = l N = 2 N = 3 N = 4: 


95% 
99% 
99.9% 


0.05 0.0087 0.0018 0.00043 
0.01 0.0013 0.00022 0.000043 
0.001 0.000098 0.000013 0.0000021 



Table 7.1: Probability product limits 



By construction, < rj < 1. For N = 1, rj will simply have a uniform distribution, i.e. 

f^(r,) = l^ ifO<7?<l, 
^ 1 otherwise. 

Thus in the general case, rj will be a product of N independent uniformly distributed random 
variables. The calculation of the probability distribution for r] is straightforward, and can 
be found in a number of standard texts. The result is 

/.(,) = -/.(-in,)^-!!^'-'"")"- 

"'''1 to otherwise. 

7.2.3 The consistency probability 

The probability P{r] < fj), the probability of getting as bad agreement as we do or worse, 
is simply the cumulative distribution function Fj^^r)), and the integral can be carried out 
analytically for any N: 

F,{fi) ^ Piv <V)= r Uu)du = f,e{f,) ^ t}^^ (7.1) 

where 9 is the Heaviside step function, and i^jj(^) = 1 for f} > 1. Since the product of N 
numbers between zero and one tends to zero as iV — ^ oo, it is no surprise that 

Fr,{fi) ^ e{f,)fie-''''' = e{fi) 

as AT — DO, i.e. that /x(?7) — ^{v)- The function is plotted in Figure 7.1, and the 

values of r) for which = 0.05, 0.01 and 0.001, respectively, are given in Table 7.1 for 

a few AT-values. For example, if three experimental results give a goodness-of-fit parameter 
fj = 0.0002 for some model, then this model is ruled out at a confidence level of 99%. Thus 
if the model where true and the experiments where repeated in very many different horizon 
volumes of the universe, such a low goodness-of-fit value would be obtained less than 1% of 
the time. 

7.2.4 Ruling out whole classes of models 

If we wish to use the above formalism to test a whole set of models, then we need to solve an 
optimization problem to find the one model in the set for which the consistency probability 
is maximized. For instance, if the family of models under consideration is standard n = 1, 
r = 0.5 CDM (sec Section 7.3), then the only free parameter is the overall normalization 
constant A. Thus we can write the consistency probability as p{A), and use some numerical 
method to find the normalization A^ for which p{A) is maximized. After this, the statistical 
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interpretation is clear: if the experiments under consideration are carried out in an ensemble 
of CDM universes, as extreme results as those observed will only be obtained at most a 
fraction p{A^) of the time, whatever the true normalization constant is. Precisely this case 
will be treated in the next section. For the slightly wider class of models consisting of CDM 
power spectra with arbitrary A, n and F, the resulting optimization problem would be a 
three-dimensional one, and the maximal consistency probability would necessarily satisfy 

p{A^,n^,T^) >p(A*, 1,0.5) =p{A^). 

An even more general class of models is the set of all models where the random fields are 
Gaussian, i.e. allowing completely arbitrary power spectra P. In section 4, we will show that 
the resulting infinite-dimensional optimization problem can in be reduced to a succession of 
two finite-dimensional ones. 



7.3 Cold Dark Matter Confronts SP91, COBE and Lauer- 
Postman 

As an example of an application of the formalism presented in the previous section, we will 
now test the standard cold dark matter (CDM) model of structure formation for consistency 
with the SP91 CBR experiment and the Lauer-Postman bulk flow experiment. 

Let El be the Lauer-Postman (LP for short) measurement of bulk flows of galaxies in a 
150/i""^Mpc sphere (Lauerfe Postman, 1994). Let E2 be the 1991 South Pole CBR anisotropy 
experiment, SP91 for short (Gaier et al. 1992). Let E3 be the COBE DMR experiment 
(Smoot et al. 1992). All of these experiments probe scales that are well described by linear 
perturbation theory, and so as long as the initial fluctuation are Gaussian, the expected 
results of the experiments can be expressed simply as integrals over the power spectrum of 
the matter perturbation: 

{si) = j Wi{k)P{k)dk. 

Here Sgp and Sco6e are the mean-square tempcratTire fluctuations measured by the exper- 
iments, and Sip = {v/cf' is the squared bulk flow. The corresponding window functions 
^ip^ ^sp and Wcobe are derived in Appendix B, and plotted in Figure 7.2. These window 
functions assume that the initial perturbations were adiabatic, that $7 = 1, and that re- 
combination happened in the standard way, x.e. a last-scattering surface at 2; 1000. The 
SP91 window function is to be interpreted as a lower limit to the true window function, 
as it includes contributions only from the Sachs- Wolfe effect, not from Doppler motions or 
intrinsic density fluctuations of the surface of last scattering. 

Now let us turn to the probability distributions for the random variables sip, Sgp and 
Sc- The standard CDM model with power-law initial fluctuations oc A;" predicts a power 
spectrum that is well fltted by (Bond k, Efstathiou 1984) 



(l + [ag + (6(?)i-5 + (cg)2]i-i3) 
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n 


r 


LP 




SP91 


COBE 




LP/SP91 


LP/COBE 


1 


0.5 


9.2 X 10" 


-7 


n — 

1.6 X 10"** 


2.0 X 10" 


-8 


56.7 


45.1 


0.7 


0.5 


1.7 X 10" 


-6 


2.9 X 10"^ 


5.1 X 10" 


-8 


57.2 


32.7 


2 


0.5 


1.4 X 10" 


_7 


2.6 X 10~^ 


1.2 X 10" 


.9 


53.9 


112.2 


1 


0.1 


1.4 X 10" 


-6 


2.3 X 10-^ 


4.3 X 10- 


-8 


57.9 


31.9 


1 


10 


2.3 X 10- 


-7 


4.1 X 10-9 


4.6 X 10" 


-9 


55.9 


49.6 



Table 7.2: Expected r.m.s. signals for CDM power spectrum with A = (Ih ^Mpc)^ 

where a = 6.4, b = 3.0, c = 1.7 and q = (lh~^Mpc)k /T . For the simplest model, T = h, but 
certain additional complications such as a non-zero cosmological constant A and a non-zero 
fraction fli, of hot dark matter can be fitted with reasonable accuracy by other values of T 
(Efstathiou, Bond & White 1992). Thus the model has three free parameters: n, T and the 
overall normalization A. Integrating the power spectrum against the three window functions 
yields the values of Cj given in Table 7.2. The two rightmost columns contain the quotients 
cip/cgp and cip/cc, respectively. As can be seen, the dependence on F is quite weak, and the 
quotient cip/cgp is quite insensitive to the spectral index n as well. Let us for definiteness 
assume the canonical values n = 1 and F = 0.5 in what follows. 

These values q would be the average values of the probability distributions for Sgp and 
Sip if there where no experimental noise. We will now model the full probability distributions 
of the three experiments, including the contribution from experimental noise. 

For a bulk flow experiment, the three components Vx, Vy and Vz of the velocity vector v 
are expected to be independent Gaussian random variables with zero mean, and 

(|v|'>=q^. 

However, this is not quite the random variable sip that we measure, because of errors in 
distance estimation, etc. Denoting the difference between the observed and true bulk velocity 

vectors by e, let us assume that the three components of e arc identically distributed and 
independent Gaussian random variables. This should be a good approximation, since even 
if the errors for individual galaxies are not, the errors in the average velocity e will be 
approximately Gaussian by the Central Limit Theorem. Thus the velocity vector that we 
measure, v + e, is also Gaussian, being the sum of two Gaussians. The variable that we 
actually measure is s^p = |v -|- ep, so 

sip = ^{cip + Vip)xl, 

where has a chi-squared distribution with three degrees of freedom, and Vip is the variance 
due to experimental noise, i.e. the average variance that would be detected even if the true 
power spectrum were P{k) = 0. The fact that the expectation value of the detected signal 
Sip (which is usually referred to as the uncorrected signal in the literature) exceeds the true 
signal cip is usually referred to as error bias (LP; Strauss, Ccn & Ostriker 1993 — hereafter 
SCO). Error bias is ubiquitous to all experiments of the type discussed in this chapter, 
including GBR experiments, since the measured quantity is positive definite and the noise 



134 



CHAPTER 7: POWER SPECTRUM INDEP. CONSTRAINTS 



errors contribute squared. In the literature, experimentally detected signals are usually 
quoted after error bias has been corrected for, i.e. after the noise has been subtracted from 
the uncorrected signal in For LP, the uncorrected signal is 807 km/s, whereas the signal 
quoted after error bias correction is 689 km/s. 

For the special case of the LP experiment, detailed probability distributions have been 
computed using Monte-Carlo simulations (LP, SCO), which incorporate such experiment- 
specific complications as sampling errors, asymmetry in the error ellipsoid, etc. To be used 
here, such simulations would need to be carried out for each value of c/p under consideration. 
Since the purpose of this section is merely to give an example of the test formalism, the 
above-mentioned x^-approximation will be quite sufficient for our needs. 

For the SP91 nine-point scan, the nine true values ATi/T arc expected to be Gaussian 
random variables that to a good approximation are independent. They have zero mean, and 

{\ATi/Tf) = csp. 

Denoting the difference between the actual and observed values by 5i , we make the standard 
assumption that these nine quantities are identically distributed and independent Gaussian 
random variables. Thus the temperature fluctuation that we measure at each point, ATj/T+ 
(5, is again Gaussian, being the sum of two Gaussians. The variable that we actually measure 
is 




where Xg has a chi-squared distribution with nine degrees of freedom, and Vsp is the variance 
due to experimental noise, the error bias, i.e. the average variance that would be detected 
even if the true power spectrum were P{k) = 0. 

We will use only the signal from highest of the four frequency channels, which is the 
one likely to be the least affected by galactic contamination. Again, although Monte-Carlo 
simulations would be needed to obtain the exact probability distributions, we will use the 
simple x^-approximation here. In this case, the main experiment-specific complication is 
the reported gradient removal, which is a non-linear operation and thus does not simply 
lead to a x^-distribution with fewer degrees of freedom. 

The amplitude of the COBE signal can be characterized by the variance in AT/T on 
an angular scale of 10°. This number can be estimated from the COBE data set as Sc = 
cr^Qo = ((11.0±1.8) X 10~^)^ (Smoot et al. 1992). The uncertainty in this quantity is purely 
due to instrument noise, and contains no allowance for cosmic variance. We must fold in 
the contribution due to cosmic variance in order to determine the probability distribution 
for Sc- We determined this probability distribution by performing Monte-Carlo simulations 
of the COBE experiment. We made simulated COBE maps with a variety of power spectra 
(including power laws with indices ranging from to 3, as well as delta-function power 
spectra of the sort described in Section 7.2). We included instrumental noise in the maps, 
and excluded all points within 20° of the Galactic plane. By estimating Sc from each map, 
we were able to construct a probability distribution corresponding to each power spectrum. 
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In all cases, the first three moments of the distribution were well approximated by 

' IJ-l = {Sc) = Cc, 

< IJL2 = (s2) - (sc)^ < 0.063c2 + 1.44 X 10-21, (7.2) 
= (s^) = 0.009c3. 

Furthermore, in all cases the probability distributions were well modeled by chi-squared 
distributions with the number of degrees of freedom, mean, and offset chosen to reproduce 
these three moments. Note that the magnitude of the cosmic variance depends on the shape 
of the power spectrum as well as its amplitude. The inequality in the above expression 
for represents the largest cosmic variance of any of the power spectra we tested. Since 
we wish to set conservative limits on models, we will henceforth assume that the cosmic 
variance is given by this worst-case value. Thus we are assuming that the random variable 
(sc — sq)I/S.s has a chi-squared distribution with 5 degrees of freedom, where 

'so = Ml - 2M2//^3, 

< As =Ai3/4Ai2, (7.3) 

The results obtained using these three probability distributions are summarized in Ta- 
bles 7.3, 7.4 and 7.5. In 7.3 and 7.4, N = 2, and the question asked is whether LP is 
consistent with COBE and SP91, respectively. In Table 7.5, N = 3, and we test all three 
experiments for consistency simultaneously. In each case, the optimum normalization (pro- 
portional to the entries labeled "Signal") is different, chosen such that the consistency 
probability for the experiments under consideration is maximized. As can be seen, the last 
two tests rule out CDM at a significance level of 95%, i.e. predict that in an ensemble of 
universes, results as extreme as those we observe would be obtained less than 5% of the 
time. Note that using both COBE and SP91 to constrain LP yields a rejection that is no 
stronger than that obtained when ignoring COBE. In the latter case, the best fit is indeed 
that with no cosmological power at all, which agrees well with the observation of SCO that 
sampling variance would lead LP to detect a sizable bulk flow (before correcting for error 
bias) even if there where none. 

7.4 Allowing Arbitrary Power Spectra 

In this section, we will derive the mathematical formalism for testing results from multiple 
experiments for consistency, without making any assumptions whatsoever about the power 
spectrum. This approach was pioneered by Juszkiewicz et al. (1987) for the case N = 2. 
Here we generalize the results to the case of arbitrary N. Despite the fact that the original 
optimization problem is infinite-dimensional, the necessary calculations will be seen to be of 
a numerically straightforward type, the case of N independent constraints leading to nothing 
more involved than numerically solving a system of N coupled non-linear equations. After 
showing this, we will discuss some inequalities that provide both a good approximation of 
the exact results and a useful qualitative understanding of them. 



CHAPTER 7: POWER SPECTRUM INDEP. CONSTRAINTS 





TjP CORR Combined 


Noise 

Si PTi al 

Noise+Signal 
Detected 


420km/s 9.8 /xK 
169 km /s 33 8uK 
453km/s 35.2 /xK 
807km/s 35.2 /xK 


fj 

P{r] < fj) 


0.046 1.00 0.046 
0.046 1.00 0.19 



Table 7.3: Are LP and COBE consistent with CDM? 





LP SP91 Combined 


Noise 
Signal 

Noise+Signal 
Detected 


420km/s 26.4 

Okm/s 0/xK 
420km/s 26.4 /xK 
807km/s 19.9 /xK 


f] 

P{r] < fj) 


0.023 0.35 0.0079 
0.023 0.35 0.046 



Table 7.4: Are LP and SP91 consistent with CDM? 





LP SP91 COBE Combined 


Noise 
Signal 

Noise+Signal 
Detected 


420km/s 26.4 /xK 9.8 /xK 
168km/s 26.9 /xK 33.8 /xK 
452km/s 37.7 35.1 
807km/s 19.9^K 35.2 


fj 

P{r] < fj) 


0.046 0.039 0.97 0.0017 
0.046 0.039 0.97 0.046 



Table 7.5: Arc LP, SP91 and COBE ah consistent with CDM? 



Tables 7.3, 7.4 and 7.5 show the consistency probability calculations. The first 
line in each tabic gives the experimental noise, i.e. the detection that wonld be 
expected in the absence of any cosmological signal. The second line is the best- 
fit value for the cosmological signal c, the value that maximizes the combined 
consistency probability in the lower right corner of the table. The third line 
contains the expected value of an experimental detection, and is the sum in 
quadrature of the two preceding lines. The fourth line gives the goodness-of-fit 
parameter for each of the experiments, i.e. the probability that they would yield 
results at least as extreme as they did. The rightmost number is the combined 
goodness-of-fit parameter, which is the product of the others. The last line 
contains the consistency probabilities, the probabilities of obtaining goodness- 
of-fit parameters at least as low as those on the preceding line. 
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7.4.1 The optimization problem 

Let us consider N = n + 1 experiments numbered 0, 1, n that probe the cosmological 
power spectrum P{k). We will think of each experiment as measuring some weighted average 
of the power spectrum, and characterize an experiment Ei by its window function Wi{k) as 
before. 

Purely hypothetically, suppose we that we had repeated the same experiments in many 
different locations in the universe, and for all practical purposes knew the quantities ci, c„ 
exactly. Then for which power spectrum P{k) would cq be maximized, and what would this 
maximum be? If we experimentally determined cq to be larger than this maximum value, 
our results would be inconsistent, and we would be forced to conclude that something was 
fundamentally wrong cither with our theory or with one of the experiments. In this section, 
we will solve this hypothetical problem. After this, it will be seen that the real problem, 
including cosmic variance and experimental noise, can be solved in almost exactly the same 
way. 

The extremal power spectrum we are looking for is the solution to the following linear 
variational problem: 
Maximize 

/ P{k)Wo{k)dk (7.4) 
Jo 

subject to the constraints that 

/o~ P{k)Wi{k)dk = Ci for = 1, n, 
P{k) > for all A; > 0. 

This is the infinite-dimensional analogue of the so called linear programming problem, 
and its solution is quite analogous to the finite-dimensional case. In geometrical terms, we 
think of each power spectrum as a point in the infinite dimensional vector space of power 
spectra (tempered distributions on the positive real line, to be precise), and limit ourselves to 
the subset Q of points where all the above constraints are satisfied. We have a linear function 
on this space, and we seek the point within the subset Q where this function is maximized. 
We know that a differentiable functional on a bounded region takes its maximum either 
at an interior point, at which its gradient will vanish, or at a boundary point. In linear 
optimization problems like the one above, the gradient (here the variation with respect to 
P, which is simply the function Wq) is simply a constant, and will never vanish. Thus 
any maximum will always be attained at a boundary point. Moreover, from the theory 
of linear programming, we know that if there are n linear constraint equations, then the 
optimum point will be a point where all but at most n of the coordinates are zero. It is 
straightforward to generalize this result to our infinite-dimensional case, where each fixed k 
specifies a "coordinate" P{k), and the result is that the solution to the variational problem 
is of the form 

Pik) = J2pAk-ki). 

i=l 

This reduces the optimization problem from an infinite-dimensional one to a 2n-dimensional 
one, where only the constants pi and ki remain to be determined: 
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Maximize 

n 

J2PjWo{kj) (7.5) 

subject to the constraints that 

f Y.]=iPjWiikj) = Ci for z = 1, ...,n, 

[ Pi > for 1 = 1, n. 

This problem is readily solved using the method of Lagrange multipliers: defining the La- 
grangian 

n n n 

j=l 1=1 [j=l 

and requiring that all derivatives vanish leaves the following set of 3n equations to determine 
the 3n unknowns pi, ki and A^: 



Pi 



0, 

0, 
0. 



Introducing matrix notation by defining the fcj-dependent quantities Aij = Wj{ki), Bij = 
Wj{ki), (li = WQ{ki) and hi = Wo(A;j) brings out the structure of these equations more 
clearly: If pi ^ 0, then 

= a, 
B\ =b, 
. A^Tp = c. 

If A and B are invertible, then eliminating A from the first two equations yields the following 
system of n equations to be solved for the n unknowns k\, kn- 



A-^aL = B-^h. 



(7.6) 



Although this system is typically coupled and non-linear and out of reach of analytical 
solutions for realistic window functions, solving it numerically is quite straightforward. A 
useful feature is that once this system is solved, a, b, A and B are mere constants, and the 
other unknowns are simply given by matrix inversion: 

= A-^a 



A 



Since the non-linear system (7.6) may have more than one solution, all solutions should be 
substituted back into (7.5) to determine which one is the global maximum. Furthermore, to 
make statements about the solution to our original optimization problem (7.4), we need to 
consider also the case where one or more of the n variables pi,...,Pn vanish. If exactly m of 
them are non-vanishing, then without loss of generality, we may assume that these are the 
first m of the n variables. Thus we need to solve the maximization problem (7.5) separately 
for the cases where P{k) is composed of n delta functions, n—1 delta functions, etc., all the 
way down to the case where P{k) is single delta function. These solutions should then be 
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substituted back into (7.4) to determine which is the global maximum sought in our original 
problem. Thus the solutions depend on the window functions Wj and the signals Cj in the 
following way: 

• From the window functions alone, we can determine a discrete and usually finite 
number of candidate wavenumbers k where delta functions can be placed. 

• The actual signals Cj enter only in determining the coefficients of the delta functions 
in the sum, i.e. in determining what amount of power should be hidden at the various 
candidate wavenumbers. 

If we have found an the optimal solution, then a small change in the signal vector c will 
typically result in a small change in p and no change at all in the number of delta functions 
in P{k) or their location. If c is changed by a large enough amount, the delta functions 
may suddenly jump and/or change in number as a different solution of (7.5) takes over as 
global optimum or one of the coefficients pi becomes negative, the latter causing the local 
optimum to be rejected for constraint violation. Thus within certain limits, we get the 
extremely simple result that for the optimal power spectrum P{k), 



Thus within these limits, cq depends linearly on the observed signal strengths Cj. This is 
exactly analogous to what happens in linear programming problems. 

7.4.2 A useful inequality 

Before proceeding further, we will attempt to provide a more intuitive understanding of the 

results of the previous section, and show how to determine how complicated a calculation 
is justified. For the special case of only a single constraint, i.e. n = 1, we obtain simply 
P{k) = pi6{k — ki), where ki is given by 



Then we see that for n = 1, ki is simply the wavenumber for which the function fi is 
maximized, and that the maximum signal possible is simply cq = fi{ki). Thus the maximum 

signal in experiment that is consistent with the constraint from experiment i is obtained 
when the power is concentrated where the function fi is large. In other words, if we want 
to explain a high signal cq in the face of low signals in several constraining experiments, 
then the best place to hide the necessary power from the i*^ experiment is where fi takes 
its maximum. These functions are plotted in Figure 7.4 for the experiments discussed in 
Section 7.3, the optimization problem being the search for the maximum LP signal that is 
consistent with the constraints from SP91 and COBE. For illustrative purposes, we here 




w;,iki)Wi(ki) = wi{ki)Woiki). 



For the case of n constraints, let us define the functions 



Wojk) 
Wi{k) 
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assume that Csp and Cc are known exactly, and given by the detected signals Ssp ~ 19.9 /xK 
and §1^^ 33.8 /xK (we will give a proper treatment of cosmic variance and noise in the 
following section). Using the n = 1 constraint for each constraining experiment separately, 
the smallest of the functions thus sets an upper limit to the allowed signal cq — Csp- Thus 
the limit is given by the highest point in the hatched region in Figure 7.4, i. e. 

CO < c^^L = sup mm fi{k). 
k » 

Wc sec that using the SP91 constraint alone, the LP signal would be maximized if all power 

were at /c ~ (940Mpc)^^. Since this flagrantly violates the COBE constraint, the best place 

to hide the power is instead at ~ (lOOMpc)"^. 

By using the above formalism to impose all the constraints at once, the allowed signal 

obviously becomes lower. If the constraints are equalities rather than inequalities, then this 

stronger limit can never lie below the value at (fc* « 250Mpc)~^, where fsp{K) = fc{k*), 

since this is the signal that would result from a power spectrum of the form P{k) oc 5{k — k^). 

Thus for the particular window functions in our example, where the constraint from the 

n = 2 calculation cannot be more than a factor /sp(80Mpc)//sp(250Mpc) 1.05 stronger 

than the simple n = 1 limits, the latter are so close to the true optimum that they are quite 

sufficient for our purposes. If the constraints are upper limits rather than equalities, then 

the limit on cq is more relaxed, and is always the uppermost point in the hatched region, 
(1) 

I. e. Crnax ■ 

7.4.3 Including noise and cosmic variance 

To correctly handle cosmic variance and instrumental noise, we need to use the formalism 
developed in Section 7.2. Thus given the probability distributions for the various experi- 
mental results Si, we wish to find the power spectrum for which the consistency probability r] 
is maximized. This optimization problem, in which all experiments are treated on an equal 
footing, will be seen to lead directly to the asymmetric case above where the signal in one 
is maximized given constraints from the others. For definiteness, we will continue using the 
example with the LP, SP91 and COBE experiments. As seen in Section 7.2, the source of 
the low consistency probabilities is that sip is quite high when compared to Sgp and Sc- Thus 
it is fairly obvious that for the power spectrum that maximizes the consistency probability, 
we will have sip > {sip), whereas Ssp < {ssp) and Sc < (sc), so we can neglect power spectra 
that do not have this property. Let us first restrict ourselves to the subset of these power 
spectra for which cip = D and Cc = E, where D and E are some constants. Then these 
power spectra all predict the same probability distributions for sip and Sc- The consistency 
probability ry is clearly maximized by the power spectrum that maximizes (s^p), and this 
will be a linear combination of one or two delta functions as shown in Section 7.4.1. The 
key point is that since the locations of these delta functions are independent of D and E 
(within the range discussed in 4.1), the infinite- dimensional optimization problem reduces 
to the following two simple steps: 
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1. Solve for the optimal number of delta functions m and their locations ki as described 
in Section 7.4.1 

2. Find the m coefficients pi for which the power spectrum P{k) = J2i^iPi^{k ~ ^i) 
maximizes the consistency probability. 

7.4.4 Power spectrum independent constraints on LP, SP91 and COBE 

When applying the above consistency test to the LP, SP91 and COBE experiments, we 
obtain exactly the same consistency probability as in Table 7.4. The reason for this is that 
the optimal normalization turns out to be zero. This will obviously change if the LP error 
bars become smaller in the future. Thus dropping the CDM assumption does not improve 
the situation at all, which indicates that main source of the inconsistency must be something 
other than the CDM model. 

In anticipation of future developments, consistency probabilities were also computed for 
a number of cases with less noise in the LP experiment. Comparing only LP and SP91, the 
optimum power spectrum has a delta function aX k ^ (941Mpc)'~^. When including all three 
experiments, treating the COBE and SP91 constraints as upper limits, the optimum power 
spectrum has a single delta function at A; (79Mpc)~^, so the addition of COBE strengthens 
the constraint only slightly, due to the flatness of fgp in Figure 7.4. Interestingly, for all 
these cases with smaller LP error bars, consistency probabilities were found to be almost as 
low when allowing arbitrary power spectra as for the CDM case. This is again attributable 
to the flatness of fgp, since weighted averages of a flat function are fairly independent of the 
shape of the weight function (here the power spectrum). 

7.5 Discussion 

We have developed a formalism for testing multiple cosmological experiments for consistency. 
As an example of an application, we have used it to place constraints on bulk flows of 
galaxies using the COBE and SP91 measurements of fluctuations in the cosmic microwave 
background. It was found that taken at face value, the recent detection by Lauer and 
Postman of a bulk flow of 689 km/s on scales of 150/i~^Mpc is inconsistent with SP91 
within the framework of a CDM model, at a signiflcance level of about 95%. However, 
interestingly, this cannot be due solely to the CDM assumption, since the LP result was 
shown to be inconsistent with COBE and SP91 at the same significance level even when no 
assumptions whatsoever were made about the power spectrum. This leaves four possibilities: 

1. The window functions are not accurate. 

2. Something is wrong with the quoted signals or error bars for at least one of the 
experiments, 

3. The observed fluctuations cannot be explained within the framework of gravitational 
instability and the Sachs- Wolfe effect. 
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4. The random fields are not Gaussian, 

Case (1) could be attributed to a number of effects: If $7 7^ 1, then both the calculation of 
the Sachs- Wolfe effect (which determines Wgp and Wcobe) cind the growth of velocity pertur- 
bations (which determines Wip) are altered. If the universe became reionized early enough 
to rescatter a significant fraction of all CBR photons, then small scale CBR anisotropics 
were suppressed, which would lower Wgp. A quantitative treatment of these two cases will 
be postponed to future work. Other possible causes of (1) include a significant fraction of 
the density perturbations being isocurvature (entropy) perturbations or tensor-mode pertur- 
bations (gravity waves). Apart from these uncertainties, we have made several simplifying 
assumptions about the window functions for LP and SP91. To obtain more accurate consis- 
tency probabilities than those derived in this chapter, a more accurate LP window function 
should be used that incorporates the discreteness and the asymmetry of the sample of Abell 
clusters used. This can either be done analytically (Feldman &; Watkins 1993) or circum- 
vented altogether by performing Monte-Carlo simulations like those of LP or SCO, but for 
the whole family of power spectra under consideration. 

As to case (2), there has been considerable debate about both the LP and the SP91 
experiments. A recent Monte-Carlo Simulation of LP by SCO basically confirms the large 
error bars quoted by LP. As is evident from the flatness of the LP curve in Figure 7.3, it will 
be impossible to make very strong statements about inconsistency until future experiments 
produce smaller error bars. With the SP91 experiment, a source of concern is the validity 
of using only the highest of the four frequency channels to place limits, even though it is 
fairly clear that the other three channels suffer from problems with galactic contamination. 
The situation is made more disturbing by the fact that a measurement by the balloon- 
borne MAX experiment (Gundersen et al. 1993) has produced detections of degree-scale 
fluctuations that that are higher than those seen by SP91, and also higher than another 
MAX measurement (Meinhold et al. 1993). Other recent experiments that have detected 
greater fluctuations include ARGO (de Bernardis et al. 1993), PYTHON (Dragovan et al. 
1993) and MS AM (Cheng et al. 1993). On the other hand, SP91 has been used only as an 
upper limit in our treatment, by including only the Sachs- Wolfe effect and neglecting both 
Doppler contributions from peculiar motions of the surface of last scattering and intrinsic 
density fluctuations at the recombination epoch. If these effects (which unfortunately depend 
strongly on parameters such as h and Clb) where included, the resulting constraints would 
be stronger. 

Case (3) might be expected if the universe underwent a late-time phase transition, since 
this could generate new large-scale fluctuations in an entirely non-gravitational manner. 

In the light of the many caveats in categories (1) and (2), the apparent inconsistency 
between LP and SP91 (Jaffe et al. 1993) is hardly a source of major concern at the present 
time, and it does not appear necessary to invoke (3) or (4). However, we expect the testing 
formalism developed in this chapter to be able to provide many useful constraints in the 
future, as more experimental data is accumulated and error bars become smaller. 
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The cumulative probability distribution for the goodness-of-fit pa- 
rameter 77 is plotted for a few different n- values. 
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Figure 7.2: Window functions. 

The window functions of the Lauer-Postman bulk flow measurement 
(LP), the South Pole 1991 nine point scan (SP91), and the COBE 
DMR 10° pixel r.m.s. measurement (COBE) are plotted as a function 
of comoving wavenumber k. 
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True LP power [km/s] 



Figure 7.3: Consistency probabilities. 

The probability that the the Lauer-Postman bulk flow measurement 
(LP), the COBE DMR experiment and the South Pole 1991 nine 
point scan (SP91) are consistent with CDM is plotted as a function 
of the normalization of the power spectrum. The normalization is 
expressed in terms of the expected bulk flow in a LP measurement. 
The dashed line is the product of these three probabilities, and takes 
a maximum for a normalization corresponding to 168 km/s. 
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Figure 7.4: The best places to hide power. 

The functions fgp (soHd line) and fcobe (dashed line) are plotted as 
a function of wavenumber k. The shaded region, i.e. the area lying 
beneath both curves, constitutes the LP bulk flows that would be 
consistent with the SP91 and COBE experiments using the = 1 
constraints only, when the power spectrum is a single delta function 
located at k. 



Appendix A 

The Efficiency Parameter 



In this appendix, we will discuss the parameter fion, and see that it rarely drops below 
30%. All the notation is that of Chapter 3. We will first discuss the thermal evolution of 
intergalactic hydrogen exposed to a strong UV flux, and then use the results to write down a 
differential equation for the volume fraction of the universe that is ionized, subject to point 
sources of UV radiation that switch on at different times. We will see that photoionization 
is so efficient within the ionized regions of the IGM that quite a simple equation can be 
given for the expansion of the ionized regions. 

The evolution of IGM exposed to ionizing radiation has been discussed by many authors. 
Important early work includes that of Arons & McCray (1970), Bergeron & Salpeter (1970) 
and Arons &; Wingert (1972). The main novelty of the treatment that follows is that whereas 
previous treatments focus on late {z < 5) epochs, when various simplifying approximations 
can be made because the recombination and Compton rates are low, we are mainly interested 
in the case 50 < z < 150. Wc show that IGM exposed to a strong UV flux rapidly approaches 
a quasistatic equilibrium state, where it is almost fully ionized and the temperature is 
such that photoionization heating exactly balances Compton cooling. This simplifies the 
calculations dramatically, since the entire thermal history of the IGM can be summarized 
by a single function x{^)i the volume fraction that is ionized. Thus a fraction xi^) is ionized 
and hot (with a temperature that depends only on z, not on when it became ionized), and 
a fraction 1 — xi^) is neutral and cold. 

In the first section, we justify this approximation. In the second section, we derive a 
differential equation for the time-evolution of % as well as a useful analytic estimate of fion- 



A.l Intergalactic Stromgren Spheres 

Let X denote the ionization fraction in a small, homogeneous volume of intergalactic hy- 
drogen, i.e. X = uhii / {nni + nnii)- {x is not to be confused with x, the volume fraction 
in ionized bubbles.) When this IGM is at temperature T, exposed to a density of 77 UV 
photons per proton, the ionization fraction x evolves as follows: 



dx 



l + z 



d{-z) VI + 



\pi{l - x) + Xcix{l - x)- 



X 



(A.l) 
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where Hq^{1 + z) ^ times the rates per baryon for photoionization, coUisional ionization 
and recombination are given by 

f Xpi « 1.04 X 10^2 [hnbCTis] V, 
Xci « 2.03 X W^hClbTl^^e-'^^-^/^*' (A.2) 



A^ec ~ 0.717/iJ^6T4 11.808 - 0.51nT4 + O.I87T4' 

and T4 = r/lO^K. ais is the thermally averaged photoionization cross section in units of 
lO^^^cm^, and has been computed in Table 3.3 for various spectra using the differential 
cross section from Osterbrock (1974). The coUisional ionization rate is from Black (1981). 
The recombination rate is the total rate to all hydrogenic levels (Seaton 1959). 

Below we will see that in the ionized Stromgren bubbles that will appear around the 
galaxies or quasars, the photoionization rate is so much greater than the other rates that 
to a good approximation, equation (A.l) can be replaced by the following simple model for 
the IGM: 

• It is completely ionized {x = 1). 

• When a neutral hydrogen atom is formed through recombination, it is instantly pho- 
toionizcd again. 

Thus the only unknown parameter is the IGM temperature, which determines the re- 
combination rate, which in turn equals the photoionization rate and thus determines the 
rate of heating. 

Let us investigate when this model is valid. Near the perimeter of an ionized Stromgren 
sphere of radius r surrounding a galaxy, the number of UV photons per proton is roughly 

V 



1/3 



cn 



where Suv is the rate at which UV photons leave the galaxy. For an 05 star, the photon 
flux above the Lyman limit is approximately 3.1 x W^^s^^ (Spitzer 1968), so if each solar 
masses of baryons in a galaxy leads to production of a UV flux equivalent to that of an 05 
star, then 

77 > 0.77 ^„ t'^' (A.3) 

inside the sphere, where r\ = r/lMpc and Mq = M/10^Mq. When a fraction fg of all 
matter has formed galaxies of a typical total (baryonic and dark) mass M, then in the 
absence of strong clustering, the typical separation between two galaxies is 



R 



M ^ /15kpc\ / Me 



fsp) "^vi+JUvJ 



where Mq = M/IO^Mq. Thus r continues to increase until r k, R^ and spheres from 
neighboring galaxies begin to overlap. We are interested in the regime where z < 150. 
Substituting this and equation (A.3) into (A.2), we see that Xpi ^ Xd and Xpi ^ A^ec 
for any reasonable parameter values. Hence we can neglect coUisional ionization in equa- 
tion (A.l). Since Xpi 3> 1, the photoionization timescale is much shorter than the Hubble 
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timescale, so equation (A.l) will quickly approach a quasistatic equilibrium solution where 
the recombination rate equals the photoionization rate, i.e. 

^ -. ^rec ^ -. 
X -L — ~7 -L. 

In conclusion, the simple x = 1 model is valid for all parameter values in our regime of 
interest. 

When a hydrogen atom gets ionized, the photoelectron acquires an average kinetic energy 
of IkT*, where T* is defined by ^kT* = (Euv) - 13.6eV, and is the average energy of 

the ionizing photons (see Table 3.3). 

Since the timescale for Coulomb collisions is much shorter than any other timescales 
involved, the electrons and protons rapidly thermalize, and we can always assume that 
their velocity distribution is Maxwellian, corresponding to some well-defined temperature 
T. Thus shortly after the hydrogen gets photoionized, after the electrons have transferred 
half of their energy to the protons, the plasma temperature is T = ^T*. 

The net effect of a recombination and subsequent photoionization is to remove the ki- 
netic energy of the captured electron, say ^kTrjreciT), from the gas and replace it with 
|/cT*, the kinetic energy of the new photoelectron. Since the recombination cross section is 
approximately proportional to u^^, slower electrons are more likely to get captured. Hence 
the mean energy of the captured electrons is slightly lower than ffeT, i.e. r]rec{T) is slightly 
less than unity (Osterbrock 1974). We compute rjreciT) using Seaton (1959). The complica- 
tion that r]rec{T) ^ 1 turns out to be of only marginal importance: rjreci^O^^) ^ 0.8, which 
only raises the equilibrium temperatures calculated below by a few percent. 

The higher the recombination rate, the faster this effect will tend to push the temperature 
up towards T* . The two dominant cooling effects are Compton drag against the microwave 
background photons and cooling due to the adiabatic expansion of the universe. Line cooling 
from collisional excitations can be neglected, since the neutral fraction 1 — a; 0. Combining 
these effects, we obtain the evolution equation for the IGM inside of a Stromgren bubble: 



d{-z) 1 + z v^r+Tv 

where 



iT)T] 



(A.4) 



is (1 + z)~^ times the Compton cooling rate per Hubble time and T^-i,,. = T,.i,r q{1 + z). We 
have taken T^brfl ~ 2.726K (Mather et al. 1994). The factor of ^ in front of the Xrec term is 
due to the fact that the photoelectrons share their acquired energy with the protons. The 
average energy of the ionizing photons is given by the spectrum -P(z^) as (Euv) = h{v), where 



!^v-^P{v)a{v)du 

Here the photoionization cross section is given by Osterbrock (1974). Note that, in 

contrast to certain nebula calculations where all photons get absorbed sooner or later, the 
spectrum should be weighted by the photoionization cross section. This is because most 
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photons never get absorbed in the Stromgren regions (only in the transition layer), and 
all that is relevant is the energy distribution of those photons that do. -P(z^) is the energy 
distribution (W/Hz), not the number distribution which is proportional to P{u)/u. 

The spectral parameters (Euv) and T* are given in Table 3.3 for some selected spectra. A 
Planck spectrum -P(i^) oc u""'/ (^e^^^^'^ — gives quite a good prediction of T* for stars with 
surface temperatures below 30, OOOK. For very hot stars, more realistic spectra (Vacca, 1993) 
have a sharp break at the Lyman limit, and fall off much slower above it, thus giving higher 
values of T* . As seen in Table 3.3, an extremely metal poor star of surface temperature 
50, OOOK gives roughly the same T* as QSO radiation. The only stars that are likely to be 
relevant to early photoionization scenarios arc hot and short-lived ones, since the universe 
is only about 10^ years old at z = 100, and fainter stars would be unable to inject enough 
energy in so short a time. Conceivably, less massive stars could play a dominant role later 
on, thus lowering T*. However, since they radiate such a small fraction of their energy 
above the Lyman limit, very large numbers would be needed, which could be difficult to 
reconcile with the absence of observations of Population III stars today. If black holes are 
the dominant UV source, the stellar spectra of Table 3.3 are obviously irrelevant. A power 
law spectrum P(i') oc with a = 1 fits observed QSO spectra rather well in the vicinity 
of the Lyman limit (Cheney Sz Rowan-Robinson 1981; O'Brien et al. 1988), and is also 
consistent with the standard model for black hole accretion. 

Numerical solutions to equation (A.4) are shown in Figure 3.9, and it is seen that the 
temperature evolution separates into three distinct phases. In the first phase, the IGM is 
outside of the Stromgren regions, unexposed to UV radiation, and remains cold and neutral. 
In the second phase, the IGM suddenly becomes ionized, and its temperature instantly rises 
to ^T*. After this, Compton cooling rapidly reduces the temperature to a quasi-equilibrium 
value of a few thousand K. After this, in the third phase, T changes only quite slowly, and 
is approximately given by setting the expression in square brackets in equation (A.4) equal 
to zero. This quasi-cquilibrium temperature is typically many times lower than T* , since 
Compton cooling is so efficient at the high redshifts involved. 

A. 2 The Expansion of Stromgren Regions 

This rapid approach to quasi-equilibrium, where the IGM "loses its memory" of how long 
ago it became part of a Stromgren region, enables us to construct a very simple model for 
the ionization history of the universe. At redshift z, a volume fraction ^(z) of the universe 
is completely ionized and typically has a temperature of a few thousand K. The ionized 
part need not consist of non-overlapping spheres; it can have any topology whatsoever. The 
remainder is cold and neutral. 

Between the ionized and neutral regions is a relatively thin transition layer, where the 
IGM becomes photoionizcd and its temperature adjusts to the quasistatic value as in Fig- 
ure 3.9. As this IGM becomes part of the hot and ionized volume, the transition layer 
moves, and x{^) increases"^. 

^We are tacitly assuming that the UV luminosity of the galaxy that creates each Stromgren sphere never 
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As long as X < 1, all UV photons produced are absorbed instantly to a good approxi- 
mation. Thus the rate at which UV photons are released is the sum of the rate at which 
they are used to counterbalance recombinations inside the hot bubbles and the rate at which 
they are used to break new ground, to increase %. Thus 

Suvpp - a [1 )nx + , 

where a^'^\T) is the total recombination rate to all hydrogenic levels except the ground 
state^. Changing the independent variable to redshift and using equation (3.6), we find that 
this becomes 



dx I ■> (2) 1 + 2; _ _ 2 f fuvpp 



d{-z) ' "-^'^i + Wz v^Vi + ^. 



l + z 

1 ~\~ ^vii 



(A.5) 



Here 



\fl K, O.rUhnbT^ [l.04 - 0.51nr4 + 0.19T4 



1/2- 



is Hq^{1+z)~^ times the total recombination rate per baryon to all hydrogenic levels except 
the ground state. The fit is to the data of Spitzer (1968) and is accurate to within 2% for 
30K < T < 64, OOOK. Arel is to be evaluated at the quasi-equilibrium temperature T{z) 
discussed above. 

Using the values in Table 3.3 for the pessimistic, middle-of-the-road and optimistic esti- 
mates, the parameter fuvpp equals roughly 4, 190 and 24,000, respectively. 

In the absence of photon waste through recombination, equation equation (A.5) would 
have the solution x*{z) = fuvppfs{z), so the ionization efficiency is 

fioniz) = X{z)/X*iz). 

Since equation (A.5) is linear in x and the initial data is % = at some redshift, it is readily 
seen that the solution x{z) is proportional to fuvpp, the constant in front of the source term. 
Combining these last two observations, we see that fion is independent of fuvpp and hence 
independent of the poorly known parameters fbumi fuv and f^sc- 

Plots of fion{z) from numerical solutions of equation (A.5) are shown in Figure 3.10 
for various parameter values, and it is seen that the dependence on z is generally quite 
weak. Let us make use of this fact by substituting the Ansatz xiz) = fion{z) fuvpp fs{z) 
into equation (A.5), and setting floniz) ~ 0. Using equation (3.6) and an asymptotic 



decreases. Although obviously untrue, this is in fact an excellent approximation, since these early dwarf- 
galaxies correspond to perturbations far out in the Gaussian tail. Since fs{z) grows so dramatically as the 
redshift decreases and we move from five sigma to four sigma to three sigma, etc., almost all galaxies in 
existence at a given redshift are in fact very young, so that older ones that have begun to dim can be safely 
neglected. 

■^A rcionization directly to the ground state produces a UV photon that usually propagates uninterrupted 
through the highly ionized Stromgren region, and then ionizes another atom in the transition layer between 
the expanding Stromgren region and its cold and neutral surrounding. Thus recombinations directly to the 
ground state were included in the above calculation of the quasi-equilibrium temperature of the Stromgren 
bubbles, since the resulting UV photons could be considered lost from the latter. Here, on the contrary, 
recombinations directly to the ground state should not be included, since the UV photons they produce are 
not wasted from an energetics point of view. 
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approximation for the error function, we obtain 



1 + 0.48AreHl + Zyiry/^/T+n^ 

independent of fuvpp, which agrees to within 10% with the numerical solutions for all rea- 
sonable parameter values. This expression highlights the connection between /jon and the 
thermal evolution of the Stromgren bubbles: Essentially, the higher the quasi-static tem- 

(2) 

perature, the lower the recombination rate Xi-ec, and the higher fion becomes. 

The value of fion relevant to computing the ionization redshift is obviously that where 
-2 = Zion- As we have seen, Zion typically lies between 2zmr and Sz^r- Substituting T ^ 
2, 500K into the expression for A^ec and taking J7o ~ 1 and z = Zion ~ ^-^Zmr, the above 
reduces to 

f 1 

•'"^^ l + 0.8/lf^6(l + Z„,)3/2' 

SO we see that fion will be of order unity unless z^ir ^ 15 or Mlb 3> 0.03. 



Appendix B 

Comparing Goodness-of-fit 
Parameters 

In this appendix, we compare the performances of the probabihty product and the hkeUhood 
product as goodness-of-fit parameters. The notation is that of Chapter 7. 

First of all, what do we mean by a goodness-of-fit parameter rj being good? That rj leads 
to the correct model being ruled out at some confidence level x cannot be held against it 
— by definition of significance level, this happens a fraction (1 — x) of the time. Rather, 
the conventional criterion for rating goodncss-of-fit parameters is rejection power: given 
a model and a set of observations, one rj is said to be more powerful than another if it 
rejects the model at a higher level of significance. An example of a very stupid goodness- 
of-fit parameter, which in fact has the lowest rejection power possible, is a random variable 
77 drawn from a uniform distribution on [0, 1], thus containing no information whatsoever 
about the model or the observed data. Use of this parameter will reject the model at 95% 
confidence only 5% of the time, even if the data is blatantly inconsistent with the model. 

B.l Both rji and rjp can be "fooled"... 

Given a random variable Sj with probability distribution /j, we define the corresponding 
random variable for likelihood by Lj = fi{si) / fmax-, where wc chose the normalization con- 
stant frnax = DiaXx- fi{x) SO that we will always have < < 1. Thus for A'' experiments, 
the likelihood product 

N 

will also be a random variable on the interval [0, 1]. 

It is easy to construct examples where either the probability product rjp (as defined in 
Section 7.2.1) or the likelihood product r// give very low rejection power. The Achilles' 
heel of the probability product is multimodal distributions, where values near the mean are 
rather unlikely. For example, suppose N = 1 and we have a double-humped distribution 
such as 

/i(s) oc s^e-^'. 

If we observe a value si = 0, then rji would reject the model with 100% confidence whereas the 
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probability product fails miserably, rejecting with 0% confidence since Si equals the mean. 
The likelihood product, on the other hand, has the weakness that the highest likelihood 
may be attained far out in the tail of the distribution. Suppose for instance that N = 5 and 
we have triangle distributions 

fi{s) = 2s 

on the interval < s < 1. If all five observed values Si lie between 0.999 and 1, something 
is clearly wrong with the model. Since > 0.004^, using equation (7.1) indeed rejects 
the model at a confidence level exceeding 1 - ^^(0.004^) k, 99.999997%. Yet the likelihood 
product is near its maximum value: -qi > 0.999^ 0.995, which gives virtually no rejection 
power at all. 

B.2 ...but they usually give similar results. 

It is important to note that neither of the two examples above were particularly physical. 
The random variables arising from cosmic variance have Gaussian or chi-squared distribu- 
tions, and the same tends to holds for the various experimental noise distributions. Thus the 
probability distributions to which our goodness-of-fit parameter is applied in this chapter 
are unimodal (which eliminates the first example) and taper off to zero smoothly (which 
eliminates the second example). Hence for cosmological applications, goodness-of-fit pa- 
rameters should not be rated by their performance with such pathological distributions, but 
rather by their rejection power when applied to continuous, unimodal distributions. We will 
now compare the performance of r]p and r/; for a few such cases. 
For a symmetric exponential distribution 

fi{s) = \e-\^\ 

it is easy to see that the likelihood Li has a uniform distribution. This means that rjp and 
r/; will have identical distributions, for arbitrary A'^. It is straightforward to show that the 
same holds for symmetric triangle distributions 

h{s) = i-\s\. 

A third case where rjp and rji give identical results is when N = 1 and / is any smooth 
unimodal function. 

For a Gaussian distribution 

/,(.) = (2vr)-V2e-V2^ 

we have 

P(Lj<x) = 2erfc[(-21nx)^/2 . 

Although the probability distribution of r]i appears not to be expressible in terms of elemen- 
tary functions for arbitrary n, it is easy to show that r^i has a uniform distribution for the 
special case N = 2. Thus the likelihood product gives rejection at a confidence level of 

Pirn < m) = e-(^?+^i)/^ 
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Comparing this with the corresponding confidence level based on rjp shows a remarkable 
agreement between the two methods. Within the disc sf + < 4, over which P{ri < -fj) 
varies with many orders of magnitude, the two methods never differ by more than a factor 
of two. Thus, at worst, one may yield a confidence level of say 99.98% where the other yields 
99.99%. The probability product is stronger in slightly more than half of the (si, S2)-plane, 
roughly for regions that are more than 20° away from any of the coordinate axes. 

In conclusion, we have seen that for unimodal, continuous probability distributions, the 
likelihood product and the probability product tend to give fairly similar — in a few special 
cases even identical — results. Thus chosing one parameter over the other is more a matter of 
personal preference than something that is likely to seriously affect any scientific conclusions. 
There is however one important practical consideration: The probability distribution of rji 
depends on the probability distributions of the random variables Sj. This means that, apart 
from a few fortuitous special cases such as described above, it can generally not be calculated 
analytically. Rather, it must be computed numerically, through numerical convolution or 
Monte Carlo simulation. The probability distribution for rjp, on the other hand, is always 
known analytically, as given by equation (7.1), so the probability product is considerably 
simpler to use. 
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Appendix C 

Window Functions 



All the notation in this appendix is defined in Chapter 7. The results of CMB anisotropy 
experiments can be conveniently described by expanding the temperature fluctuation in 
spherical harmonics: 

1=2 m=-l 

(The monopole and dipole anisotropics have been removed from the above expression, since 
they are unmeasurable.) If the fluctuations are Gaussian, then each coefficient is an 
independent Gaussian random variable with zero mean (Bond &; Efstathiou 1987). The 
statistical properties of the fluctuations are then completely specified by the variances of 
these quantities 



(The fact that the variances are independent of m is an immediate consequence of spherical 
symmetry.) Different CMB experiments are sensitive to different linear combinations of the 
Q's: 

oo 

S = J2FiCi, (C.l) 

1=2 

where S is the ensemble-averaged mean-square signal in a particular experiment, and the 
"filter function" Fi specifies the sensitivity of the experiment on different angular scales. 
The filter functions for COBE and SP91 are 

j^icobe) ^ {21+1) ^^2(;+|)2 

pisp) =4e-?('+i)^EL=-.i^o^(«m), 

where Hq is a Struve function. Cc = 4.25° and 0"^ = 0.70° are the r.m.s. beamwidths for 
the two experiments, and a = 1.5° is the amplitude of the beam chop (Bond et al. 1991; 
Dodelson & Jubas 1993; White et al. 1993). 

For Sachs- Wolfe fluctuations in a spatially flat Universe with the standard ionization 
history, the angular power spectrum Ci is related to the power spectrum of the matter 
fluctuations in the following way (Peebles 1984; Bond &: Efstathiou 1987): 



Ci = — dkP{k)jt{k). (C.3) 
tttq Jo 
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Here tq is the conformal time at the present epoch, and 

jl{k) = J ji{kr)ViT)dT, 

where ji is a spherical Bessel function. The visibihty function V is the probabihty distri- 
bution for the conformal time at which a random CMB photon was last scattered. ji{k) 
is therefore the average of jiikr) over the last scattering surface. We have used the V of 
Padmanabhan (1993). 

We can combine equations (C.l), (C.2), and (C.3) to get the window functions for the 
two experiments: 

Wcobe = ^ E£2 j7(^)e-"'^'"'^^'(2^ + 1) 



The mean-square bulk flow inside of a sphere of radius a is (see, e.g., Kolb & Turner 
1990) 

2\ _ f^j 18 jKka) 



However, we must make two corrections to this result before applying it to the LP data. 
This formula applies to a measurement of the bulk flow within a sphere with an infinitely 
sharp boundary. In reality, errors in measuring distances cause the boundary of the spherical 
region to be somewhat fuzzy. If we assume that distance measurements are subject to a 
fractional error e, then the window function must be multiplied by e^'^^'^")^. We have taken 
e = 0.16, the average value quoted by LP. It should be noted that this value varies from 
galaxy to galaxy in the LP sample, due to the distance estimation technique used, and that 
a more accurate window function that reflects the discrete locations of the Abell clusters 
used in the survey should take this into account. 

The second correction has to do with the behavior of the window function at small 
k. Equation (C.4) applies to the velocity relative to the rest frame of the Universe. The 
velocity measured by LP is with respect to the CMB rest frame. If there is an intrinsic 
CMB dipole anisotropy, then these two reference frames differ. Therefore, we must include 
in equation (C.l), a term corresponding to the intrinsic CMB dipole. This correction was 
first noticed by Gorski (1991). After applying both of these corrections, the LP window 
function is 



^ tt'^Tq \ ka kro 
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